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The radioactive isotope La, known to be produced by 
deuteron and slow neutron bombardments of lanthanum 
is found to have a half-life of 40.0+0.3 hours. This isotope 
has now been produced by the reaction Ce™(n, p)La™. 
Evidence for the occurrence of a d, y reaction has been 
obtained through the formation of La’ by the reaction 
Ba™8(d, y)La™. The 40-hour isotope decays with the 
emission of 1.41++0.05-Mev electrons and 2.00+0.05-Mev 
gamma-rays. Deuteron and also proton bombardments of 


barium produce a new radioactive isotope of lanthanum 
of half-life 17.5+0.5 hours. This activity decays by the 
process of K-electron capture with the emission of x-rays 
identified as the characteristic K radiation of barium. 
The presence of a 0.88+0.1-Mev gamma-ray of low in- 
tensity indicates that the resulting barium nucleus may 
also be left in an excited state. This 17.5-hour activity is 
assigned to La'?, 





INTRODUCTION 


RTIFICIAL radioactivity in lanthanum was 
first induced by a slow neutron activation 
of the oxide, from which a 45.6+4.8-hour ac- 
tivity, assigned to La™®, was obtained.' The 
determinations of the half-life of the lanthanum 
activity produced by different nuclear processes 
and described by several investigators vary to a 
rather large extent, namely, from 31 hours to 46 
hours.?~® 
It is the purpose of this paper to investigate 
more accurately the La‘ period and in addition 
to describe a new 17.5-hour period found in 
lanthanum after activation of barium with 
deuterons and protons. 


1 J. K. Marsh and S. Sugden, Nature 136, 102 (1935). 

2M. L. Pool and L. L. Quill, Phys. Rev. 53, 437 (1938). 

*Q. Hahn and F. Strassmann, Naturwiss. 27, 11 (1939). 

‘OQ. Hahn and F. Strassmann, Naturwiss. 28, 54 (1940). 

®’W. G. Chlopin, M. A. Passwik-Chlopin, and N. F. 
Wolkov, Nature 144, 595 (1936). 

* G, N. Glasoe and J. Steigman, Phys. Rev. 58, 1 (1940). 
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EXPERIMENTAL 


All bombardments were carried out in the 
42-inch cyclotron of The Ohio State University. 
The bombarding particles used were approxi- 
mately 20-Mev alpha-particles, 10-Mev deu- 
terons, 5-Mev protons, slow neutrons produced 
by the Be+d reaction, and fast neutrons from 
the Li+d reaction. Chemical separations were 
made after all activations except those with slow 
neutrons. 

In all bombardments of lanthanum and cerium, 
Hilger H.S. brand oxides were used while, in the 
barium activations, Hilger ‘“‘Specpure”’ barium 
chloride was used. A spectroscopic analysis ac- 
companying the material revealed the following 
maximum impurities: 

Lanthanum oxide, 99.9 percent, Lab No. 6404 

—very slight traces of Na, Mg, Ca, Bi, Sm; 
Gd, 0.015 percent; Nd, 0.1 percent. 

Cerium oxide, 99.99 percent, Lab No. 6401— 

traces of Mg and Gd. 
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Barium chloride, 99.98 percent, Lab No. 8442 
—Ca, 0.0015 percent; Fe, 0.005 percent; 
Mg, less than 0.001 percent; Na, about 0.001 
percent; Pb, 0.01 percent; and Sr, absent. 
The other source of barium used was a special 

grade of the nitrate prepared by a modification of 
the method used in preparing barium for atomic 
weight determinations. Spectroscopic analysis of 
this material revealed that it was of a high degree 
of purity. Hilger cesium chloride was also bom- 
barded in the course of the work. 

Measurements of the radioactivity were made 
on a Wulf unifilar electrometer with an ionization 
chamber containing two atmospheres of freon, or, 
in some cases, on another Wulf unifilar electrome- 
ter with an ionization chamber having a 0.25-mil 
aluminum window and containing air at atmos- 
pheric pressure. The latter chamber was used in 
the detection of soft beta-radiation. 

For the study of the emitted radiations, a 
Wilson cloud chamber of 15-cm diameter and a 
180° focusing type magnetic spectrometer with a 
radius of curvature of 16 cm were available. 


THE 40-HOUR PERIOD 


A very intense 40.0-hour lanthanum activity 
emitting electrons and gamma-rays, was pro- 
duced by the deuteron bombardment of the 
lanthanum oxide. The procedure used for the 
purification of lanthanum after activation was 





Doys after Bombardment 


Fic. 1. Decay curve for lanthanum fraction from La+d. 
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identical with that described for praseodymium.’ 
Since the decay curve of the deuteron-activated 
lanthanum, shown in Fig. 1, was linear through 
11 half-lives, a very reliable evaluation of the 
period could be made. From an average of three 
such determinations, the half-life was found to be 
40.0+0.3 hours. 

Activation of lanthanum with slow neutrons 
gave a value of 39.8 hours for the period of La™®, 
Although the activity produced in this case was 
weak compared to that produced by the deuteron 
bombardment, it yet was easily followed through 
several half-lives. 
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Fic. 2. Gamma-ray absorption in lead of 40-hour 
lanthanum. 
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RADIOACTIVE ISOTOPES OF 


After-cerium was bombarded with fast neutrons 
a lanthanum fraction was separated from the 
cerium. The activity of the lanthanum fraction 
also decayed with a 40-hour period. This activity 
was likewise weak compared to that obtained by 
deuteron bombardment. The new reaction, 
Ce'°(n, p)La'° seemed therefore to be established. 

An energy of 2.0 Mev for the gamma-rays 
emitted by the 40-hour period has been de- 
termined by absorption in lead. Measurements, 
as shown in Fig. 2, made with a very strong 
source gave a single absorption coefficient of 
0.507 +0.005 cm~, corresponding to an energy of 
2.00+0.05 Mev. 

From aluminum absorption measurements, 
shown in Fig. 3, and from a beta-ray spectrum, 
shown in Fig. 4, taken with the magnetic 
spectrometer, the maximum energy of the beta- 
particles was estimated to be 1.41+0.05 Mev. A 
tail of high energy Compton recoil electrons pro- 
duced by the 2.0-Mev gamma-rays was observed 
in both cases. 

In Fig. 5, a section of the periodic table is 
reproduced showing the assignments of and the 
reactions producing the various radioactive 
isotopes in the lanthanum region. 


THE 17.5-HOUR PERIOD 


A new radioactive isotope of lanthanum with a 
half-life of 17.5 hours has been found to be pro- 
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Fic. 3. Beta-ray absorption in aluminum of 40-hour 


lanthanum. 


LANTHANUM 69 


duced by bombardments of barium with 10-Mev 
deuterons. From the lanthanum fractions of these 
bombardments, complex decay curves were ob- 
tained; a portion of one is shown in Fig. 6. The 
longer period was observed to be 40 hours. Upon 
subtraction of this activity from the original 
curve, a period of 17.5+0.5 hours resulted. 
Proton bombardment of barium gave only a 
17.5-hour activity in the lanthanum fraction. 
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Fic. 4. Beta-ray spectrum of 40-hour lanthanum. 


In order to observe the 17.5-hour period a 
careful separation of barium and lanthanum is 
necessary because of the very strong beta-radia- 
tion emitted by Ba"® of 85.6-minute period. It has 
been found that even 0.01 percent of this activity 
immediately after bombardment, if left in the 
lanthanum fraction, is sufficient to mask the 
correct evaluation of the 17.5-hour period. 

The chemical procedures were arranged in the 
following manner: Bombarded barium chloride, 
which was contained in a platinum holder, was 
dissolved in nitric acid with the addition of 
lanthanum oxide. The radioactive chlorine was 
mainly removed by evaporation and the nitrate 
solution was.converted into a chloride solution. 
Added platinic acid was reduced and platinum 
was filtered off. 

Separation of barium and lanthanum was first 
made with ammonium hydroxide. Precipitated 
lanthanum was dissolved and reprecipitated as 
an oxalate in presence of stabilizers for active 
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impurities other than those of the rare-earth 
group. Separate purification of lanthanum was 
performed from phosphoric acid to remove 
radioactive phosphorus, which was present be- 
cause of the Cl*5(n, a) P* reaction. 


K-Electron Capture 


The radiation from the 17.5-hour period con- 
sisted entirely of x-rays and gamma-rays as 
evidenced in the aluminum absorption curve of 
Fig. 7 and in the copper absorption curve of 
Fig. 8. Because of the presence of a low intensity 
gamma-ray, a gradual bending of these curves 
was observed. Upon subtraction of the activity of 
the gamma-ray in each case, new curves were 
obtained, from which were determined mass 
absorption coefficients of 0.86 cm?/g in aluminum 
and of 9.11 cm?/g in copper. These coefficients 
corresponded, respectively, to x-ray wave-lengths 
of 0.367A and 0.389A, or to an average wave- 
length of 0.378+0.011A. The close agreement of 
this wave-length with the estimated mean wave- 
length, 0.376A, of the K series for barium, sug- 
gested at once that the radioactive lanthanum 
might be decaying by the process of K-electron 
capture. 

The deuteron bombardment of barium could, 
however, give rise to a radioactive isotope of the 
stable La™®, which would then decay to the 
ground state by the emission of a gamma-ray. A 
possibility therefore exists that the x-rays ob- 
served might arise from the internal conversion 


Fic. 5. Isotopic chart of the lanthanum region. 












of this gamma-ray. As the K-radiation produced 
in this case would be that of lanthanum while in 
the case of K capture it would be that of barium, 
a more refined measurement of the wave-length 
of these x-rays should make it possible to 
distinguish between the two processes. 








Critical Absorption 





The most widely used method for the determi- 
nation of the wave-length of the x-rays emitted 
by radioactive substances is provided by critical 
absorption. From the data in Table I, it may be 
seen that tellurium should absorb the x-radiation 
from both barium and lanthanum strongly, but 
that iodine should absorb only about 20 percent 
of the barium x-rays as compared with 74 percent 
of the lanthanum x-rays. In addition, if the 
K radiation were that of lanthanum, the iodine 
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Fic. 6. Decay curve for lanthanum fraction from Ba+d. 
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absorption curve should reveal the presence of the 
two principal lines of the K series, one of which 
(Ka;) would be absorbed and the other (Kaz), 
transmitted. 

Filters of powdered tellurium metal and of 
sodium iodide were therefore prepared for the 
purpose of determining the mass absorption 
coefficients of the x-rays in the two substances. 
The main source of error in measurements of this 
kind lies in the difficulty in preparing thin filters 
of uniform thickness.* An estimate of this error, 
for an extreme case, was made by measuring the 
absorption coefficients for the x-rays emitted by 
the 17.5-hour isotope with tin foils of 9.8 mg/cm? 
and with filters of 30-mesh tin powder. The 
coefficient obtained with the foils was 28.9 cm?/g, 
which is very nearly the true value, while with 
the powder, it was found to be only 12.6 cm?/g. 

The x-rays from the 17.5-hour activity had a 
mass absorption coefficient, in tellurium, of 18 
cm?/g. A correction, which was estimated from 
the tin absorption measurements just described, 
was made necessary by the fact that the thin 
filters of tellurium metal used were not homo- 
geneous. The value of the absorption coefficient 
in tellurium then obtained was approximately 30 
cm?/g, showing that the x-rays were strongly 
absorbed by the metal. This strong absorption 
was to be expected for either barium or lanthanum 
x-rays. 

In iodine, a single mass absorption coefficient 


TABLE I. Data on K series lines of barium and lanthanum 








K series line Kaz Kai Kp KB: 





Fraction of total 

intensity 0.26 0.53 0.16 0.05 
Ba wave-length (A) 0.38899 0.38443 0.34022 0.33222 
La wave-length (A) 0.37466 0.37004 0.32726 0.31966 


Critical absorption limits: Te, 0.38926A; I, 0.37244A. 











of 8.5 cm?/g was found for the x-rays from the 
17.5-hour period. Correction for inhomogeneity 
was relatively small in this case. Barium x-rays 
would give a simple absorption curve in iodine, 
corresponding to a mass absorption coefficient 
estimated to be approximately 11.7 cm?/g. The 
observed value, 8.5 cm?/g, is in satisfactory 


* A. H. Compton and S. K. Allison, X-Rays in Theory 
and Experiment (D. Van Nostrand Company, Inc., New 
York, 1935), second edition, p. 517. 
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Fic. 7. Aluminum absorption curve for 17.5-hour 
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agreement with this estimated coefficient. As the 
iodine was relatively transparent to the x-rays, 
the possibility that the K radiation was that of 
lanthanum seemed definitely excluded. 

The evidence from these critical absorption 
measurements, therefore, leads to the conclusion 
that the x-rays are those of barium and that, 
consequently the 17.5-hour isotope decays by the 
process of K-electron capture. 


The 0.88-Mev Gamma-Ray 


The gamma-ray observed in the aluminum and 
copper absorption measurements was also found 
to decay with a half-life of 17.5 hours. From lead 
absorption measurements, the energy of this 
penetrating radiation was determined to be 
0.88+0.1 Mev. Should this gamma-ray be in- 
ternally converted, electrons of approximately 
0.84 Mev would be produced. As no such electrons 
were observed, the possibility that the x-rays 
might arise from the internal conversion of this 
gamma-ray has been excluded independently of 
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the indentification of the K x-rays as those of 
barium. 

No positrons were observed in any of the cloud- 
chamber photographs taken of the 17.5-hour 
activity. This would indicate that K-electron 
capture is a very much more probable type of 
disintegration for this isotope than positron 
emission. The absence of annihilation radiation 
also confirmed the failure to observe positrons. 

The number of x-rays relative to the number of 
gamma-rays is about 50 which means that in 
about two percent of the capture processes the 
barium nucleus is left in an excited state. 


Assignment 


The assignment of the 17.5-hour period has 
been made to La'*’ as illustrated in Fig. 5. 
Bombardments of cesium with alpha-particles 
and of lanthanum with fast neutrons failed to 
produce any activity in the lanthanum fractions 
other than a small amount of the 40-hour activity 
in the latter case. It seems therefore unlikely that 
the 17.5-hour activity can be due to either La'* 
or La'"*®, It is improbable that the lighter and less 
abundant isotopes of barium would produce an 
activity as intense as was obtained. In view of 
these considerations, the activity was assigned to 
La'*’, which is probably produced by the reactions 
Ba'*6(d, n)La'*? and Ba"7(p, n)La'’. 


EVIDENCE FOR A d, y REACTION 


A weak 40-hour period was observed in the 
lanthanum fractions from the deuteron bom- 
bardments of barium. The decay curve for one 
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such fraction, shown in Fig. 6, has already been 
considered. It was found that the ratio of the 
intensity of the total radiation for this 40-hour 
activity to that of the gamma-radiation alone 
was of the same order of magnitude as that ob- 
tained previously for the 40-hour activity of La™, 

The activity of the 17.5-hour lanthanum 
calculated for an infinite bombardment, relative 
to that calculated for this 40-hour period re- 
mained practically constant in each of the several 
deuteron activations of barium. As the spectro- 
scopic analysis for both sources of barium used in 
the bombardments showed the complete absence 
of lanthanum, there was no possibility of pro- 
ducing the 40-hour activity by a d, p reaction 
from a lanthanum impurity. It therefore 
seems reasonably certain that the reaction 
Ba'8(d, y)La™® does take place. 

There are only a few cases where an activity 
can be produced unequivocally by a d, y reaction. 
In general, the distribution of the stable isotopes 
is such that the more probable d,m or d, 2n 
reactions will occur with a consequent masking 
of any contribution made by a d, y reaction. 

The probability that the d,y reaction will J 
occur in this particular case, relative to the 
occurrence of the d, m reaction was calculated to 
be 1 to 1300. 

It is a pleasure to acknowledge the support 
received from Mr. Julius F. Stone and from The 
Ohio State University Development Fund. Ac- 
knowledgment is also made of the valuable 
assistance obtained through W.P.A. cyclotron 
project No. 65-1-42-89. 
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A Simple Method for the Investigation of Secondary Electrons Excited by y-Rays 
and the Interference of These Electrons with Measurements of Primary 


6-Ray Spectra 


LisE MEITNER 


Forskningsinstitutet for Fysik, Stockholm, Sweden 


(Received December 17, 1942) 


It has been shown that the special directional distribution of the Compton electrons can 
be used to obtain their absorption curve and therefore their upper limit of energy without 
employing coincidence measurements. This fact has to be taken into account when inspecting 
primary §-spectra by absorption measurements. It has been found that the complexity of 
the primary 8-ray spectrum of Fe** reported by Livingood and Seaborg is apparently due to 
the interference of secondary electrons. Some problems connected with the directional distri- 
bution of the Compton electrons and concerning the efficiency of y-counters are discussed. 








N a paper! dealing with the 8- and y-rays 
emitted by Sc*® we had found that under 
certain experimental conditions the secondary 
electrons excited by the y-rays may be observed 
quite easily without using coincidence arrange- 
ment and it was pointed out that this fact has 
to be taken into account, when inspecting 
primary §-ray spectra by absorption measure- 
ments. It was also briefly mentioned that, in 
case only light elements such as Al are used, 
the presence of the secondary electrons is due 
to the special directional distribution of the 
Compton electrons, the greater part of these 
electrons moving at large angles relative to the 
direction of the y-rays by which they are 
produced. So if the radioactive source is placed 
at some distance from a thin-walled aluminum 
counter the Compton electrons excited near the 
source and reaching the counter may be due in 
a large part to y-rays which do not enter the 
counter at all. If the solid angle which the 
source subtends at the counter, is much smaller 
than z, a larger number of these electrons will 
enter into the counter than would correspond to 
the saturation number excited by the y-rays 
impinging on the counter. Hence if absorbing 
foils of Al are inserted directly in front of the 
counter, one can obtain the absorption curve of 
these secondary electrons and in this way 
determine the energy of the respective y-rays. 
Of course, one can produce a surplus of 
secondary electrons as well by placing a lead 
absorber in front of the source and measure the 


1 L. Meitner, Arkiv f. Mat. Astron. och Fysik, Stockholm 
28B, No. 14 (1942). 





absorption of the secondary electrons excited in 
the lead by inserting the absorbing Al foils in 
front of the counter. In this case one gets in 
addition to the surplus of Compton electrons an 
excess of photoelectrons as no comparable photo- 
effect takes place either in the wall of the counter 
or in the foils inserted in front of the counter. 

The correctness of these statements has been 
proved by a series of experiments, a detailed 
report of which is to be published in the issues 
of the Academy of Sciences of Stockholm. 

Here I shall mention only some of the main 
points. 

To prove that actually the directional distribu- 
tion of the Compton electrons can provide an 
excess of secondary electrons entering the counter 
the following experiment was carried out. A 
thick layer of metallic thorium emitting practi- 
cally only the 8- and y-rays belonging to Th C 
+C” was placed at 18-cm distance from a 
thin-walled aluminum counter. The source was 
covered with aluminum plates of 4-mm total 
thickness in order to cut off the primary §-rays 
and to get at the same time secondary electrons 
excited in the Al by the y-rays. Two sets of 
measurements were made: first the Al plates 
placed in front of the source were 1.23.2 cm, 
secondly the plates had an area of 5X10 cm. 

In the first case the activity found amounted 
to 184/min. When 4.1-mm Al (corresponding to 
the maximum range of the secondary electrons) 
was inserted directly in front of the counter, 
in order to absorb all the secondary electrons 
coming from the source, the activity obtained 
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was 146/min. In the second case the activity 
found with no absorber was 216/min., with 4.1- 
mm Al in front of the counter the rate was 
147/min. The equality of the two values found 
with 4.1-mm Al in front of the counter proves 
that the surplus of activity obtained with the 
larger Al plates in front of the source and no 
absorber in front of the counter is actually due 
to additional excited Compton electrons (as is 
to be expected from their directional distribution) 
and is not produced by scattered y-rays. 

As to the entering of the secondary electrons 
into the measurements of the primary electrons 
there are two cases to be distinguished. 

(1) The primary §-rays of the substance to be 
investigated have an upper limit of energy much 
lower than the maximum energy of the secondary 
electrons coupled with the emitted y-rays. If 
the source is a thin layer, one will find, with the 
absorbers inserted closely in front of the source, 
that beyond the end of the primary #-ray 
spectrum an increase of the activity takes place 
because of the increasing number of excited 
electrons with increasing thickness of the alumi- 
num absorbers. With the absorbers inserted in 
front of the counter the interference of the 
secondary electrons will be much smaller, be- 
cause the y-rays are not impinging at large angles 
upon the counter or upon the absorbing foils. 

If the source is a thick layer, the absorption 
curve obtained with the absorbing foils in front 
of the source will, beyond the range of the 
primary 8-rays, exhibit a constant value, as far, 
of course, as the absorption of the y-rays is 
negligible and as far as the total thickness of the 
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inserted absorbers is small in comparison with 
the distance between source and _ counter. 
Otherwise the surplus of Compton electrons 
reaching the counter decreases with increasing 
thickness of the absorbing foils, since more and 
more of these electrons are produced nearer to 
the counter. If the absorbing foils are placed in 
front of the counter, one gets in addition to the 
absorption curve of the primary #-rays the 
absorption curve of the secondary electrons. 

Therefore in the case of a thick layer the 
curve for the primary §8-rays is least disturbed 
by secondary electrons, if the absorbing foils 
are inserted as closely as possible in front of the 
source. At any rate if the two curves obtained 
with the absorbing foils in front of the source in 
comparison with those in front of the counter 
display differences in the ranges, one can con- 
clude that there are secondary electrons entering 
into the measurements. 

(2) The maximum energy of the secondary 
electrons is smaller than the upper limit of the 
primary 6-ray spectrum. In this case the second- 
ary electrons do not interfere with the determi- 
nation of the upper energy limit for the primary 
B-rays. By absorbing the latter in Al foils 
inserted in front of the source, one can obtain 
the absorption curve of the secondary electrons 
(excited in the foils close to the source) by 
putting absorbing foils in front of the counter. 

The following figures give some examples of 
the procedure described. 

Figure 1 relates to experiments carried out 
with Sc**. For the curves a and b we used a 
small and rather thin layer of Sc*® (as Sc,03) 
supported by an Al sheet of 0.1-mm thickness. 
The source was placed at 5.5-cm distance from 
the Al counter. To exclude as far as possible 
scattering from the surroundings the counter 
was mounted on a plate of ebonite and no other 
material but wood was employed in the neighbor- 
hood of source and counter. 

Curve a was obtained with the absorbing foils 
inserted in front of the source, curve 6 represents 
the measurements with the absorbers put in 
front of the counter. The absorbing foils were in 
both cases just sufficiently large to cover the 
source or the counter. 

Curve a shows, beyond the end of the primary 
B-ray spectrum (R=0.24 mm Al, we take into 
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INVESTIGATION OF SECONDARY ELECTRONS 


account the thickness of the counterwall), an 
increase of the intensity due to the formation 
of secondary electrons by the strong y-line of 
1.5-Mev energy. 

Curve b exhibits in addition to the steep slope 
corresponding to the absorption of the primary 
g-rays a slight indication of a much flatter slope 
caused by a small number of additional secondary 
electrons. This number of secondary electrons 
excited in the thin layer of Sc2,O; itself did by 
chance nearly correspond to the saturation 
number of Compton electrons excited in the 
wall of the counter by the impinging y-rays. It 
may be mentioned that without the interference 
of the secondary electrons curve a should begin 
with a steeper slope than curve }, because the 
scattering of the (primary) f-rays is more 
important in the former case than in the latter. 

From the difference between the terminal 
values of both curves one can compute the 
excess of Compton electrons produced in Al 
placed closely to the source over the number 
produced in the Al when inserted in front of the 
counter. It amounted to 57 percent. 

Curve c exemplifies the determination of the 
range of secondary electrons. The source was a 
sample of Sc** [as Sc2(OH); ] which filled a glass 
tube of 0.8-mm thickness, surrounded by 0.5-mm 
Al. It was placed at 9-cm distance from the 
counter and the absorbing foils were inserted 
in front of the counter. All the electrons are, of 
course, of secondary origin. From the curve c 
(the dotted line represents the absorption of the 
y-rays in Al) one deduces a range of 2.0 mm in 
Al, corresponding to a maximum energy of the 
electrons of 1.3 Mev, in very good agreement 
with the value to be expected from the y-energy 
of 1.5 Mev. 
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Fic. 2. Absorption curve of secondary electrons of the 
gamma-rays of ThC”. 
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Fic. 3. Absorption curves of primary and secondary 
electrons of Fe®. 


In order to establish that this simple method 
can be used as well when more than one y-line 
is present, experiments were carried out with 
the three y-lines of 0.537-, 0.592-, and 2.65-Mev 
energy emitted by ThC”. Here the primary 
8-particles have an upper limit of energy larger, 
in fact comparable with the maximum energy of 
the secondary electrons excited by the y-rays. 
The source, a thick sample of metallic Th was 
arranged at 18-cm distance from the counter, 
so that the rays were impinging nearly parallel 
upon the counter. The source was covered with 
0.75-mm Pb+4-mm Al. The photoelectrons 
excited in the lead are absorbed in the Al and 
only Compton electrons excited in the Al reach 
the counter. Lead was chosen in order to absorb 
the y-line of 0.24 Mev belonging to Th B and 
emerging at a small rate from the thick Th 
sample. 

The curve of Fig. 2 relates to these measure- 
ments. It shows clearly a form different from 
that of curve c of Fig. 1, exhibiting two groups 
of secondary electrons with ranges of about 
0.36-mm and of 4.2-mm Al. The energy corre- 
sponding to a range of 0.36 mm Al can only be 
estimated as lying between 0.4 and 0.45 Mev. 
The energy corresponding to a range of 4.2-mm 
Al is 2.40 Mev. The values to be expected from 
the y-energies of 0.592 and 2.65 Mev are 0.41 
and 2.42 Mev, respectively. The agreement of 
the measurements with these values is therefore 
rather good. 

Measurements with Pb in front of the source 
instead of Al gave practically the same results. 
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THE PRIMARY §-SPECTRUM OF Fe** 


By investigating the B- and y-rays of Fe®® we 
found an example where the complexity of the 
primary §-spectrum reported by Livingood and 
Seaborg? apparently has to be explained by the 
interference of secondary electrons. The authors 
found by inspecting the §-ray spectrum that 
most of the primary §-particles have a maximum 
range of 0.09 g/cm? Al while a small number 
extends to 0.35 g/cm? Al. In addition there is 
emitted a y-line of about 1-Mev energy. 

By carrying out experiments of the above 
described procedvre with Fe*® the curves of the 
Fig. 3 were obtained. 

The Fe®® was a thick layer of metallic iron. 
It was supported by a sheet of wood and placed 
at 8-cm distance from the counter. For the 
measurements given in curve c¢ of Fig. 3 the 
source was covered with 4-mm Pb, because 
there was present a weak y-radiation of about 
0.25 Mev which I could not identify. It may be 
mentioned that for y-rays of 1-Mev energy the 
energy of the photoelectrons excited in the K 
level of Pb is practically the same as the maxi- 
mum energy of the Compton electrons. The 
absorbing Al foils were, of course, inserted in 
this set of measurements in front of the counter. 

From the curve a, which was obtained with 
the absorbers in front of the source, one deduces 
for the primary f-particles a range of 0.36-mm 
Al corresponding to 0.098 g/cm?, a value some- 
what larger than that one found by Livingood 
and Seaborg. There is no indication of a second 
8-group. 

Curve b exhibits in addition to the absorption 
of the primary §-rays a flatter slope due to the 
presence of the secondary electrons, the range 
of which is 1.1-mm Al or 0.30 g/cm*. The same 
value is found in curve c, where, of course, only 
secondary electrons enter into the measurements. 
So the conclusion that Fe®® emits a single 
primary §-group is supported also by the 
measurements given in curve Cc. 

The fact that the surplus of Compton electrons 
to be expected from their directional distribution 
is found in the described experiments, proves 
that the original directional distribution is not 
much disturbed by the scattering of the electrons 
when passing through the aluminum foils. It is, 


- J. Livingood and G. T. Seaborg, Phys. Rev. 54, 51 
(1938). 
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however, possible to investigate this point in a 
different way. Since no Compton electrons are 
excited at directions forming angles larger than 
a/2 with the incident y-radiation, the back wall 
of the counter will have a noticeable influence 
only if it gives rise to an appreciable scattering 
of the Compton electrons. Hence the following 
experiment was carried out with y-rays of about 
1-Mev energy. 

The wall of the counter facing the source was 
covered with 2-mm Al in order to insure that 
no secondary electron coming from the source 
entered into the counter. The intensity of the 
radiation was measured first when the back 
wall of the counter was not screened by any 
foils, secondly, when different foils of Al were 
placed closely to the back wall. If the Compton 
electrons excited in the back wall are noticeably 
scattered in Al, a large increase of the intensity 
is to be expected with increasing thickness of the 
foils inserted behind the back wall of the counter, 

The values obtained with no Al, 0.1- and 0.5- 
mm Al, respectively, were Jo=52/min., Io, 
= 53.7/min., and Jo,;=56.6/min. The increase is 
very small proving that the scattering in Al is 
not important and does not change the original 
directional distribution at a large rate. 

If instead of the Al foils a Pb foil of 0.11-mm 
thickness corresponding to 0.4-mm Al was 
inserted closely to the back wall of the counter, 
the activity found amounted to 67/min. This 
large increase is due partly to the photoelectrons 
excited in the lead having originally an isotropic 
distribution, partly to the fact that the scatter- 
ing of the (Compton) electrons is much more 
important in lead than in aluminum. 

These results are of some significance when 
computing for different y-energies the efficiency 
of y-counters made of different materials. As 
long as a y-radiation does not noticeably excite 
photoelectrons in the counter material and as 
long as the scattering of the (Compton) electrons 
is not important, the back wall of the counter 
does not contribute to the activity at the same 
rate as the front wall, as it was supposed usually. 

The author wishes to express her sincere 
gratitude to Professor N. Bohr and Professor J. 
Jacobsen (Copenhagen) for the irradiation of 
the substances used in this work. 

Acknowledgment is also made for a grant 
given by the Academy of Sciences of Stockholm. 
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The charge of radium E §8-particles was measured by first measuring the charge carried 
per second by a magnetically analyzed beam and then, after the beam had been reduced by 
a known fraction by the decay of the source, by counting the number of particles per second 
in the beam. The average value obtained was —4.809+0.03;x 10-" e.s.u. in agreement with 
the accepted value of the charge of the electron. Also a half-life of 4.999 days for radium E 
was observed. No dependence of either the charge or the half-life upon Hp was detected. 





INTRODUCTION 


HE method of this experiment? is to measure 
the current carried by a beam of radium E 
8-particles emerging from a 8-ray spectrometer 
by means of a Faraday collector and a vacuum 
tube electrometer and then to replace the col- 
lector by a Geiger counter and to count the num- 
ber of particles emerging per second. This is a 
modification of the method which has been used 
to measure the charge of the alpha-particle.® 
Until recent developments in Geiger counter 
technique it has not been possible to determine 
the charge of the §-particle in this way, but 
instead the first experiment of this type on 
B-particles‘ was performed to measure the effi- 
ciency of the counter for 8-rays. 

Because it was necessary to have a current 
larger than 5X10-'® ampere (3X10* electronic 
charges per second) for accurate current measure- 
ments while it was not possible to measure accu- 
rately counting rates greater than 300 per sec., 
for the counting measurements it was necessary 
to reduce the strength of the beam by a known 
fraction by allowing the source to decay for 
nearly ten half-lives. The amount of the decay 
could not be calculated accurately over such a 
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long period of time from the known half-life of 
radium E because of impurities of radium D. 
Therefore the spectrometer was equipped with a 
second exit slit and set of baffles similar to the 
first ones through which the beam could be sent 
by reversing the magnetic field into an ionization 
chamber which could be connected to the vacuum 
tube electrometer. (See Fig. 1.) At the times of 
the collector current and counting measurements 
the corresponding ionization currents (Js and Jy, 
respectively) were observed, these being measures 
of the relative strengths of the source. Then if e 
is the charge of the 8-particle, Jy the collector 
current, and n’ the number of counts per second, 


Ip=n'e(Is/Iw) or e=(Ir/Is)/(n'/Iw) (1) 


(Ir, Is, Iw, and n’ are all measured at the same 
value of the magnetic field). 

The collector current and both ionization cur- 
rents were measured by means of an electrometer 
tube and the null method of Townsend,' in which 
the grid is maintained at zero potential by 
inducing upon it through a condenser a charge 
equal in magnitude but opposite in sign to that 
being measured. The current is given by the 
product of the condenser capacity C and the rate 
of change of potential V’ supplied to the con- 
denser such that the grid is maintained at zero 
potential. A condenser C, of small but very 
accurately known capacity was used for measure- 
ments of the small currents Jp and Jw, and a 
large condenser C; was used for the ‘‘strong”’ 
ionization current Js. Hence Eq. (1) becomes 


Ci Vp'/ V3’ 
mo( 2) 2 
C3 n'/ Vw’ 


5 J. S. Townsend, Phil Mag. 6, 603 (1903); H. G. S. 
Mosely, Proc. Roy. Soc. 87, 230 (1912). 


(2) 
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where the subscripts F, S, and W refer, as before, 
to the collector current and to the “‘strong”’ and 
“weak” ionization currents, respectively. It is to 
be noted that the absolute value of C; is not 
involved but only its ratio to C;. 


APPARATUS 


The double 8-ray spectrometer, Faraday col- 
lector, and the ionization chamber are illustrated 
in Fig. 1. The sides of the spectrometer were 
made of Armco iron blocks 18 cm X12 cm X 1 in. 
which became the actual pole pieces when the 
apparatus was placed in position. These were 
held 6 cm apart by a brass frame which contained 
the entrance slit (4X15 mm), the two exit slits 
(2X40 mm), and the aluminum baffle system. 
These were designed to give as large an intensity 
as possible and still restrict the emerging beam to 
within 30 degrees to the normal of the exit slit so 
that the beam could easily be trapped in either 
the collector or counter. The radius of the 
trajectory was about 3 cm. The spectrometer and 
Faraday collector were evacuated to a pressure 
less than 10-* mm Hg by means of a diffusion 
pump. 

The magnetic field was measured by means of 
a flip coil and fluxmeter, the calibration of which 
could be checked with a mutual inductance. Be- 
cause of uncertainty in the average radius, the 
absolute values of Hp given herein may be in 
error as much as 10 or 15 percent, but the results 
of this experiment depend only upon the accu- 
racy of reproducing the magnetic field which is 
estimated to be about } percent. As most of the 
measurements were made near the maximum of 
the spectrum, small errors in reproducing the field 
produced little effect upon the results. 

The Faraday collector F (see Fig. 1) consisted 
of a rectangular box which, except for the bottom, 
was made of aluminum. In the bottom (,’;”’ brass) 
was cut a slit 9 mm X48 mm. The inside surfaces 
were covered by graphite or soot. This box was 
suspended from the roof of a larger box by amber 
insulators so that its bottom was less than 2 mm 
above the exit slit of the spectrometer. The outer 
box served as an electrical shield, and it was made 
air-tight with wax so that it could be evacuated. 
Rough estimates of the efficiency of collecting 
8-particles entering the collector based upon con- 
sideration of the solid angle subtended by the slit 
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Fic. 1. Double 8-ray spectrometer with Faraday collector 
and ionization chamber. 





at the opposite wall and of the reflection coeff- 
cient of the surfaces® indicate that it was at least 
99.5 percent. 

The ionization chamber (see Fig. 1) consisted 
of a copper cylinder 2’ in diameter and 18 cm 
long closed at its lower end by a copper gauze 
with a 3” copper rod mounted coaxially. These 
were supported from the interior of an air-tight 
brass cylindrical case by amber insulators. For 
admitting the beam there was in the bottom of 
this case a hole 13’ in diameter covered by an 
0.1-mm aluminum foil which was held between 
two flat brass rings with Apiezon Q as a gasket. 
The foil could withstand a pressure difference of 


6B. F. J. Schonland, Proc. Roy. Soc. A108, 187 (1925). 
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at least four atmospheres, and when the chamber 
was mounted in place with the spectrometer 
previously evacuated, it could be evacuated for 
filling. The pressure of the gas (argon) was meas- 
ured by the aid of an open ended mercury 
manometer 2 meters long. The concentration of 
the gas (P/T) was maintained between 7.00 and 
8.00 mm Hg/degree absolute (about 3 atmos- 
pheres). All measurements were corrected to a 
standard concentration of 7.50 mm Hg/degree 
absolute by use of some empirical correction 
factors found by determining the dependence of 
the ionization current due to a constant source 
upon gas concentration. The voltage upon the 
chamber was supplied by a 750-volt bank of “‘B”’ 
batteries. With the strongest source ever used, a 
reduction in the voltage to 400 volts caused the 
current to decrease by a factor of 0.997 +0.003 as 
a probable error, indicating that the systematic 
error in the final results due to imperfect satura- 
tion was very small, not more than 0.2 percent. 
The background current was about 410-5 
ampere. The largest source of experimental error 
in the “‘weak source’? measurements was due to 
fluctuations in this background current, mainly 
due to alpha-particle contamination of the cham- 
ber. Numerous precautions were taken to reduce 
this contamination to a minimum. The probable 
error in a determination of the background 
lasting for 450 seconds was 3.4 X 10~* ampere. 
The Faraday collector and ionization currents 
were measured by means of a vacuum tube 
electrometer employing a Western Electric 
D-96475 tube in a modified Barth circuit.”? In 
order to reduce the grid (background) current, 
the entire grid circuit was placed in an evacuated 
shield. This current varied between +6 and 
+12X10-'* ampere, and the largest source of 
experimental error in the Faraday collector meas- 
urements was due to fluctuations in it. The 
probable error in a determination of the grid 
current which required 450 seconds was 5.8 X 10-17 
ampere and seemed to be independent of the 
magnitude of the current. To obtain the most 
stable operating conditions, the electrometer 
circuit ran 24 hours a day, and the storage battery 
which supplied the power was charged by a 
trickle charger with an Acme-Delta voltage 
stabilizer on its input. 
7D. B. Penick, Rev. Sci. Inst. 6, 115 (1935). 
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This electrometer was operated according to 
the null method of Townsend,® as mentioned in 
the introduction. An especial advantage of this 
method is that the capacity of the grid and its 
circuit are not directly involved. 

The three cylindrical condensers used with this 
electrometer were equipped with special shielding! 
so that the capacity between the lead wires was 
eliminated. The smallest of these condensers 
(referred to as C,) was a primary standard of 
11.77+0.03 micromicrofarads. The largest con- 
denser (referred to as C3) was designed to have 
high capacity rather than to have high accuracy 
as a primary standard. Only the ratio of its 
capacity to that of C, is involved in the final 
result, as is explained in the introduction. This 
ratio was obtained by using both condensers to 
measure the same steady current and was equal 
to the reciprocal of the ratio of the corresponding 
rates of change of potential. In this way a value 
of 364.3 micromicrofarads was found for the 
capacity of C3. C,; was also compared by this 
same method with another primary standard C; 
and was found to be in excellent agreement with 
it. The value of C; found from its geometry was 
in good agreement with the calibration by the 
National Bureau of Standards. 

The counters were made of 1” diameter brass 
tubing. On one side was milled a flat place, in the 
center of which was cut a slit 48 mm X6 mm over 
which was waxed a 4.65 mg/cm? aluminum foil. 
Into the ends were waxed closed glass tubes sup- 
porting the counter wire of 3-mil tungsten which 
was kept taut by a spring. To obtain a more 
uniform field the ends of the wire were shielded 
by }” glass tubes which extended into the 
cylinder for 1’. These counters were normally 
filled by a mixture of 9 cm Hg of argon and 1 cm 
of alcohol and were self-extinguishing.* Only 
Apiezon wax was used as this absorbed alcohol 
very slowly, if at all. In the principal experiments 
no counter was used for more than a total of six 
million counts or later than ten days after its 
previous filling. Usually before filling the wires 
were heated in vacuum by passing a current of 
one ampere for several periods of a second or two. 
Only those counters with an essentially flat region 
in their plateaus were used ; otherwise the meas- 
ured counting rate could have no unambiguous 
8 A. Trost, Zeits. f. Physik 105, 399 (1937). 
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connection with the true number of particles 
entering the counter. The absolute efficiency of 
these counters will be discussed later. The circuit 
consisted of a two-stage amplifier followed in turn 
by a pulse sign discriminator,’ a vacuum scale of 
eight,® a thyratron scale of eight,!° and a Cenco 
recorder. The method which was used to correct 
for losses in the counting rate due to the re- 
solving time of the counter has been described 
elsewhere and will not be discussed here."' The 
radium E sources (practically free of radium D 
but containing radium F) were deposited upon 
40-mil nickel wires from a solution containing 
about 7 millicuries of radium D and its products, 
which had been prepared previously in this 
laboratory from old radon tubes. The deposit was 
restricted to a portion of the wire 1 cm long by 
covering other portions with collodion. For work 
with the Faraday collector the source strength 
was usually 2 to 3 millicuries. The procedure used 
in depositing the sources was the usual one” but 
it should be mentioned that the process ap- 
parently depends critically upon the constitution 
of the wire. Pure nickel from Baker and Company 
of Newark, New Jersey gave very good results, 
but some wire from another manufacturer gave 
hardly any yield at all. 

No difficulty was caused by the presence of 
radium F (polonium) in the source because the 
alpha-rays were prevented from entering the 
detecting devices by the magnetic spectrometer, 
and the gamma-rays were too weak to be ob- 
served when the source and counter were 
mounted in their normal positions. 


FOIL CORRECTIONS; INFLUENCE OF SCATTERING 
AND OF SECONDARY ELECTRONS 


The slit on the counter and the exit slit of the 
spectrometer were each covered by a4.65-mg/cm? 
aluminum foil. A certain fraction of the beam is 
lost when the counter is used which is not lost 
with the Faraday collector (unless an equivalent 
amount of foils is placed in front of the col- 
lector). Therefore, if 2’ of formulae (1) and (2) is 

*H. Lifschutz and J. L. Lawson, Rev. Sci. Inst. 9, 83 
(1938) ; H. Lifschutz, Rev. Sci. Inst. 10, 21 (1939). 

10 J, Giarratana, Rev. Sci. Inst. 8, 390 (1937). 

u Y, Beers, Rev. Sci. Inst. 13, 72 (1942). 

12Q. Erbacher and K. Phillip, Zeits. f. Physik 51, 309 


(1928); W. F. Libby and D. D. Lee, Phys. Rev. 55, 252 
(1939). 
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Fic. 2. Foil correction data. The counting rate is 
plotted against increasing foil thickness as extra foils are 
placed in front of the counter. The Geiger counter data are 
extrapolated to zero thickness (solid lines). Also is indi- 
cated (dashed lines) the behavior at small thicknesses as 
inferred from the Faraday collector measurements. 








to be construed to be the observed counting rate 
(corrected, however, for resolving time losses), its 
value must be multiplied by a correction factor J, 

One value of J (hereafter called J,) can be 
found by placing additional foils between the 
counter and the spectrometer slit and measuring 
the corresponding counting rates due to a con- 
stant source. In Fig. 2 are plotted such data for 
the values of Hp which were of interest in the § 
present work. Each point represents the average 
of two values involving 60,000 or more counts 
each. These data are well represented by straight 
lines which have been extrapolated to zero foil 
thickness (solid lines). The correction factor J, is 
found by dividing the extrapolated value of the 
counting rate for zero foil thickness by the value 
for 2 foils. 

Another value (hereafter denoted Jr) of this 
correction factor can be found by measuring the 
current from a constant source, first with no foils 
and secondly with two foils (connected electrically 
to ground) in front of the Faraday collector. Jp 
is given by dividing the first of these quantities 
by the second. Values of Jr (somewhat less 
accurate than the corresponding J, values) were 
obtained at four of the same values of Hp. These 
turned out to be consistently larger than the 
corresponding J, values, the discrepancies being 
too large to be accounted for by errors in the 












































measurements. It can be inferred that the be- 
havior at small foil thicknesses must be according 
to the broken line curves of Fig. 2 and not ac- 
cording to the extrapolations of the straight lines. 

There are two possible causes for a discrepancy 
between the values found by these two methods. 
If there are any secondaries ejected from the foil, 
these will be recorded separately by the collector 
but not by the counter, and J; will be smaller 
than J,. On the other hand, if there are some 
very low energy particles emerging from the 
spectrometer after having been scattered therein, 
these are completely stopped in the first one or 
two foils, and will cause J; to be larger than J,. 
The latter process gives the observed sign of the 
discrepancy, and it must be concluded that this 
is the predominant one. Therefore, it can be con- 
cluded that the value found with the collector 
(Jr) is the more reliable. Instead of attempting 
to measure this quantity accurately, however, 
there was employed in the final experiment the 
much simpler procedure of placing in front of the 
collector during the measurements of Jr a foil 
thickness equal to that used with the counter. 
Then the measured values of Jr and n’ were used 
for the calculation of the charge, and no numerical 
correction was necessary. 

Evidence for scattering was found by photo- 
graphing the beam at a known distance from the 
exit slit. A slight indication that the beam was 
not confined entirely to within thirty degrees to 
the normal was detected. Particles emerging at 
large angles cause the measured value of the 
charge to be too large, for while they fail to enter 
the counter and are not recorded by it, they may 
be recorded by the collector even though they 
strike the outside surface. 

Secondary electrons, if ejected from the foil 
into the collector, make the collector current too 
large, while secondaries ejected from the collector 
produce a reverse effect. Because of the large 
stray magnetic field, it is not to be expected that 
either effect is large, and the number of second- 
aries leaving the collector must be relatively very 
small because of solid angle considerations. To 
investigate the possibility of errors due to these 
sources an experiment was performed in which 
the foils were insulated from the spectrometer 
and in which various potentials were applied to 
them. The application of a positive potential 
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gave no detectable effect, while negative po- 
tentials of 45 and 200 volts both gave an increase 
in current of about one percent, slightly greater 
than the probable error, which indicated that 
there might be secondaries of this amount being 
emitted from the foil in excess of those emitted 
from the collector. 

As the pressure was normally 10 mm Hg or 
less it was unlikely that a large error could have 
been caused by the production of ion pairs. A 
rough calculation based on the specific ionization 
of 8-particles, the length of path, and the pres- 
sure indicated that the probability that a 
B-particle produced one ion pair was about 0.01 
percent. Also an experiment was performed in 
which the pressure was raised to 5X 10-* mm Hg 
with no appreciable change in collector current. 


THE ABSOLUTE EFFICIENCY OF COUNTERS 


The absolute efficiency of a counter has been 
shown to obey the following theoretical formula 
based upon an analogy with the mean free path 
theory of statistical mechanics. (The mean free 
path for producing a primary ion pair is the 
reciprocal of the primary specific ionization.) 


E=100(1—e-**), (3) 
K= (wsdTo/PoT), (4) 


s=number of primary ion pairs produced per 
unit path at standard pressure Py and tempera- 
ture To; P and T are the pressure and tempera- 
ture, respectively, of the gas in the counter; 
d=path of the particle in the counter ; # = proba- 
bility that a primary ion pair causes the counter 
to discharge. 

Some workers"*-'5 have assumed that w was 
equal to unity and have used the measured 
efficiencies of Geiger counters for determining the 
primary specific ionizations of the gas, while 
other workers'* find that w must be considerably 
less than unity if their data are to be consistent 
with values of s obtained by independent 
methods. All are agreed, however, upon a de- 
pendence of the efficiency on pressure of the type 
given by Eq. (3). 

1% W. E. Danforth and W. E. Ramsey, Phys. Rev. 49, 


854 (1936). 
1M. Cosyns, Bull Tech. Ass. CA to, § 173-265 (1936). 


where 


18 T, Graf, J. de phys. et rad. [7], 10, 513 (1939). 
-16J. DeVries and G. J. Sizoo, Physica 6, 593 (1939) 
and references quoted therein. 















If one assumes w is unity and uses the following 
values: s=29.4 per cm at Po=76 cm Hg (see 
reference 14), d=1.9 cm, P=9 cm Hg, and 
T=Tp> (assumed), the efficiency of the counters 
used in this work is calculated to be 99.94 
percent. 

For checking the absolute efficiency experi- 
mentally a method of comparing counters having 
different pressures was used. If two counters are 
used to measure the same beam of particles, the 
measured counting rates are proportional to their 
efficiencies. If E and E, refer, respectively, to the 
efficiencies of two counters which have pressures 
of P and aP but which are identical otherwise, 


E,/E=(1—e-*%?)/(1—e-?). (5) 


In practice e~*? is a very small quantity and if a 
is roughly 2 or greater, this ratio becomes 
1+e-*? approximately. 

Several comparisons were made of the effi- 
ciencies of counters having 1.1 cm Hg alcohol 
vapor in addition to 4.7, 9, and 12 cm argon, 
respectively. While these did not agree within the 
probable error there was no systematic depend- 
ence of efficiency upon pressure. The largest 
discrepancy was observed in one experiment in 
which a counter containing 9 cm argon gave a 
counting rate 2 percent greater than the others. 
Even if such a large change in efficiency could be 
shown to be attributed to changing the argon 
pressure from 4.7 cm Hg (P) to 9 cm Hg (aP), 
e-XP must be not larger than 2 percent and the 
efficiencies of the higher pressure counters must 
be at least 99.9 percent except for the ‘inherent 
error’ (that is, discrepancies between the coun- 
ters which cannot be attributed to a change in 
pressure). 

The “inherent error’ was further investigated 
by comparing a number of counters all having 
the same pressure (9 cm Hg of argon and 1.1 
cm Hg of alcohol), and these were found to 
disagree by amounts large enough to account 
entirely for the discrepancies observed with the 
counters having different pressures. Comparison 
of the five counters used in the final experiment 
of the charge of the 8-particle showed that a 
probable “‘inherent’’ error of 1.2 percent should 
be assigned to each counter. This result was in 
agreement with earlier experiments of this type. 
Therefore, such errors are included in the sum- 
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mary of errors. (See Table III.) Otherwise, the 
counters are considered to have perfect efficiency, 

The explanation of this “inherent” error is by 
no means clear, and this matter could bear 
further investigation. Undoubtedly some of it is 
connected with the fact that no counter has a 
perfectly ‘‘flat” plateau and that there is some 
uncertainty in selecting the correct operating 
voltage. A small amount of it may be due to 
temperature effects in the counter and in the 
voltmeter used to measure the counter voltage, 


PROCEDURE 


Altogether three separate sets of experiments 
upon the charge of the 8-particle were carried 
out. The first two of these, especially the first set 
which was very rough, are to be regarded as 
preliminary experiments. The procedure of the 
final experiment, now to be described, was based 
upon the experience gained from the other two, 

Separate sets of measurements were taken at 
each of five values of Hp, and therefore five 
nearly independent values were obtained. A 
source of radium E having an original strength of 
roughly 3 millicuries was used. The Faraday 
collector measurements (with foils in front of the 
collector) were carried out during the first six 
days. The individual values at each Hp were 
corrected to their average time by the use of an 
assumed half-life of 5.00 days, and the corrected 
values were then averaged together. The ioniza- 
tion current was measured at intervals of about 
every three days until the currents were too 
small to be measured conveniently with the large 
condenser (C3), 25 days after the start of the 
experiment. These values were corrected for 
decay to the average time of the corresponding 
set of collector measurements and then were 
averaged together. Various values of the half- 
life, differing only very slightly from 5.00 days, 
were tried until one was found which gave a 
minimum probable error in the average value, 
and the corresponding average value of Js was 
divided into the average value of Jr to give the 
quantity Ir/Is [see Eq. (1) ]. 

While this quantity was obtained by de- 
termining average values of Ir and Js and by 
calculating their quotient, the analogous ‘‘weak 
source” quantity, n’/Iw, was measured directly 
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by making practically simultaneous observations 
upon ’ and Jw, and the individual values of 
n'/Iw were averaged for calculating the final 
result. In obtaining the individual values, first 
the beam was sent into the counter and a 
reading was made upon it (by using at least 
75,000 counts) while simultaneously the back- 
ground current of the ionization chamber was 
measured. Then the field was reversed and the 
beam was measured twice by using the ionization 
chamber while at the same time the background 
of the counter was determined. Then the field 
was reversed and measurements upon the beam 
with the counter and upon the background of the 
chamber, similar to the first ones, were taken. 
The average value of the counting rate was 
divided by the average value of the beam current 
(after these had been corrected for background) 
to give n’/Iw. In this way corrections for decay, 
which are unreliable with such an old source 
because of the possible presence of impurities, 
were usually avoided, though they were made 
when necessary (never greater than 0.2 percent). 
These measurements were started when the 
counting rate had decreased to 300 counts per 
second, some forty days after the experiment was 
started, and were continued for several days. At 
each of the five values of Hp a separate Geiger 
counter was used. The resolving time of each 
counter was measured three times: at the be- 
ginning, in the middle, and at the end of ‘“‘weak 
source” measurements. The results of the 
resolving time measurements are published 
elsewhere."! 

The largest source of experimental error in the 
collector current and ‘‘weak source”’ ionization 
current measurements was due to fluctuations in 
the background current. If the error caused by 
these fluctuations is assumed to be inversely 


TABLE I. Results of second experiment. 








Charge of 8-particle 


Counter foil Collector foil 
corrections used corrections used 


— 1.627 —- X10-* coulomb 
—1.63, —1.58 
— 1.677 — 1.58; 
—1.685 —1.60, 
—1.69, —1.57, 


— 1.66, — 1.58. 10-?® coulomb 
— 4,99 —4.75 X10" e.s.u. 








Unweighted average 
Or 








proportional: to the square root of the time 
duration of the reading,'’ it can be easily shown! 
that one obtains the greatest accuracy by dividing 
the available measuring time equally between the 
beam and background measurements. It can also 
be shown! that if all the individual values of the 
same quantity are to have the same relative 
accuracy, the time duration of each reading 
should be made inversely proportional to the 
square of the strength of the beam. (Therefore, 
there is little advantage to be gained in carrying 
on measurements on such quantities for more 
than one half-life.) These principles were followed 
throughout these measurements. 

The first preliminary experiment was so rough 
that it is sufficient to say that it gave a result of 
the correct order of magnitude. The procedure of 
the second preliminary experiment is different 
from that of the final experiment in several 
respects. The foils were not placed in front of the 
collector, and therefore a numerical correction 
had to be made. Also only two counters were 
used, but each was used at all values of Hp. As a 
weaker source was used and fewer observations 
were made, the results were not as accurate. 


RESULTS 


The results of the second preliminary experi- 
ment are given in Table I. While the foil cor- 
rections determined with the Faraday collector 
are the more reliable, values of the charge based 
upon both sets of corrections are given for 
comparison. It is to be noted that the values 
obtained with the counter foil corrections increase 
with Hp, while those employing the collector foil 
corrections show no such dependence. The prob- 
able error in the average value is estimated to be 
less than 1.5 percent. 

The results of the final experiment are shown 
in Table II and a summary of the experimental 
error of this experiment is given in Table III. 

The first three items of Table III were found 

17 This assumption is based upon an anal with the 
theory of counting (with a Geiger counter). If B counts 
are observed in a time 7, the counting rate N=@=B/T. 
The mean error in B counts is B4 and the mean error in 
N is B}/T =(N/T)*. Therefore, as N is constant, the error 
is inversely proportional to the square root of the time. 
This assumption is probably more valid in the case of the 
ionization current measurement, for the fluctuations in 
the background are mainly due to the alpha-particle 


contamination of the chamber than in the case of the 
collector current. 
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TABLE II. Results of final experiment. 














n’/Iw Charge =/F/Is 

d Ip Ig Ir/Is X10° (reciprocal coulombs n’/Iw 

Hp (107% amp.) (10-2 amp.) (pure number) X 1018) 10~!* coulomb 
2200 —9.35s —4.11, 2.275 — 1.42, — 1.59, 
2400 — 8.742 —4.64, 1.88. — 1.18; —1.592 
2650 — 8.785 —5.37,4 1.63; — 1.01; —1.61, 
2900 — 8.316 — 5.59 1.48, —0.92, — 1.60, 
3400 —7.012 — 5.30, 1,32, —0.79; — 1.662 

Average — 1.61; 10-" coulomb 


— 4.83, 10-" e.s.u. 


Or 
If the value at Hp 3400 gauss cm is excluded (see text) the average value is — 1.60; x 10- coulomb or — 4.809 10~" e.s.y 








TABLE III. Summary of experimental errors of final experiment. 








Probable error at Hp 








Origin 2200 2400 2650 2900 3400 
Collector current 0.15% 0.07% 0.05% 0.18% 0.23% 
“Strong source’”’ ionization current 0.18% 0.11% 0.11% 0.11% 0.15% 
“‘Weak source’ measurements (n’/Iw) 0.34% 0.55% 0.23% 0.39% 0.30%, 
Error due to error in resolving time 0.15% 0.10% 0.07% 0.17% 0.16% 
inherent error of counter 1.2% 1.2% 1.2% 1.2% 1.2% 
Total (square root of sum of squares) 1.3% 1.3% 1.2% 1.3% 1.3% 





= 





by application of the theory of error to the data. 
The fourth item was determined by the method 
described in another paper." It is to be noted 
that the fifth item, the ‘‘inherent error’’ of the 
counter, is so much larger than the others that 
these can be nearly neglected and therefore the 
five values are essentially equally accurate. The 
probable error in the average value due to these 
items is found by dividing 1.3 percent by the 
square root of 5, giving 0.58 percent. These are 
not, however, “the only experimental errors in- 
volved. Others to be included [see Eq. (2) ] are: 
error in ratio of the condenser capacities 0.1 
percent; error in calibration of small condenser 
0.3 percent ; and the error in the calibration of the 
potentiometer 0.1 percent. The total experi- 
mental error in the average value is 0.7 percent. 

It is debatable whether the Hp 3400 value 
should be included because of the large dis- 
crepancy between it and the other values. If it is 
omitted, the average deviation between the 
remaining individual values is much smaller than 
one has a right to expect on the basis of Table ITI, 
while its inclusion leads to an average deviation 
of just the right size. It should be remarked, 
however, that the application of the theory of 
error to so few values is hardly justified. The 
discrepancy is probably largely due to the par- 





ticular counter used for this value, which gave a 
lower value than the other counters when they 
were all used to measure the same beam. (This 
counter was No. 8 of Table I of the paper giving 
the resolving time data.") If this value is omitted, 
the average becomes — 1.601 X10~'* coulomb or 
— 4.80) 10-" e.s.u. with a probable error of 0.7 
percent, which because of the reason just stated 
is to be considered as the final average of the 
present measurements. 

In addition to these experimental errors, there 
were some systematic errors which are listed 
in Table IV. These have all been discussed 
previously. 


DECAY OF RADIUM E 


In the final experiment the decay of the source 
was followed until 50 days after the start. The 
half-lives used to correct the “strong source” 
ionization were found to predict the activity 
accurately except at the very end when a slight 
amount of radium D content became apparent. 
Because of this impurity these half-lives were 0.1 
percent too large. The corrected values are shown 
in Table V. The probable error of each is about 
0.2 percent, and they may be 0.2 percent too 
large because of a systematic error caused by 
imperfect saturation of the ionization chamber. 
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CONCLUSIONS AND DISCUSSION 


(1) The charge of the -particle is found to be 
— 4,80 9+0.03;x10~- e.s.u., independent of the 
portion of the momentum spectrum from which 
the particles are selected. This value is in agree- 
ment within the uncertainty of the measurements 
with the charge of the electron of — 4.8025; X 10-” 
e.s.u.!8 

(2) The half-life of radium E is found to be 
4.99, days with a probable error of 0.1 percent 
and is found to be independent of the portion of 
the momentum spectrum from which the particles 
are selected. This value is in fair agreement with 
that of 4.975 days found by Curtiss.” 

(3) Geiger-Mueller counters can be used for 
absolute measurements with probable ‘‘inherent”’ 
errors of about one percent. This error was the 
largest single error in this work. Before any 
attempt is made to perform another experiment 
of this type in which counters are used the cause 
of this error should be investigated in the hope of 
building counters which agree with each other 
more accurately. 

The electron multiplier suggests itself as a 
possible alternative to the Geiger counter in such 
an experiment. As the resolving time” is much 
smaller, it might be possible to make the counting 
measurements with the same strength of beam as 
is used for the collector measurements. Another 
advantage is that it requires no foils as it operates 
in a high vacuum. However, much development 
work must be done on this instrument before it 
can be used in such an experiment. 

The ninety-degree type of 6-ray spectrometer 
is much better suited to this work than the 180- 


TABLE IV. Systematic error.* 














Estimated 
Origin upper limit Sign 
Imperfect efficiency of 
araday collector 0.5% Negative 

Scattering — Positive 
Secondary electrons 1% Positive 
Imperfect saturation of 

ionization chamber 0.2% Positive 








* Positive and negative systematic errors are defined, respectively, 
as those which give too large and too small a magnitude of the charge 
of the £-particle. 


%R. T. Birge, Rev. Mod. Phys. 13, 233 (1941). 

19 L. F. Curtiss, Phys. Rev. 30, 539 (1927). 

*Z. Bray, Rev. Sci. Inst. 12, 127 (1941); James S. 
Allen, Rev. Sci. Inst. 12, 484 (1941). 
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degree type, which was used in the present 
experiment, for, according to Stephens,”' the 
ninety-degree type brings the beam to a focus at 
a point outside the magnetic field. Therefore, the 
focal point can be made to be at the center of the 
detecting instruments (counter, collector, and 
ionization chamber) rather than at a point in 
front of their entrance slits. The source should 


TABLE V. Half-life of radium E. 














Half-life 

Hp (days) 
2200 4.975 
2400 4.99, 
2650 4.99, 
2900 5.00 
3400 4.98, 
Average 4.995 








be mounted at the bottom of the spectrometer 
rather than at the top so that active deposits 
could not fall into the spectrometer and con- 
taminate it. 

The method developed in the present work 
could be used for measuring the charge of the 
positive 8-particle (positron), but there would be 
some difficulties not encountered in the present 
work. There would be an inherent error due to 
the annihilation gamma-rays. As it would be 
very difficult to prepare a source entirely free of 
negative particles, one would have to equip the 
apparatus with reliable shutters for measuring 
the background, while in the present experiment 
the background was measured by reversing the 
field, for the shutters which were included in the 
apparatus were not found reliable. 
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The proton energy lost in an inelastic collision with a 
nucleus reappears as an excitation energy of the nucleus. 
The measurement of the energy loss of inelastically scat- 
tered protons is a convenient method of determining 
the energies of nuclear excitation levels. Protons scattered 
at an angle of 135° by a thin scattering foil are detected 
by a proportional counter feeding a scaling circuit which 
is biased to count only the large pulses resulting from 
slowly moving protons. Aluminum stopping foils intro- 


duced between the scatterer and counter are used to slow 
down the scattered protons so that protons of a particular 
energy group move through the counter slowly and are 
counted. A plot of the number of counted protons against 
stopping foil thickness gives a curve in which peaks appear 
for different proton energy groups. Well-defined peaks 
corresponding to inelastically scattered protons are ob. 
tained for a number of substances and excitation energies 
are computed. 





INTRODUCTION 


F a proton approaches a nucleus closely, any 
one of several things may happen. The 
proton may be scattered elastically so that the 
total kinetic energy after collision is the same 
as before. On the other hand, it may be scattered 
inelastically so that the nucleus is left in an 
excited state. Another possibility would be 
nuclear transformation induced by the proton. 
An example of this would be the (f, a) reaction. 
Inelastic scattering can be thought of as a 
(p, p) reaction in which the proton enters the 
nucleus, and later a proton leaves it with a lower 
kinetic energy than that of the original proton, 
thus leaving the nucleus in an excited state. 
It has been shown that in general a (p,m) 
reaction is more likely to occur than a (p, p) 
reaction.! In the case of a (p, p) reaction, the 
proton must penetrate the Coulomb barrier in 
order to leave the nucleus. This is not the case 
for the neutron in a (p, m) reaction. Thus, if a 
(p,m) reaction is energetically possible, one 
would expect to observe very little inelastic 
scattering of protons. The substances most 
. suitable for the observation of inelastic scattering 
then are those which have high (, m) thresholds. 
Heavy elements are, in general, unsuitable 
for inelastic scattering experiments because of 
the difficulty of penetration of the Coulomb 
barrier. In general, the substances most favorable 
for inelastic scattering are light elements with 
high (p, 2) thresholds. 


1V. Weisskopf and D. Ewing, Phys. Rev. 57, 472 (1940). 
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Wilkins and Kuerti,? Wilkins and Wrenshall; 
and Wilkins,‘ using their photographic emulsion 
technique, observed the inelastic scattering of 
protons. From their work levels were found in 
Mg and Al. Their values agree fairly closely 
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Fic. 1. Slit system. 


2T. R. Wilkins and G. Kuerti, Phys. Rev. 57, 108 
(1940). 


*T. R. Wilkins and G. Wrenshall, Phys. Rev. 58, 758 
1940 


(1940). 
*T. R. Wilkins, Phys. Rev. 60, 365 (1941). 
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with those given here. In the case of aluminum, 
there is agreement with a known y-ray energy. 


APPARATUS | 


The University of Rochester cyclotron gives a 
proton beam of approximately 6.9-Mev energy. 
The beam, however, is not completely homo- 
geneous. To get a monochromatic beam of 
protons from the cyclotron, one may bring out 
the beam through a system of slits in the 
fringing field of the cyclotron magnet. With 
three slits it is possible to get as nearly mono- 
chromatic a beam as is desired by making the 
slits narrow enough. The system of slits used is 
shown in Fig. 1. The slit system is open to the 
cyclotron vacuum and the exit slit is covered 
by an aluminum window. 

Figure 2 shows the scattering chamber, 
monitoring ionization chamber, and proportional 
counter used in these experiments. The foil is 
mounted in the center of the chamber which is 
evacuated to eliminate scattering from air around 
the foil. A small ionization chamber is mounted 
at the rear of the scattering chamber for the 
purpose of monitoring the beam. The window 
through which the scattered protons leave the 
chamber is mounted at an angle of 135° from 
the direction of the incident protons. It is 
composed of a thin aluminum foil waxed onto a 
supporting grid. After the scattered protons 
leave this window, they pass through an alumi- 
num stopping foil and enter a proportional 
counter. 

The proportional counter was filled to a 
pressure of 20 cm with commercial argon and 
operated a scale of 32 through a pulse amplifier. 
The scaling circuit could be adjusted to respond 
to pulses of various sizes. 

There are, in general, two methods of range 
measurement. If the scaling circuit bias is 
adjusted so that all protons are counted but 
that the y-ray background is not counted, we 
speak of the ‘plateau method.” If the number 
of protons counted is plotted against absorber 
thickness, we obtain a step-like curve with a 
step at every proton energy group. If on the 
other hand, the bias is adjusted so that only the 
largest pulses are counted, we speak of the “‘peak 
method.”’ Only very slowly moving protons can 
produce pulses which are large enough to be 








counted, and a plot of number of protons against 
number of absorbers gives a curve with peaks 
corresponding to the proton energy groups. This 
is illustrated in Fig. 3, which is a plot of number 
of recorded proton pulses against thickness of 
stopping foils for various scaling circuit biases. 
A remote control system employing selsyn 
motors was used for changing the absorber foils. 
The absorber foils were cemented over holes 
around the edge of a large aluminum disk. This 























Fic. 2. Scattering chamber and counter, 


disk was coaxial with a gear which was driven 
by a one-tooth pinion mounted on the selsyn 
motor shaft. Thus one revolution of the selsyn 
motor changed from one absorber foil to the next. 

The absorber foils were calibrated by com- 
paring their stopping power with that of air, 
by the use of the proportional counter as a 
detector. The foils were built up in multiple 
layers of approximately } mil of aluminum foil. 
It was found that each of these }-mil absorbers 
was equivalent to 1.10 cm of air at 15°C and 
760 mm of mercury. Proton energies were 
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Fic. 3. Effect of scaling circuit bias on scattering curves 
of aluminum. 


calculated from their ranges in aluminum and 
Livingston and Bethe’s curves. 


RESULTS 


Figure 4 is a plot of protons scattered by a 
Pt foil. As would be expected, there are no 
inelastic peaks observable. The background to 
the left of the elastic peak is probably due to 
coincidence between two small proton pulses 
giving the effect of a large one. The ordinate in 
this and the following curves if multiplied by 32 
gives the total number of protons counted for 
any one point. 

Aluminum 

Figure 5 shows the scattering curve for a 
0.35-mil Al foil. As can be seen, there are four 
peaks visible. The energy of the elastic peak was 
calculated from the known energy of the beam, 
correcting for the recoil of the nucleus. The beam 
energy was measured as 6.9 Mev by its range 
in air. Thus, the elastic peak for Al? when 
protons are scattered at an angle of 135°, 
becomes 6.08 Mev. 


5M. S. Livingston and H. A. Bethe, Rev. Mod. Phys. 
9, 268 (1937). 
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TABLE I. Al? scattering. 
Elastic peak Inelastic peaks 
Run energy Mev I II III IV 
1 6.08 5.28 4.22 3.68 — 
2 6.08 5.21 4.13 3.48 — 
3 6.08 5.24 4.21 3.46 —- 
4 6.08 5.31 4.23 3.51 8 
5 6.08 5.27 4.10 3.44 = 
6 6.08 5.26 4.16 3.54 
7 6.08 5.26 4.22 3.65 ~_ 
8 6.08 5.28 4.15 3.60 - 
Average E 6.08 5.26 4.18 3.55 28 
6.08—E=A 0 0.82 1.90 2.53 3.3 
Excitation energy Mev 0.87 2.03 2.70 3.5 
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Using the calculated value for the energy of 
the protons of the elastic peak, the air equiva- 
lence of our stopping foils, and Livingston and 
Bethe’s range energy curves for air, we were 
able to compute the energies of the protons 
corresponding to the inelastic peaks. Table I 
shows a collection of the results of a number of 
scattering experiments. The first column con- 
tains the calculated values of the elastic peak. 
The other columns are the observed energies of 
the inelastic peaks. In two cases we are able to 
observe a fourth inelastic peak. The value for 
the energy of the fourth peak is rather question- 
able. From the observed energies of the elastic 
and inelastic peaks, it is possible to compute 
the nuclear excitation energies corresponding to 
these peaks. In this computation, correction was 
made for the difference in recoil energy of the 
nucleus in the case of elastic and inelastic 
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Fic. 6. Chromium scattering curve. 


TABLE II. Cr® scattering. 














Elastic peak Inelastic peaks 
Run energy Mev I II III 
1 6.46 5.15 3.70 — 
2 6.46 5.21 3.74 2. 
Average 6.46 5.18 3.72 2.9 
we 0 1.28 2.74 3.6 
Excitation energy Mev 1.32 2.84 3.7 








INELASTIC SCATTERING OF PROTONS 
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Fic. 7. Magnesium scattering curves. 


TABLE III. Mg scattering. 














Elastic ak Inelastic peaks 

Run energy Mev I II Ill 
1 5.98 4.77 3.45 2.38 
2 5.98 4.73 3.38 2.40 
3 5.98 4.75 3.5 — 
Average 5.98 4.75 3.44 2.39 
A 0 1.23 2.54 3.59 
Excitation energy Mev 1.32 2.74 3.88 











scattering. The excitation energies for the nucleus 
are shown at the bottom of the table. 


Chromium 


Figure 6 shows a scattering curve for Cr®. 
As in the previous case, this curve shows both 
elastic and inelastic peaks. The foil was made by 
electroplating Cr on copper. The copper base 
was then dissolved away with concentrated 
nitric acid. We are very grateful to Dr. Flagg of 
the Chemistry Department for making these 
foils for us. The results of the Cr® runs are 
shown in Table IT. 


Magnesium 


Figure 7 shows a scattering curve for Mg. 
It is to be noted that the first inelastic peak is 
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Fic. 8. Sulphur scattering curve. 


almost as large as the elastic peak. There are 
two other observable peaks. The scattering foils 
used were 0.4 and 0.6 mil thick. They were ob- 
tained by rolling magnesium ribbon under oil. 
The results are shown in Table III. There is 
some doubt that all the excitation levels can be 
assigned to Mg™ since Mg*® and Mg?® have 
abundances of about 11 percent. 


Sulphur 


Figure 8 shows a scattering curve for S®. 
Sulphur foils were prepared by rolling plastic 
sulphur between sheets of tantalum. 

The foils obtained are quite small and a special 
technique was necessary to hold them in the 
beam. A pyroxylin film was laid down on water 
and the sulphur foil dropped on top before the 
film had hardened. This would break the film 


R. H. DICKE AND J. 
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so that when it hardened it would adhere only 
to the edges of the foil. The film and enclosed 
foil were then removed on a wire frame which 
could be mounted in the scattering chamber, 
The results are given in Table IV. 

A number of other elements were investigated, 
among them were Cu, Ni, Zn, A, Ag, N, and QO, 
No conclusive evidence for inelastic peaks was 
obtained. In Zn there seemed to be faint evidence 
of inelastic peak structure, but this evidence 
was not considered good enough to warrant 
inclusion in these data. The reason for these 
negative results is probably found in the strong 
competition of other processes as (p, ) or (p, y), 

It would have been desirable to try a number 
of other elements. Among the desirable elements 
are Si, Na, Ne, P, and Cl. These should almost 
certainly yield results because of the high 
instability of the product nucleus of a (pf, 2) 
reaction in each case. They were not tried 


TABLE IV. S* scattering. 








Inelastic peaks 








Elastic peak 
Run energy Mev I III 
1 6.2 4.35 2.57 
2 6.2 4.32 2.57 
Average 6.2 4.34 2.57 
A 0 1.86 3.63 
Excitation energy Mev 2.25 4.34 








because of the difficulty of preparing good 
scattering samples. 

The authors wish to express their appreciation 
to the members of the University of Rochester 
Physics Department for helpful advice and 
criticism. 
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The recoil electron spectrum of the gamma-rays from As’* has been examined in a magnetic 


spectrograph. The spectrum was found to be composed of two electron groups corresponding 
to quantum energies of 0.83+0.02 and 1.94+0.04 Mev, the intensity ratio being 3.8 : 1. 
These results suggest excitation levels at 0.8 Mev and 2.8 Mev in the Se’® residual nucleus. 
No contributions to the spectrum of gamma-rays of energy greater than 1.94 Mev were 


observable. 









INTRODUCTION 


HE nuclear radiation of As’* has been the 

subject of a number of recent experi- 
ments,!~-> and numerous conjectures have been 
made as to the energy levels of the Se’® residual 
nucleus resulting from beta-decay by As’*&. A 
gamma-ray spectrograph® has been employed in 
examining the energies and relative intensities 
of the gamma-rays arising from the de-excitation 
of Se*76, 


PROCEDURE AND RESULTS 


As,O; of chemical purity 99.95 percent was 
activated by slow neutrons, the reaction being 
As**—n—vy. As’ is the only stable isotope of 
arsenic. Of the known impurities in the chemi- 
cally pure As,O3, none can form a radioactive 
isotope when irradiated by neutrons with a 
half-period which differs from that of As’® by 
less than a factor of two. The gamma-ray activity 
of the irradiated As,O; was noted from time to 
time by means of an electroscope while the 
momentum distribution of the Compton recoils 
of the gamma-rays was being observed repeatedly 
in the gamma-ray spectrograph. The decay 
curve was followed for about 100 hours, and the 
activity was observed to decay with a half-life 
of 26.9+0.3 hr. in agreement with previously 
determined values'~> running from 26 to 27 hr. 
The decay as observed in the electroscope is 


1P. Harteck, F. Knauer, and W. Schaeffer, Zeits. f. 
Physik 109, 153 (1938). 

?F. Norling, Zeits. f. Physik 111, 158 (1938). 
a 939) Schaeffer and P. Harteck, Zeits. f. Physik 113, 287 

*A.C. G, Mitchell, L. M. Langer, and P. W. McDaniel, 
Phys. Rev. 57, 1107 (1940). 

5G. L. Weil, Phys. Rev. 62, 229 (1942). 

°C, E. Mandeville, Phys. Rev. 62, 309 (1942). 
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plotted in Fig. 1. Before being placed in the 
spectrograph, the As,O; was compressed into 
small cylindrical pellets so as to obtain a con- 
siderable concentration of activity, the effective 
width of the source remaining small. 

A typical curve of the momentum distribution 
of the Compton secondaries of the gamma-rays 
from As" is given in Fig. 2. It is there apparent 
that the spectrum has two strong components. 
The end points of the electron groups correspond 
to gamma-ray energies of 0.83+0.02 Mev and 
1.94+0.04 Mev, respectively. The ratio of the 
intensity of the gamma-ray at 0.83 Mev to 
that of the harder one is about 3.8 :1. This 
momentum distribution was obtained at a 
number of different times over the 100-hr. period 
in which the gamma-ray activity was likewise 
being observed in the electroscope as previously 
mentioned. The shape of the curve was found 
to remain unchanged with time, suggesting the 
presence of only one half-life, in agreement with 
the decay curve of the electroscope. Moreover, 
ordinates corresponding to the same Hp value 
were observed to decrease in magnitude in 
accordance with a half-period of 27+1.0 hr., 
associating the gamma-rays with the decay of 
As” with certainty. A marked advantage of 
this method over the cloud-chamber technique 
when sources of sufficient strength are available 
is that a curve such as that of Fig. 2 may be 
completed in 3 or 4 hr. with high statistical 
accuracy whereas cloud-chamber methods re- 
quire a much longer period. 


DISCUSSION OF RESULTS 


Harteck, Knauer, and Schaeffer! have reported 
gamma-rays at 1.50, 2.15, and 3.1 Mev using 



























to produce the rest masses of a positron-electron 
pair, and it is well known that absorption 


methods give only the value of the energy of the 
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Fic. 1. Gamma-ray activity of AsO; irradiated by slow 
neutrons. The half-life taken from the slope of this curve 
r erin) hr. The activity is therefore ascribed to As”® 

= aX 


hardest gamma-ray present with any accuracy. 
No gamma-rays were observable in the neighbor- 
hood of 3 Mev. Statistical accuracy was very 
poor in that region, however, so that the possi- 
bility of the presence of radiation of low intensity 
at about 3 Mev is not excluded by this experi- 
ment. The quantum energy and intensity of the 
gamma-ray at 0.83 Mev, coupled with the results 
of previous experiments,'~ suggest immediately 
the level diagram for Se” given in Fig. 3. ya and 
yp are there given the values reported in this 
experiment, and yc is suggested as corresponding 
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the method of cloud-chamber pairs. Mitchell, 
Langer, and McDaniel‘ have reported a maxi- 
mum gamma-ray energy of 2.05 Mev as a 
result of absorption experiments. It is, of course, 
not surprising that the intense gamma-ray at 
0.83 Mev was not observed in the two foregoing 
experiments, since its energy is not sufficient 


- slight readjustment in the level diagram. The 


to the 3-Mev radiation observed by Harteck, 
Knauer, and Schaeffer.! 

It is to be noted that the level arrangement 
of Fig. 3 agrees in a qualitative sense with the 
B—y and y—y coincidence experiments of 
Norling? and of Mitchell, Langer, and McDaniel‘ 
in that ya is associated with the beta-ray 
spectrum of high maximum energy and two 
quanta are emitted per disintegration in cascade 
with the beta-ray spectrum of low maximum 
energy. The level diagram also eliminates the 
association of yc with a very soft beta-ray 
spectrum, in agreement with the findings of 
Schaeffer and Harteck.* 

The maximum energies of the two beta-ray 
spectra have been tentatively assumed to be 1.0 
and 3.0 Mev; these are means of the various 
observed and extrapolated values which have 
been previously reported. Any change in the 
assumed values would of course result in only a 


















suggested arrangement of levels and the relative 
intensities of y4 and yz are in approximate agree- 
ment with the value 3.5 : 1 given by Weil5 for the 
ratio of the intensity of the beta-ray distribution 
of maximum energy about 3 Mev to the intensity 
of the less probable distribution of lower maxi- 
mum energy. 
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Fic. 2. Momentum distribution of the Compton recoils 
of the gamma-rays from As’*, The decay of the ordinates 
of this curve in the spectrograph was found to correspond 
to a half-life of 27+1.0 hr. 
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ELEMENT NO. 87 


Professor A. L. Hughes and the cyclotron group 
of Washington University, St. Louis, Missouri. 


Note added in proof: This experiment has been recently 
repeated. Cacodylic acid (Kahlbaum) in the form of an 
aqueous solution has been irradiated by slow neutrons, the 
Washington University cyclotron again being the neutron 
source. H2S was passed into the solution, precipitating the 
activated As ions as a sulfide. AsCl; and HCI had been 
previously added. 

The results were essentially the same, ya and yg again 
being present in about the same relative intensities. A 
better value of ys appears, however, to be 2.00+0.04 Mev. 
Although the source was stronger and more highly concen- 
trated, statistical uncertainty was still rather large in the 
neighborhood of 3 Mev, thus preventing any conclusive 
statement with regard to presence or absence of yc. Had 
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Fic. 3. Level scheme for Se”®. 


ve been present with an intensity less than 20 percent of 
that of ys, it would have very probably escaped detection. 





PHYSICAL REVIEW VOLUME 63, 


NUMBERS 3 AND 4 


FEBRUARY 1 AND 15, 1943 


The Search for Element No. 87 
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An attempt to secure the element 87 La;,2 x-ray lines in fluorescence spectra is made with 
negative results: no trace of the line sought for was found. The fluorescent material was CsHSO, 
derived from lepidolite, in which free thallium has been found. The writer suggests that the 
free thallium is derived from unstable element 87 atoms which had been present in the lepi- 


dolite molecule. 


INTRODUCTION 


HE status of the search for element 87 was 
fully discussed by the present writer in a 
previous paper.' Since then, there has been little 
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Fic. 1. Schematic spectral diagram showing effective 
portion of continuous spectrum emitted by molybdenum 
eer anticathode which might cause fluorescence of 

o portion of fluorescent target to effect registration of 
Mo K lines on plate. 


* Now at the University of Southern California. 
1F. R. Hirsh, Jr., Phys. Rev. 51, 584 (1937). 


further work done on the subject. Hulubei*® has 
published x-ray spectrograms which reveal little 
or nothing. Moreover, the present writer’s ob- 
jections! to Hulubei’s wave-length measurements 
still hold: the wave-lengths are far from the cor- 
rect value; neither is the frequency separation of 
the 87” Lay,» lines correct as given by Hulubei. 
As Hulubei admits,* they are ‘‘mésures trés 
difficiles” and “‘raies trés faibles.”’ 

One might possibly assume, then, that there 
is no evidence for the existence of element 87. 
There is, however, ample evidence for the exist- 
ence of the unstable atom. J. A. Cranston,‘ in 
1913, found that three atoms in 105 of MsTh.(89) 
emitted alpha-particles to form element 87 (see 
Fig. 4). He used the Geiger-Nuttall method. 
G. Guében,® * of Liége, using the Geiger-Nuttall 


*H. Hulubei, Comptes rendus 205, 854 (1937). 

*H. Hulubei, Comptes rendus 209, 675 (1939), 

‘J. A. Cranston, Phil. Mag. 25, 712 (1913). 

5G. Guében, Ann. Soc. Sci. de Bruxelles 52, 66 (1932). 
* G. Guében, Ann. Soc. Sci. de Bruxelles 53, 115 (1933). 
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Fic. 2. Vertical section through spectrograph (note exaggerated scales) drawn 
from measured dimensions of spectrograph, and 





roving that both Mo K lines 


and element 87 La:,2, produced by fluorescence, should register on photographic 


plate, since alignment is correct. 


method, showed the existence of alpha-particles, 
estimated that one atom in 10° to 10° atoms of 
MsTh:z disintegrated to give alpha-particles and 
form an atom of element 87, instead of proceeding 
in the usual manner, with negative electron 
emission to form RaTh(90). He definitely showed 
that the alpha-particle activity had a half-life the 
same as that for the beta- and gamma-activity 
for MsThe. Thus the discovery of the unstable 
atom of element 87 is due to Cranston and 
Guében, with ample supporting evidence. 

Paneth’ has very recently reported the possible 
formation of element 87 by alpha-particle emis- 
sion from Ac(89), forming radioactive AcK(87) 
of atomic weight 223. Meyer, Hess, and Paneth 
(see reference 7) found the alpha-activity, which 
was also found by Mlle. Perey. They estimate 
that one atom in one hundred of Ac(89) dis- 
integrates by alpha-emission. 

The present writer, however, at the instigation 
of Dr. T. G. Kennard of Claremont Colleges, 
California, who had an element 87 concentrate 
to be tested, decided to make a final search for 
the stable atom. 


EXPERIMENTAL 


In order to conserve the concentrate for further 
investigation, it was decided to attempt to excite 
the LZ x-radiations of element 87 by fluores- 
cence. The writer has previously described the 
fluorescence x-ray tube® which was used. The 
Mo primary x-ray anticathode was operated at a 
potential of 25 kv, d.c. Since the Mo K radiation 
is excited at 20 kv, there should be a considerable 


7F. Paneth, Nature 149, 565 (1942). 
8F. R. Hirsh, Jr. and F. K. Richtmyer, Phys. Rev. 44, 
955 (1933). 










beam of Mo K radiation leaving the primary 
x-ray anticathode and incident on the fluorescent 
target. The wave-length of the Mo Kay,» lines is 
0.707A; the Ly absorption edge of element 87, 
secured by interpolation from a Moseley graph, 
is located at 0.812A. Thus the Mo K radiation 
lies on the high absorption side of the Zin edge 
of element 87, and 87, La should be strongly 
excited (La comes from the transition Ly1My). 

The fluorescent target was made of Mo, with 
a three-mm hole in the center, into which the 
concentrate for element 87 was gently melted, 
The Mo-continuous radiation at 25 kv should 
excite the Mo K radiation, from the fluorescence 
target, as reference lines. The reference lines did 
not register on the photographic plate: this is 
not unreasonable, however, as very little of the 
continuous radiation is able to excite Mo K 
radiation (see Fig. 1). The shaded portion of the 
continuous spectrum alone can excite Mo K radi- 
ation by fluorescence. The short wave-length cut- 
off of the continuous spectrum is at 25 kv. Re- 
membering the rule of 12.3 kv/A, we find the 
cut-off lies at 12.3/25+0.5A. Mo K radiation is 
excited at 20 kv : 12.3/20~0.6A. The Mo K radi- 
ation is at 0.7A approximately. Thus it will be 
seen that the Mo K reference lines should not be 
intense if they were registered. Lack of condi- 
tioning time for the x-ray tube made impossible 
the use of a higher potential for the primary 
anticathode. The axial alignment of the fluores- 
cent target was excellent: this was done by 
rotating the glass insulating sleeve while the wax 
seal was hot. Figure 2, carefully constructed from 
the spectrograph and tube dimensions, shows 
that the photographic plate sees both top and 
bottom parts of the Mo part of the fluorescent 
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target, hence the Mo K reference lines should 
register if the fluorescence were intense enough. 
X-radiation was scattered by the fluorescence 
target, however, as evidenced by the presence of 
a central beam on the photographic plate (see 
Fig. 3). In this case: 


hy=hy' + Ev’ +me*/2 ((1—6*)-*-1]+ MV?/2, 


where hy is the original photon, hy’ the scattered 
photon, Ev’ the recoil electron ejection energy, 
me?/2 [(1—6*)-*—1] the recoil electron kinetic 
energy, and MV?/2 the energy of motion im- 
parted to the atom. 

Previously, under approximately the same 
operating conditions, an intense fluorescence 
plate had been secured in 10-20 hours. The plate 
used for testing the presence or absence of ele- 
ment 87 in the concentrate was exposed for 36 
hours; the primary and fluorescence targets were 
cleaned hourly to minimize the absorption of 
tungsten deposits from the filament, which was 
run at an emission current of 15 ma, the primary 
anticathode was operated at 25-kv positive po- 
tential. No trace of element 87 La; lines was 
present on the photographic plate. 

We conclude that element 87 was not present 
in the CsHSQO, concentrate from lepidolite. Al- 
though this simple experiment shows that ele- 
ment 87, in the stable form, was not present in 
the lepidolite from which the CsHSO, was made, 
it does not prove it had not been present. In fact, 
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S7La MoKe C 


Fic. 3. S 
(scattered) 


trogram showing registration of central 
m of x-radiation, and absence of element 


87 lines. A quartz crystal was used: the Bragg angles for 
Mo, Ka (A=0.707A) and 87, La (A= 1.045A) are 5° and 7°, 
respectively (approximately). (The plate faces in the 
direction from which the scattered x-rays came.) 
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Kennard and Rambo® have reported the presence 
of free Tl in this lepidolite molecule. The writer 
wishes to suggest that it is evidence that unstable 
element 87 had been present there. MsTh;2(89) is 
derived by negative electron emission from 
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Fic. 4. Disintegration chart showing possible origin of 
Tl found in the lepidolite molecule. Atomic mass plotted 
against atomic number for the MsThi:—MsThz chain of 
disintegrations. 























MsTh,(88), which in turn is derived by an alpha- 
emission from Th(90), (see Fig. 4). MsTh2(89) 
goes down by an alpha-emission to element 87; 
another alpha-emission produces element 85; 
two more alpha-emissions would produce T1(81) 
of mass 212, and would account very nicely for 
the presence of Tl in the lepidolite molecule. 
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To extend knowledge concerning the positive and nega- 
tive point-to-plane corona, studies were made of the dis- 
charge in very pure He, Ne, and A. Neither pre-onset 
streamers nor burst pulses were observed in pure Hz and 
Ne, in accordance with present theories of these effects. 
At high fields weak streamers developed which ultimately 
led to breakdown. Adding traces of O2 to either Hz or Nz 
at once produced weak pre-onset streamers, and at higher 
concentrations produced burst pulses. Trichel pulses, 
which occur in properly triggered corona discharges in 
gases which form negative ions, are not observed in pure 
H: and N: but are immediately produced upon the addition 
of traces of O2. For these gases the potentials for negative 
corona were observed to be lower than for the positive 
corona. In the case of argon, it was found that with 
positive point potential the first streamer formed was in- 



















tense enough to cause spark breakdown, while with the 
point at a negative potential a heavy arc replaced the 
corona. Contamination of the argon with He or Ng intro- 
duced streamers for the positive corona which showed all 
the characteristics of those formed in air. This effect is 
attributed to the suppression of metastable atoms in the 
argon due to the presence of the impurity. For the negative 
corona small amounts of these gases (0.1 to 0.5 percent) 
prevented arcing and yielded a strong corona with glow 
discharge characteristics but no Trichel pulses. Addition 
of sufficient O2 causes argon to act in all respects like air 
except that the onset potentials were much lower than 
in air. A clean-up of impurities in the negative corona in 
H:2 is shown to be due to formation and absorption of water 
vapor. The etching of negative points in the discharges in 
pure gases appears to be related to the sputtering process. 








INTRODUCTION 





HE application of newer techniques to the 
study of point-to-plane coronas in air has, 
as a result of recent studies,' revealed a series of 
mechanisms active in the individual cases of the 
discharge which depend on photo-ionization and 
space charge formation for their interpretation. 
In order to elucidate the mechanisms and to 
confirm the theories put forth, it seemed desirable 
to extend the studies to pure free electron gases of 
widely different photoelectric characteristics. 
This was especially important in view of the fact 
that one of the mechanisms, streamer formation, 
appears to be responsible for the sparking mecha- 
nism in uniform gaps above pd=200 (mm-cm) 
in air. 

Investigations of the positive and negative 
point corona phenomena in the highly purified 
gases He, Ne, and A were undertaken. The re- 
sults obtained indicated the value of studies with 
the admixture of certain impurities. 

The results of W. H. Bennett? on the purifica- 
tion of Hz by the negative corona and the at- 




























1G. W. Trichel, Phys. Rev. 55, 382 (1939); L. B. Loeb, 
and A. F. Kip, J. App. Phys. 10, 3 (1939); L. B. Loeb, 
Fundamental Processes of Electrical Discharge in Gases 
(John Wiley & Sons, New York, 1939), p. 520. 

2 W. H. Bennett, Phys. Rev. 58, 992 (1940). 
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tendant phenomenon of etching of points made it 
desirable to investigate these questions. 

All measurements were made in Pyrex tubes, 
highly outgassed by being baked out, subjected 
to glow discharge in an induction furnace, and 
the pure gases then introduced into the tubes 
which were sealed off subsequent to filling. A 
cross-sectional diagram of a typical tube is shown 
in Fig. 1. Pt points, cylindrical in form with a 
hemispherical end carefully polished under the 
microscope with tin-oxide were used in the 
measurements. The platinum plates were 4 cm in 
diameter, with tungsten leads. The gap length 
ranged from 3.1 to 4.6 cm. The walls were 
shielded by evaporated thin gold films and all 
metal parts were heated to dull redness for ap- 
proximately three hours with an induction fur- 
nace. Unless otherwise stated, all work was done 
at 760-mm Hg pressure and 22°C. 

The tubes were filled with tank He, No, or A 
that had been passed over a conventional puri- 
fying train and liquid air, mercury being carefully 
excluded. In the case of N2 a W filament in the 
sealed off tube was heated to incandescence for 
20 hours to remove the last traces of Oz. 

The d.c. high potential power supply was of 
the same type as the one employed by Sanders; 

*F. H. Sanders, Phys. Rev. 41, 667 (1939), Fig. 1. 
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it was stabilized by an electronic stabilizer of the 
saturable core reactor type to better than one 
percent above 5000 volts. The discharge was ob- 
served visually by a telemicroscope (70), and 
the electrical pulses were recorded by an oscillo- 
scope! which was connected across a resistor 
between the plate electrode and ground. 


POSITIVE CORONA 


In order to analyze the behavior of the dis- 
charges in pure free electron gases and in con- 
taminated gases it is necessary to recall some of 
its salient features as found in air. 

The streamer® which characterizes the positive 
corona in air depends on the first Townsend 
coefficient, the field strength, and mainly the 
efficiency of photo-ionization in the gas. The 
positive ion space charge left behind by a 
streamer tends to choke off further streamer for- 
mation and gives rise to the laterally spread burst 
pulse corona. The crucial phenomenon, photo- 
ionization in the gas, occurs most efficiently in 
mixed gases or in gases where metastable states 
are abundantly produced. 


Experimental Results in H, 


In He with a gap of 3.1-cm corona onset oc- 
curred at 3500 volts, showing a very feeble 
localized light spot of 0.005 cm in diameter on the 
point. In air under the same geometrical con- 
ditions one observes a corona glow spread uni- 
formly over the tip of the point. In the potential 
region of corona onset, streamers were not ob- 
served, and the discharge seemed to be a suc- 
cession of electron avalanches which formed in 
the relatively localized region of highest field 
strength at the tip of the point. Photo-ionization 
in the gas seemed inadequate to spread the glow 
over the whole tip of the point. At corona onset 
no pulses of definite character could be detected 
with the oscilloscope. With a 4-stage amplifier an 
uninterrupted succession of small irregular kicks 
was observed which may be attributed to the 
avalanches approaching the point. 

At 9000 to 10,000 volts the fields were such 
that weak streamers could form which increased 


‘A. F. Kip, Phys. Rev. 55, 549 (1939), Fig. 1B. 
* A. F. Kip, Phys. Rev. 55, 545 (1939). 
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in length and intensity with increasing fields and 
at higher potentials led to a spark.® 

The addition of varying amounts of Oz» (less 
than 0.1 percent to 1 percent) to the Hz immedi- 
ately resulted in the formation of pre-corona 
onset streamers. These increased gradually in 
strength and number on raising the potential 
such that they could be seen visually as what is 
later in Table I termed a “‘streamer-corona.”’ On 
further raising the potential the streamers finally 
led to a spark. Increasing the percentage of im- 
purity again gave pre-onset and onset streamers 
which had a higher intensity than those for lower 
concentrations of O, in H». At potentials above 
corona onset further streamer formation was not 
observed, and the discharge degenerated into a 
burst pulse corona as described by Trichel and 
Kip. At still higher potentials breakdown stream- 
ers occurred. The potential range between the 
appearance of the burst pulse corona and the 






































Fic. 1. Cross section of a typical discharge chamber, 
showing the point-to-plane gap which can be viewed 
through a plane window, the auxiliary W point in the 
upper left-hand corner, and the evaporated gold films 
which served shielding purposes. 






* These streamers are later referred to as ‘‘breakdown 
streamers.” 














breakdown streamers is broadened with in- 
creasing amounts of O2 impurities in He. Subse- 
quent additions of O, led to an extension of the 
discharge area over the point. The currents in He 
contaminated with O2 were smaller than those in 
pure He (Fig. 2). 

For the interpretations of these results it seems 
profitable to recall that photo-ionization in the 
gas occurs most readily in mixed gases or in gases 
in which metastable states are abundantly pro- 
duced. The intense photo-ionization in a mixed 
gas is due to the fact that one constituent has an 
ionization potential which is lower than the higher 
states of excitation of another constituent. In a 
pure gas such as H2 one should expect very little 
photo-ionization. Since photo-ionization in the 
gas is the main mechanism to propagate a 
streamer from the high field region close to the 
point into the low field region several millimeters 
away it is not surprising to find that no streamers 
were observed at or before corona onset. The 
addition of 0.1 percent O2 or more to H: realizes 
the state of a mixed gas, the efficiency of photo- 
ionization increases with increasing amounts of 
Os, as is borne out by the fact that the intensity 
of streamers in successive mixtures also becomes 
greater. One percent Oz in He produced such 
intense streamers at corona onset that they were 
then able to leave a sufficiently large positive 
space charge behind such as to choke off further 
streamer formation by means of a considerable 
reduction in field intensity. The discharge then 
chose the way of least resistance and spread 
laterally over the point in the form of a burst 
pulse corona. To some extent this lateral ex- 
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Fic. 2. Current-voltage characteristics for positive corona 
in pure and impure Hg, gap length 3.1 cm. 
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tension of the glow is facilitated by the diffusion 
of negative ions which in turn produce triggering 
electrons for the formation of new avalanches the 
sum of which then constitutes a burst pulse. 

When the field is increased to the point where 
the positive ion space charge region is cleared in 
time to receive new initiating electrons from the 
burst pulse glow, and it is increased such as to 
counterbalance the space charges, streamers will 
again appear. Under these conditions of higher 
fields the streamers will be sufficiently vigorous 
to cross the gap, form a cathode spot, and start 
an arc- or spark-breakdown. 

It should be noted in Table I that the onset of 
an avalanche corona in pure He occurs at 3500 
volts, while the appearance of a streamer corona 
with 0.1 percent O2 in He requires 5000 volts, 
Subsequent additions of O: increased the effi- 
ciency of photo-ionization and therefore the ease 
of streamer and burst pulse formation, which 
lowered the threshold potential of corona onset. 

The smaller currents in contaminated Hz may 
be due to changes in Townsend coefficients and 
changes in ionic mobilities. A change in mobility 
would be more pronounced with O2 in H, than 
in N2 or A, owing to the great difference in ionic 
mobilities of O2 and He. The ratio of mobilities in 
H: and O; is about five. 


Experimental Results in N; 


In pure N2 the conditions were not materially 
different from those in Hz: No pre-corona onset 
streamers were observed. Corona onset occurred 
at 4800 volts in form of a localized spot on the 
positive point. The oscillographic pattern was the 
same as in the case of He, namely, a sequence of 
irregular fluctuations of the electron avalanche 
type. Increasing the potential resulted in a 
brighter corona, which extended further toward 
the cathode in form of a continuous luminous 
channel, possibly due to the action of metastable 
atoms or active Ne. The potentials at which 
various phenomena appear is increased above 
those of He. 

At 11,000 volts very faint streamers appeared, 
which give the well-known kicks in the oscillo- 
scope superimposed upon the continuous corona 
ripple. In both N2 and Hz breakdown streamers 
could go across the gap at relatively low po- 
tentials without causing breakdown; but at 
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POINT-TO-PLANE CORONA 


TABLE I. Essential data of the positive point-to-plane corona. 























Corona 
Cur- Onset of Breakdown 
Gas Pre-onset Onset Visual Burst | rentt breakdown potential 

streamers | potential character pulses | wA streamers approx. 
H, puret none 3500 v localized* none | 4 10,000 v 18,000 v 
H.+ 0.1% O2 yes 5000 v intermittent none | — 8,000 v — 

streamer-corona 

H.+ 0.5% O» yes 4500 v streamer-corona none | 3 9,000 v — 
H.+ 1% O2 yes 4000 v uniform** small | 3 15,000 v 18,000 v 
N; puret none 4800 v localized none | 1 11,000 v 20,000 v 
N2+0.2% Oz incipient | 4500 v streamer-corona none | 1 9,300 v ~ 

streamers 
N2+1% O2 incipient | 4800 v localized, for yes 1 12,500 v 20,000 v 

streamers higher v uni- 

form glow 

Room-airt yes 5500 v uniform yes =: 15,000 v _ 
A purett none none none none | — — 3,500 v 
A+0.1% He yes 3300 v localized none | 0.5 as — 
A+0.5% H: yes 3500 v localized small | 0.5 — — 
A+1% H: yes 3450 v localized small | 0.5 — — 
A+0.3% N:2 none none none none | — 3,500 v 3,500 v 
A+0.5% Neo yes 3600 v streamer-corona small | — 4,600 v 4,600 v 
A+1% N2 incipient | 3500 v localized small | 0.5 6,500 v 6,500 v 

streamers 
A+0.1% O2 yes 3800 v streamer-corona none | — 4,100 v 4,100 v 
A+0.4% O2 yes 4500 v streamer corona small | 0.6 —~ 7,000 v 
A+1% O2 yes 4000 v uniform small | 0.6 — 15,000 v 











** Uniform glow, caused by the burst-pulse corona. 
+ Gap length 3.1 cm. 
tt Gap length 4.6 cm. 


higher potentials the number of streamers in- 
creased considerably until in Nz at 20,000 volts 
they became strong enough to initiate spark 
breakdown. 

The effect of O2 impurities in Ne (less than 0.1 
percent up to 1.5 percent) was studied. The 
sequence of streamer and burst pulse phenomena 
closely duplicated those in H». The visible corona 
also changed from a localized glow in pure N: to 
a uniformly spread glow with increased amounts 
of Oz. The magnitude of currents did not change 
materially by addition of O2 to Ne. 

Although metastable states and activated 
forms of atoms exist, their abundance seems to 
be too small to increase the probability of photo- 
ionization with respect to that of He to any 
appreciable extent, since no streamers were 
formed in the region of corona onset. The first 
Townsend coefficient takes on appreciable values 
in Ne at distinctly higher values of the field 
strength as compared to He, which accounts for 
the differences in starting potentials. The com- 
plete absence of pulse formation in the pre-onset 
and onset region of both pure Hz and Nz means 
that no Geiger counter regime exists in those 


* Localized glow, due to electron avalanches approaching the point in the highest field region. 


¢ Current at a potential which is about 20 percent above the corona onset potential. 










gases and confirms the conclusion that Geiger 
counter action is determined by space charge 
action.’ 


Experimental Results in A 


In argon a spark crossed a gap of 4.6 cm at 
approximately 3500 volts. Neither pulses of any 
kind nor visible corona occurred. Since both H: 
and Ne are known to quench the metastable 
states of argon by collisions of the second kind, 
varying amounts of these two gases were tried in 
argon. When 0.3 percent N2 was admitted, a very 
few strong streamers occurred at 3500 volts 
followed by a spark at the same potential. An 
increase of N: to 0.5 percent gave streamers in 
the corona onset region which disappeared above 
4100 volts and a localized corona with occasional 
burst pulses was observed. A spark, preceded by 
a few strong breakdown streamers, broke down 
the gap of 4600 volts. Increasing the percentage 
of Ne in argon to values above 0.5 shifted the 
breakdown potentials to higher values. Argon 


7C. G. Montgomery, and D. D. Montgomery, J. Frank. 
Inst. 231, 447 (1941). 
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Fic. 3. Current-voltage characteristics for positive 
corona in pure and impure argon. The arrows indicate the 
upper limit of the current-voltage curve and the sparking 
threshold for a specific mixture of O: in A. 
contaminated with Hz showed essentially the 
same sequence of phenomena. 

The addition of O2 to argon again resulted in 
the formation of pre-corona onset and onset 
streamers without causing breakdown, while 
larger amounts of O2 shifted the breakdown to 
higher potentials. The only difference was that 
above onset streamers still existed in the presence 
of a burst pulse corona (Fig. 3). 

The current-voltage characteristics did not 
change materially for different amounts of Ne or 
H; added to argon and were of the same type as 
those shown in Fig. 3 for O2 in argon. 

In argon the large production of metastable 
atoms distinctly enhances photoelectric ioniza- 
tion in the gas. In addition the first Townsend 
coefficient reaches effective values at very low 
values of the field strength, because electron 
impacts are highly elastic. For these reasons the 
first streamer was strong enough to cause a 
spark. The quenching of metastable states by He 
and Ne means a reduction of the efficiency of 
photo-ionization ; as a consequence we were able 
to observe streamers which were stronger than 
commonly observed though not strong enough to 
cause an immediate breakdown. This is further 
supported by the fact that larger amounts of 
impurities destroy more of the metastable states, 
thus inhibiting formation of breakdown streamers 
at lower potentials and shifting breakdown to 
successively higher potentials. 

Apparently the destructive effect of O2 on 
metastable states is so great that streamers are so 
weakened at low potentials that they do not 
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succeed in choking themselves off. With higher O, 
concentrations streamers were observed together 
with a burst pulse corona at higher potentials 
which indicates some choking of the streamers 
due to increased space charge formation. The 
presence of negative ions produced by O2 appears 
to help spread the glow over the tip of the point 
at higher fields as distinguished from the case of 
H2 and N¢ in argon. 


Low Pressure Phenomena in H2, N2, and A 


These gases were also studied at pressures of 
300 and 100 mm Hg. The sequence of phenomena 
was the same as that at atmospheric pressure 
except that corresponding onset potentials were 
lowered as the pressure was lowered and that the 
visual appearance of the corona tended to become 
more diffuse (Fig. 8). 

It is interesting to note that despite the greater 
diffusion at lower pressures, streamers can still 
form and propagate. This differs from streamer 
formation in sparking in a plane parallel gap in 
air? and must be ascribed to the high fields 
around the point which makes streamer forma- 
tion still possible when it cannot take place in a 
plane parallel gap. The current-voltage charac- 
teristics were such that the slope of the curves 
increased with decreasing pressure due to the 
higher mobilities at lower pressures (Fig. 4) as 
observed by J. Zeleny® in gases of less purity. 
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Fic. 4. Current-voltage characteristics for itive 
corona in pure N2 at 3 pressures. I: p=100 mm Hg, II: 
p=300 mm Hg, III: p=760 mm Hg. The end points of 
curves J and J/ indicate spark breakdown, curve JJJ ends 
at approximately 23 kv with a spark gap length 4 cm. 


8 L. B. Loeb, Mechanism of the Electric Spark (Stanford 
University Press, 1941). 

%J. Zeleny, Phys. Rev. 25, 304 (1907); J. Zeleny, J. 
Frank. Inst. 232, 23 (1941). 
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POINT-TO-PLANE 


NEGATIVE CORONA 


The intermittent relaxation oscillator-like 
character of the negative corona has been as- 
cribed to the formation of negative ion space 
charges'® due to electron attachment to O2 mole- 
cules.'! In the very complex sequence of ioniza- 
tion phenomena in the concentrated field region 
owing to the differences in mobilities of positive 
and negative ions the negative ion space charge 
cloud acts to choke off the discharge until it is 
dissipated several millimeters into the gap. Thus 
we should expect to find the regular Trichel 
pulses appearing in all gases forming negative 
ions, the pulses becoming regular if the necessary 
triggering electrons” can be provided, either by 
field emission from rough points, fine dust 
specks,'* thermionic emission, or photoelectrons. 
In the absence of electron attachment we should 
not expect to observe discharges interrupted by 
space charges since the high electron mobilities 
would preclude the formation of negative space 
charges. 

For the same conditions as for the positive 
point the following experimental results were 
obtained. 


Experimental Results in H, 


In pure H2 negative corona onset occurred at 
2600 volts, which is notably lower than the onset 
for the positive corona. The discharge began as a 
continuous discharge with no Trichel pulses or 
other discontinuities and appeared to be spread 
at onset uniformly over the tip of the point. This 
form of discharge was exceedingly unstable and 
sooner or later contracted to a concentrated point 
discharge in the region of highest field strength. 
The change was accompanied by a change in 
current by a factor of 10%, it occurred more 
rapidly at lower pressures and more rapidly if the 
potential was raised slightly above the onset 
value. At 2800 volts the discharge appeared in 
the concentrated form within a few seconds, 
showing no indications of the Trichel mechanism 
and had all the characteristics of a glow discharge 
with Crookes’ and Faraday’s dark spaces, nega- 





0G. W. Trichel, Phys. Rev. 54, 1078 (1938). 
"N. E. Bradbury, Phys. Rev. 44, 883 (1933). 
2 L. B. Loeb, A. F. Kip, G. Hudson, and W. H. Bennett, 
Phys. Rev. 60, 714 (1941). 
8H. Paetow .Zeits. f. Physik 111, 770 (1939). 
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tive glow, and positive column. The onset of this 
discharge only is registered in Table II. At lower 
pressures or higher potentials the Crookes’ dark 
space became distinctly visible. In the region 
from onset of the high current corona to break- 
down two interesting facts were observed : 

1. The concentrated discharge above onset 
tended to move continually over the surface of 
the point. 

2. The form of the corona changed frequently 
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Fic. 5. Current-voltage characteristics for negative 
corona in pure and impure Hp. I: pure He, III: H2+0.7 
percent Oz, II: after clean-up run on auxiliary W electrode 
as discussed in the paragraph on clean-up corona. Lower 
end of the fully drawn curves marks onset of steady 
corona, dotted lines mark intermittent corona. Gap length 
3.1 cm. 


and unsystematically from the concentrated glow 
to one which was uniformly spread over the tip 
of the point without change in current. 

While the localized glow if it was not changing 
positions did not reveal any pulses, the uniform 
glow at the higher currents showed a succession 
of irregular inductive pulses of no definite charac- 
ter in the oscilloscope. On microscopic examina- 
tion of a point run on a negative corona in form 
of a uniform glow with high currents, it was 
found that the whole tip over which the glow was 
visible was pitted with innumerable holes. 

The introduction of the smallest amounts of O» 
immediately produced Trichel pulses, the onset 
of corona was 4200 volts, the currents were 
lowered to approximately one-half of their value 
in pure H, (Fig. 5). On increasing the amount of 
Oy, the onset of the steady corona was preceded 
by a region of interrupted Trichel pulses which 
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TABLE II, Essential data on the negative point-to-plane corona. 

















Localized corona 
Breakdown 
Gas Onset Onset Visual Trichel potential Peculiarities 
potential current wA character pulses approx. 

Hot pure 2800 v 140 unsteady glow* none — eeveation uniform + localized 
glow 

H.+0.1% O2 | 4200 v 70 unsteady glow ves — same as above 

H.+0.6% O2 | 5800 v 2 unsteady glow yes — same as above 

Not pure 3700 v 90 steady glow** none 14,000 v 

N2+0.1% O2 | 4300 v 50 steady glow yes — 

N2+1% Oz 5500 v 0.1 intermittent*** yes -- 5,000 to 7,000 v intermittent 
corona 

Room-air 8000 v 0.1 intermittent yes = 8,000 to 11,000 v intermittent 
corona 

Att pure — _— — _ 3,500 v no corona 

A+0.3% Nz | 2500 v 200 steady glow none 4,100 v 

A+1% Ne 3900 v 250 steady glow none — 

A+0.1% He | 3200v 400 unsteady glow none 4,100 vy 

A+1% i, 3300 v 300 unsteady glow none — 

A+0.1% O2 | 3000 v 1 steady glow yes — 

A+04% O2 | 3000 v 0.1 intermittent yes 24,000 3,000 to 6,000 v intermittent 
corona 

A+1% O2 3300 v 0.1 steady glow yes — 








* See experimental! data in He. 
** See experimental data in Ne. 
*** Due to lack of triggering electrons. 
t+ Gap length 3.1 cm. 
tt Gap length 4.6 cm. 


began at 5800 volts for 0.6 percent O2 and was 
caused by a lack of triggering electrons.'? The 
steady corona began at 7500 volts. With in- 
creasing potential the frequency of the Trichel 
pulses at onset increased. Larger amounts of O, 
required higher breakdown potentials and caused 
lower currents. 

The difference in onset potentials for the posi- 
tive and negative corona in pure He is very 
distinct from those in air. In air the efficiency of 
photo-ionization is such that the discharge 
mechanism operates through electrons generated 
in the gas by photo-ionization in preference to the 
liberation of electrons from the cathode by 
positive ion bombardment. Since photo-ioniza- 
tion also determines streamer formation from 
positive points, we can explain the nearly identical 
onset potentials for both positive and negative 
corona in air by the fact that one mechanism is 
dominating in both cases. In the case of pure He 
with inefficient photo-ionization in the gas the 
mechanism for maintaining the discharge from 
the positive point is different from that main- 
taining the negative corona as will be discussed 
more in detail in the paragraph on Positive and 
Negative Currents. Consequently the onset 
potentials for both signs of corona differ con- 


siderably. 


Because negative ions are not formed in pure 
Hz no Trichel pulses were observed. 

The two above-mentioned facts on the be- 
havior of the high current corona between onset 
and breakdown can be explained in the following 
way : In the event of insufficient photo-ionization 
secondary electron liberation would occur from 
the negative point by positive ion bombardment. 
The work function of the point and its con- 
ditioning would greatly influence electron emis- 
sion, and the discharge would be sensitive to such 
alterations. Extreme variations of work function 
from 6 volt to 4 volt have been observed" be- 
tween O-coated Pt surfaces and He-coated Pt 
surfaces. The continuous bombardment of Pt by 
positive ions effectively alters the work function 
even in a given gas. This instability of the surface 
would explain the movements of the concentrated 
discharge over the surface, and, together with the 
microscopically observed pitting, it is suggested 
that the discharge form of a uniform glow with 
high currents is caused by the action of one or 
more simultaneous localized spots, constantly 
changing their positions on the tip where every 
change was accompanied by an inductive kick in 
the oscilloscope. These changes as judged from 


“C. W. Oatley, Proc. Phys. Soc. 51, 318 (1939). 
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the approximate frequency of the pulses occurred 
so rapidly, of the order of several thousand for 
the time of persistence of vision, that the obser- 
vations gave the illusion of a uniform glow. The 
validity of this viewpoint also seemed to be indi- 
cated by the fact that the current did not change 
during those transitions, being much too high for 
a Trichel mechanism, and that the uniform glow 
faintly showed a cathode glow, Faraday dark 
space, and positive column which extended 
laterally over the whole tip. 

The interpretation of the very first mentioned 
transition of a low current uniform glow to the 
concentrated discharge with increased current, 
which was a continuous phenomenon, must lie in 
the conditioning of some spot in the high field 
region by positive ion bombardment leading to 
the formation of an adequate cathode spot. 


Experimental Results in N, 


The sequence of phenomena in pure N2 is 
essentially the same as in pure He. The onset 
potential of the localized corona was 3700 volts, 
which is higher than in H, but still lower than the 
onset of the positive corona in Ne. A schedule of 
the potentials at which various phenomena occur 
is given in Table II. As in He the discharge began 
as a uniform glow at 3600 volts (not mentioned in 
Table II), it contracted in time visibly to the 
concentrated discharge with resultant increase 
in current by a factor of 10°. The gradual transi- 





1000 





Sat) a 


| 


3 yi is @§ gg g 


—_> KV 


Fic. 6. Current voltage characteristics for negative 
corona in pure and impure Noe. I: pure No, Il: Ne+0.8 
percent Oz», III: after clean-up run, refers to the paragraph 
on clean-up corona. Gap length 3.1 cm. 
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tion of the discharge could be observed visually 
and the current increase accompanying this visual 
change noted. The localized corona differed from 
the case of He only in remaining steadily at one 
position instead of moving over the surface. 
The effects of small amounts of O:2 ranging 
from 0.1 percent to 1.5 percent introduced into 
pure N» were also studied. The smallest percent- 
age of O2 resulted immediately in Trichel pulses. 
Their choking action caused the currents to drop 
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Fic. 7. Current-voltage characteristics for negative 
corona in pure and impure A. I: pure A, II: A+ <0.1 
percent Os, III: A+-0.1 percent Oo, IV: A+-0.5 percent Oz, 
V: A+1.5 percent O2. The low potential end of the curves 
marks the “offset” potential of the steady discharge. Gap 
length 4.6 cm. 





to approximately one-half of the magnitude in 
the pure gas (Fig. 6). The behavior in general was 
similar to that in He. 

It is obvious that in the conditioning of the 
point surface there must be two mechanisms that 
work : One occurs in both Hz and N: and appears 
to be a conditioning of the rested’® point by 
positive ion bombardment giving increased elec- 
tron emission. The other appears in H2 only and 
results in a de-activation of the point by strong 
H.-ion bombardment. 

The lower onset in H, as compared to N; is to 
be expected from the difference in value of the 
first Townsend coefficient in the two gases. 

8 The initial uniform glow which transformed to the 
concentrated high current form did not appear on a 
surface after a discharge had been freshly interrupted. It 


apparently took a rest period of a few seconds to reform 
the gaseous film causing the uniform glow at low currents. 
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Fic. 8. Visual forms of corona discharges under indicated 
conditions in the onset region. Pressures in mm Hg. The 
forms in the second row show the uniformly spread dis- 
charge at low currents which contracts sooner or later to 
the concentrated forms of discharge with high currents 
which can be seen in the last row. 









Experimental Results in A 





For the same reasons as were mentioned in the 
discussions of the positive corona in argon we 
found that no matter whether the gap distance 
was 3.1 or 4.6 cm an arc discharge set in at 
approximately 3500 volts. It was continuous 
across the gap and no current fluctuations or 
pulses were observed. The current was of the 
order of magnitude of 10* microamperes. Once 
the arc had kindled, it persisted even with 
decreased potentials to a definite “‘offset’’ po- 
tential. The current voltage characteristics from 
onset to offset are shown in Fig. 7. Increasing the 
gap length had the effect of prolonging the time 
necessary for the arc to build up. This lag was of 
the order of seconds for a 3.1-cm gap and tens of 
seconds for a 4.6-cm gap. It must be pointed out 
that these results can only be achieved in argon 
of the highest purity as shown by the effects 
produced by traces of added gases. 

The addition of 0.3 percent of N» to argon did 
not result in an arc immediately, but corona 
onset occurred at 2500 volts: Once started it con- 
tinued to operate on reducing the potential to an 
offset at 1300 volts. On increasing the potential 
a few hundred volts above onset, a continuous 
transition to an arc was observed. Admitting 
larger quantities of N»2 merely shifted the arcing 
potential to higher values above corona onset. 
This may be ascribed to the quenching action of 
the metastable states by No. 

Essentially the same succession of phenomena 
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occurred on addition of Hz to argon. There is one 
notable difference with this mixture which is the 
erratic movement of the corona glow from one 
spot on the tip of the point to another in rapid 
succession. This can also be seen oscillographically 
in the production of irregular inductive pulses. 
These indicate brief interruptions of the current 
but are definitely not due toa Trichel mechanism, 
since the currents were orders of magnitude 
larger than those connected with Trichel pulses. 
This motion of the discharge is again character- 
istic of the Pt point in the presence of He as 
already stated. It does not occur with No in 
argon. Increasing the amount of impurity of 
either He or Nein argon decreases the slope of the 
current-voltage characteristics. 

Admission of varying small amounts of O, 
(less than 0.1 percent to 1 percent) immediately 
initiated regular Trichel pulses, and the break- 
down was at once raised to 10,000 volts, or more 
if the Os-content exceeded 0.3 percent. Adding 1 
percent of Os, altered the appearance of the 
corona at onset, apparently through the nearly 
complete destruction of metastable states. In 
place of the regular Trichel pulses found with 
small amounts of O: present, due to the triggering 
electrons from the few metastables still existing, 
the larger amounts of O2 caused the fluctuating 
corona observed by Hudson” on smoothly 
polished points to appear, when no triggering 
electrons from metastables were present. The 
currents in argon in the presence of Os were 
reduced below those in the presence of He or N, 
in argon by a factor of approximately 10-* as a 
result of the action of the negative space charges 
inhibiting the free flow of electrons and ions 
(Fig. 7). 

As already mentioned above, in gases with 
efficient photo-ionization we find the same onset 
potentials for both signs of corona, because in 
both cases we find that electron generation in the 
gas is the primary mechanism for the mainte- 
nance of the discharge. Pure argon confirms this 
even stronger than air with its spark breakdown 
without preceding corona for the positive point at 
3500 volts and its arc breakdown also without 
preceding corona for the negative point at 3500 
volts. 

The phenomenon of an offset markedly lower 
than the onset potential is characteristic of all 
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negative corona phenomena. In He, No, and air 
the difference in potentials amounts to a few 
hundred volts for a 3-cm gap, while in argon it 
extends over a range of more than a thousand 
volts, probably due to the action of metastable 
states. 


Low Pressure Phenomena for Pure 
Hp, N:, and A 


In the case of the negative corona a decrease 
in pressure resulted in the same sequence of 
phenomena as at atmospheric pressure. The 
corresponding potentials were decreased, as was 
to be expected, and the visual form was spread 
out (Fig. 8), emphasizing its glow discharge-like 
characteristics. The current-voltage curves again 
showed a steeper slope® for lower pressures due to 
increased mobilities (Fig. 9). 


POSITIVE AND NEGATIVE CORONA CURRENTS 


In pure gases transparent to photons it will in 
general be expected that near onset the negative 
corona currents would be several orders of 
magnitude greater than the positive corona cur- 
rents for the reason that with a negative point 
the high field is concentrated about the cathode. 
Hence the necessary secondary ionization due to 
any one of three agencies, photo-ionization in the 
gas, photoelectric effect at the cathode, secondary 
electron emission by positive ion bombardment, 
is active in a high field region and ion multipli- 
cation is effectively ensured. In the case of the 
positive point the main source of electrons avail- 
able for maintaining the discharge is the photo- 
electric effect at the distant plate. The low fields 
at the plate yield practically no electrons by 
positive ion bombardment and _ photoelectric 
emission from the plate is small. Hence, in the 
case of the positive point in Hz or No, both of 
which fulfill the condition of being transparent to 
photons, secondary electron emission is small, 
and the observed currents at onset are small. For 
the negative point in Hz and N, the positive ions 
are accelerated in the high field region toward the 
point without interruption by negative space 
charges and hit the point surface with energies 
10 to 100 times higher than the energies with 
which they would strike the negative plate during 
a positive point corona; therefore, it is assumed 
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that the second Townsend coefficient takes on 
relatively large values, thus yielding high cur- 
rents (Figs. 4 and 9, also shown by current values 
in Tables I and IT). 


CLEAN-UP CORONA 


Clean-up effects on He contaminated with O, 
had been noted earlier. Bennett? ascribed them to 
the formation of an oxide of tungsten. The points 
used in our work were chiefly Pt but for test 
purposes an auxiliary W point was mounted in 
the chamber. 

1. In the point-to-plane gap used in these 
measurements the point was of cylindrical shape 
with a hemispherical cap 0.05 cm in diameter, 
while in Bennett’s work points of much smaller 
radius of curvature were employed, which natu- 
rally could not be shaped as accurately as the 
larger points. 

2. The auxiliary W point was used so that the 
carefully shaped Pt point was not changed in 
form by the cleansing discharge, a procedure 
which was omitted by Bennett. Only under such 
carefully controlled condition can the observed 
current-voltage characteristic be completely as- 
cribed to a change in the gas composition. For it 
follows both theoretically and from Bennett's 
observations that the Trichel mechanism which 
was used to detect electron attaching impurities 
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Fic. 9. Current-voltage characteristics for negative 
corona in pure N». The low potential end of the curves 
indicates “‘offset;’’ the high potential end does not indicate 
breakdown. Gap length 4 cm. I: p= 100 mm Hg, II: p= 300 
mm Hg, III: p=760 mm Hg. 
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is not active on fine points. At the same time a 
severe etching of the fine points of Bennett was 
observed in the clean-up run. Thus it could not 
be determined whether the disappearance of the 

Trichel pulses was due to clean-up or to etching. 

The procedure in this study was as follows: 
The chamber was carefully cleaned and filled 
with the pure gas, either N2 or He. The phe- 
nomena occurring from corona-onset to spark 
breakdown were then observed in the main gap 
on the Pt point to serve as a check on previous 
results and as a control for purity. Then a small 
amount, around 0.4 percent, of O2 was introduced 
into the chamber and the phenomena were re- 
corded over the same potential interval as before. 
After this the clean-up discharge was run from 
the auxiliary W point for 8 to 12 hours at a high 
negative current. At the end of the run the con- 
dition of the gas was examined by a repetition of 
the test procedure with the clean Pt point. 

Admitting traces of O» of less than 0.1 percent 
into pure He, one obtained in the test run the 
characteristic phenomena with positive and nega- 
tive points detailed above. After the cleansing 
run with a W point no pre-onset streamers were 
observed for the positive corona and no Trichel 
pulses for the negative corona. Clean-up had 
obviously occurred and seemed to be complete. 

If the impurities exceed 0.1 percent in He, the 
negative corona currents before clean-up fell to 
very low values due to Trichel pulse formation 
and the positive corona revealed streamers. 
Cleaning the gas with the discharge on the auxil- 
iary W point restored in a large measure the 
purity of the gas; the currents, however, never 
reached the initial values in the pure gas, and 
incipient streamers were noted in the positive 
corona. It was clear that some clean-up had 
occurred, but with large amounts of Oz: it was 
never complete (Fig. 5). Use of a Pt point to 
produce the clean-up did not materially change 
the results of the clean-up compared with the W 
point, which seemed to indicate that the forma- 
tion of tungsten oxide was not the primary 
mechanism for the clean-up action. This was 
substantiated by repeating experiments in N»2 
which was contaminated by Oz. Under these con- 
ditions no clean-up was detected (Fig. 6). There- 
fore, clean-up in Hz may be ascribed primarily 
to the formation of water. If the concentration of 
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O: is low, the water formed by the discharge will 
be completely adsorbed by the baked-out and 
very hygroscopic walls of the Pyrex chamber, 
Above 0.1 percent of O2. the amount of water 
vapor generated was too great to be completely 
adsorbed. Hence, one may ascribe the residual 
pulse formation to the remaining water vapor. 


ETCHING OF THE POINT ELECTRODE 


The occurrence of etching in pure nonreactive 
or inert gases about a negative point, as noted by 
Bennett, seemed peculiar owing to the high pres- 
sures used in these experiments and owing to the 
fact that Hudson™ found pitting always associ- 
ated with oxidation. While the phenomenon of 
sputtering occurs at the cathode of many dis- 
charge tubes, the sputtering becomes important 
only with higher cathode falls and low pressures. 
It is clearly a function of the energy of the 
incoming ions. In our experiments the positive 
ions strike the cathode with an energy not much 
in excess of one electron volt. 

An attempt was made to correlate the degree 
of etching with some properties of the metals, and 
point electrodes of W, Pt, Cu, Al, and Pb were 
tested on the negative corona in gases Hz and N, 
under standardized conditions (12-hour run, con- 
stant current of 1250 wa corresponding to a cur- 
rent density of about 1 amp./cm?, gap length 
3.1 cm). The points were carefully polished with 
tin-oxide and examined under the microscope. 
Etching occurred in all cases except Al and 
manifested itself in a symmetrical crater-like pit 
with a depression below the level of the point 
surface surrounded by a heaping up of metal at 
the sides. This heaping up may be due to the fact 
that a considerable fraction of the sputtered 
particles were reflected back and redeposited in 
the neighborhood of the region of impact owing 
to the high pressure of the gas molecules present. 
The amount of displaced material could be 
estimated roughly and appeared to be greater in 
N.2 than in He. No discoloration or other indi- 
cations of chemical activity were observed. No 
simple correlation between any properties of the 
metals and the number of displaced atoms could 
be made. Etching in Al did not result in a crater- 
like pit, instead the whole surface of the tip of 
the point was roughened by a large number of fine 
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holes which were caused by a discharge constantly 
changing its position of the point because of 
alterations of the oxide-coated surface film by the 
impact of positive ions. 

The introduction of impurities made the in- 
vestigations difficult because it immediately re- 
duced the currents by a factor of more than one 
hundred, thus making the time of observation 
prohibitively long. There was no doubt, however, 
that etching also occurred in these cases. 

While not many data as to the mechanism are 
at hand, the character of the process is sufficiently 
like the sputtering phenomenon in glow dis- 
charges to make one feel that one is here dealing 
with a type of sputtering occurring at relatively 








low positive ion energies, but in distinction to 
sputtering caused by positive ions of high energies 
in low pressure glow discharges, this sputtering 
occurs with low ion energies at high current 
densities. 


ACKNOWLEDGMENTS 


The writer wishes to express his appreciation 
to Professor Leonard B. Loeb, who suggested the 
problem, guided its progress, and helped in the 
interpretation of the observed data. The writer is 
also indebted to the Research Corporation of New 
York, under whose Research Fellowship Grant at 
the University of California the main part of this 
work has been completed. 











PHYSICAL REVIEW 


VOLUME 63, NUMBERS 3 AND 4 


The Specific Primary Ionization of Cosmic Rays in Helium 


Wayne E. Hazen 


FEBRUARY 1 AND 15, 1943 





University of California, Berkeley, California 


(Received November 17, 1942) 


Values for the specific primary ionization of 46 cosmic rays in helium have been determined 


from cloud-chamber photographs. A magnetic field allowed separation into two groups: 21 
electrons with energies near the energy for minimum ionization had a mean ionization of 
7.33+0.12 per cm, and 26 mesotrons with energies greater than the energy for minimum 
ionization had a mean ionization of 7.23+0.12 per cm. The expected value for the latter 
figure, taking into account a logarithmic rise of ionization with energy, is 8.80 ions/cm. The 
experimental value for the minimum specific primary ionization in helium (corrected for the 
ionization in the vapor) is 6.5+0.1 at N.T.P. Three protons were observed in the Hp range in 
which they can easily be identified; in this Hp range three particles constituted approximately 


one percent of the radiation. 





HE primary ionization of electrons from 
radioactive sources has been determined by 
Williams and Terroux' in hydrogen and oxygen 
and by Skramstad and Loughridge’ in nitrogen 
and neon. C. T. R. Wilson* counted the primary 
ions along the tracks of electrons with energies 
between ten and thirty kev. Theory predicts a 
logarithmic increase in ionization with energy 
when the kinetic energy is large compared to the 
rest energy. Kunze‘ failed to detect any rise in 


1 E. J. Williams and F. R. Terroux, Proc. Roy. Soc. A126, 
289 (1930). 

*H. K. Skramstad and D. L. Loughridge, Phys. Rev. 
50, 677 (1936). 

°C. T. R. Wilson, Proc. Roy. Soc. A104, 192 (1923). 
*P. Kunze, Zeits. f. Physik 83, 1 (1933). 









ionization in measurements on ten cosmic-ray 
particles with energies around 10° ev. Since it is 
now known that most of the particles with this 
energy are mesotrons, the ionization should corre- 
spond to that of electrons with Hp=10°m/300xz, 
where uw and m are the masses of the mesotron and 
the electron, respectively. When this correction is 
made, Kunze’s value of 19 ions/cm in normal air 
agrees with the minimum ionization found for 
electrons. C. T. R. Wilson*® also counted primary 
ions in nine straight tracks of unknown energy 
which presumably were cosmic-ray mesotrons 
and obtained a value of 20 ions/cm. In the 
present experiment the primary ionization of both 
cosmic-ray electrons and mesotrons was measured. 
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APPARATUS 


The cloud chamber used in obtaining the 
photographs was 30 cm in diameter and 30 cm 
deep with a magnetic field of about 1000 oersteds 
from two Helmholtz coils. Illumination was pro- 
vided by four 6000-lumen medium bi-post street 
lamps which-were operated at about 70 volts for 
one second. Four sets of four-inch condensing 
lenses were adjusted to give slightly convergent 
beams of light which entered the cloud-chamber 
cylinder near the rear and passed out of the 
cylinder near the front, making an angle of 50° 
with the axis. Photographs were taken on Agfa 
Ultraspeed 35-mm film with a stereoscopic camera 
equipped with 5-cm lenses and a curtain type 
shutter. The lenses were focused for different 
distances (54 and 59 cm) in order to make use of 
a large portion of the cloud-chamber volume. 
With lens apertures of f : 8 the depth of focus of 
each lens was about 5 cm. The program of opera- 
tions involved in taking photographs was con- 
trolled by a cam system driven by an induction 
motor, the speed of which was essentially 1800 
r.p.m. A pulse of x-radiation could be sent into 
the chamber at any desired time during the 
program in order to see if the chamber was still 
sensitive. The pulse came from a Coolidge type 
tube during the application of a 15,000-volt a.c. 
pulse to the target. 


PROCEDURE 


The chamber was filled with helium gas and 
the saturated vapor from a 1:3 water-ethanol 
liquid mixture, to a total pressure of 1.2 atmos. 
Such a gas mixture gives about seven primary 
ions per centimeter of cosmic-ray path at the ex- 
panded pressure and since the magnification was 
about one-tenth, one obtains about seven droplet 
images per millimeter on the film. The mean 
image diameter was 0.003 cm so that the correc- 
tion factor for overlapping of images is 1.1. In 
using the cloud chamber to determine ionization 
one must consider the efficiency of ions as con- 
densation nuclei at the time the track in question 
was formed. In measuring primary ionization 
only post-expansion tracks can be used because a 
large cloud chamber cannot conveniently be ex- 
panded quickly enough to produce a truly sharp 
track. Expansions are therefore made at random 
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Fic. 1. Post expansion tracks of cosmic rays. There are 
two electron tracks in the left picture, a mesotron in the 
middle picture, and a proton and two electrons in the 
right picture. A pre-expansion diffuse track also appears in 
the right picture. The index marks represent 5 cm in the 
cloud chamber. 


and sharp tracks in good focus sometimes occur. 
The frequency of occurrence depends on the 
usable sensitive time of the chamber, which in 
turn may depend on the velocity of fall of the 
drops or on the actual sensitive time of the 
chamber. In the case of the present experiment, 
the maximum sensitive time for condensation on 
about 50 percent of the positive ions was 0.3 sec. 
The dependence of condensation efficiency on 
expansion ratio has been determined by Bagley‘ 
and by Nielsen* and the dependence of sensitive 
time on expansion ratio by the author.’ From 
such information the fraction of the sensitive 
time during which the condensation efficiency on 
ions remains greater than say 90 percent can be 
found. This fraction is approximately one-third 


5G. D. Bagley, Phys. Rev. 56, 851 (1939). 
®°C. E. Nielsen, Ph.D. thesis, University of California 


(1941). 
7W. E. Hazen, Rev. Sci. Inst. 13, 247 (1942). 
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and the usable sensitive time, when we allow ten 
percent for uncertainty, was therefore 0.1 sec. in 
the present experiment. The sensitive time was 
checked after every few expansions and the 
chamber was adjusted to keep the sensitive time 
greater than 0.3 sec. 

Since the drops required 0.15 sec. or more to 
grow to photographable size, under the condi- 
tions existing in this experiment, the camera 
shutter was tripped 0.25 sec. after the expansion. 
At the end of 0.25 sec. a water-alcohol drop in 
helium falls two cm/sec. and as a result the image 
moves 0.2 cm/sec. A shutter speed of 1/150 sec. 
therefore allows an image motion of only 0.0013 
cm (less than one-half the image diameter) and 
the resulting elongation is scarcely noticeable. 
With this shutter speed good photographic images 
were obtained with lens apertures of f:9; larger 
apertures increase the image size and hence lower 
the resolution. Tracks of particles that passed 
through the cloud chamber immediately after the 
expansion were not appreciably distorted by 
their fall of 3 mm which occurred before the 
photograph was taken. 

The position and angle of the tracks were de- 
termined by reprojection through the camera 
onto a screen which was adjusted so that the two 
images coincided. This method allowed the de- 
termination of track lengths to within a few 
percent. The apparent curvature of the images of 
straight lines that did not pass through the 
camera axis was found to be negligible (40-m 
radius of curvature). Reproductions of some of 
the photographs are shown in Fig. 1. 


TABLE I. Figures used in estimating mean primary 
ionization per cm for a group of cosmic rays with Hp> 10°. 
Blackett gives the spectrum in terms of energy E. Since 
the values for E were obtained from E=300Hp and since 
there are significant departures from this simple relation 
for mesotrons in the lower energy range, Blackett’s values 
of E were changed to Hp by dividing by 300. 


Fraction of 
particles 
with 
Hp >10° 


Products of 
columns one 
and two 


Fraction of 
sea level 
radiation 


Mean 


Hp group ionization 





1.64 
3.52 
2.42 
0.80 
0.42 
= 8.80 


0.213 7.7 
0.414 8.5 
0.260 9.3 
0.051 0.077 10.3 
0.024 0.036 11.7 
Expected mean ionization 


1—3X 10° 

0.3-—1X 107 

1—3X107 

0.3—1X 108 
> 105 


0.142 
0.276 
0.173 


IONIZATION OF 


COSMIC RAYS 


RESULTS 


Out of 180 photographs, 46 tracks of high 
velocity particles were obtained which were in 
sharp focus for at least several centimeters. These 
46 tracks fall into two distinct groups; one group 
consists of 21 electrons in the Hp range from 
2.5 10* to 3X10*, the other of 25 particles with 
Hp=10°. The electron group extends over an 
energy band which is roughly centered on the 
minimum and is sufficiently narrow to allow con- 
sideration of the ionization of all the particles as 
a group. 

The weighted mean for the specific primary 
ionization of the 21 electrons is 7.33 ions/cm 
with a probable error of the mean of 1.6 percent 
or a mean error of 11 percent for one observation. 
The mean statistical error 1/\/n for the indi- 
vidual observations averages 13 percent; the 
main error in the observations is therefore 
statistical. The weighted mean for the 25 par- 
ticles with Hp>10° is 7.23 ions/cm with a 
probable error of the mean of 1.7 percent or a 
mean error of 12.7 percent for one observation. 
Again the mean statistical error for individual 
observations averages 13.4 percent and the main 
error is therefore statistical. 

The close agreement between the specific 
ionization for the electrons, which are known to 
have velocities near that for minimum ionization, 
and the specific ionization for the group of 
particles with Hp2>10*, indicates either that 
nearly all of the latter particles had an [ip less 
than 6X 10° or that the primary ionization does 
not increase as rapidly with velocity as the 
theory predicts. The energy spectrum for parti- 
cles with Hp greater than 10°, as determined by 
Blackett,*® is given in Table I together with the 
ionization® predicted by assuming a minimum 


SE. J. Williams, Proc. Roy. Soc. A172, 206 (1939). 

* The value for the ionization (/) to be multiplied with 
the fractional number of particles in a given energy band 
can be approximated as follows for the higher energies. 
The fractional number of particles df with energy between 
E and E+dE is essentially df =cdE/E*, where c=constant. 
The ionization J of particles beyond the minimum is given 
approximately by J=k log E/Eo, where k and Epo are 
constants. Then [= /;*/df/ fi2df and if all the quantities 
are expressed in terms of E 

log B/Ey= £2108 £1/Eo~ E: log Ex/Eo, 5 
E.-—E, 


where k log E/Ey=T. For E:=nE; 


E=eE,/n'!\-» (e is the base of natural logarithms) 
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Fic. 2. Experimental values for the minimum primary 
ionization per centimeter in He, He, Ne, No, and Os are 
indicated by the circles (the upper circle at abscissa 2 is 
for He). The calculated value for atomic hydrogen is 
indicated by the cross. The straight line was used merely 
to estimate values for other molecules. 


ionization of 7.3 and by assuming that the pro- 
portional increase in ionization with Hp is the 
same as that for atomic hydrogen.'® Since the 
logarithmic rise obtained experimentally by 
Corson and Brode" for the probable ionization in 
nitrogen is about the same as the calculated value 
for atomic hydrogen, the second assumption just 
noted seems amply justified. By weighting each 
ionization with the corresponding fractional 
number of particles in that part of the spectrum, 
and by summing the resulting products, one 
obtains 8.80 ions/cm as the expected value for 
the mean ionization. This result is 22 percent 
larger than the observed value of 7.23 ions/cm, 
whereas the probable error of the mean is only 
1.7 percent. Fermi” has shown that polarization 
of the medium reduces the ionization, but the 
effect is negligible in air for energies less than 


and for n=3, E=1.6E; or forn= «, E=eE,. The particles 
must be considered in several energy bands because the 
energy spectrum does not follow the 1/E* distribution 
closely: the spectrum has a maximum at E= 10° and falls 
off less rapidly, than 1/E? for higher energies. 

10 E, J. Williams, Proc. Roy. Soc. A139, 163 (1933). 

11 R. B. Brode, Rev. Mod. Phys. 11, 222 (1939). 

12 E, Fermi, Phys. Rev. 57, 485 (1940). 
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10'° ev. There are so few particles with energy 
greater than 10'° ev that even if one assumes no 
increase whatsoever for the ionization in this 
range, the mean ionization should still be 8.70 
ions/cm. 

The above values are for the ionization in the 
helium, alcohol, and water atmosphere existing 
in the cloud chamber at a total pressure of 1.03 
atmos. and a temperature of 22°C. The ionization 
in the alcohol and water vapor can be estimated 
from their partial pressures and their electron 
densities. If one plots the values for the ionization 
in He, Ne, Ne, and Oz obtained by previous ex- 
periments’? as a function of the number of 
electrons per molecule the ionization of other 
elements can be estimated from the resulting 
curve (Fig. 2). Thus for C2H;OH the value for J 
is 33 or 42 ions/cm depending on whether one 
chooses to extrapolate to »=26 electrons per 
molecule or to add values of J for Cs, 3He, Ox. 
The latter method seems more reasonable and 
was adopted. The partial pressures of H2O and 
CsH,;OH in equilibrium with a 1:3 liquid mix- 
ture were obtained from Gautier and Ruark." 
The final value for He at N.T.P. is 6.5 ions/cm. 

Three of the observed tracks had ionization 
and momentum consistent with protonic mass. 
Two of these were diffuse and appreciably curved 
in the magnetic field; the third (Fig. 1) was 
straight on each side of a bend produced by a 
close nuclear encounter but a delta-ray aided in 
the identification. These three particles consti- 
tuted approximately one percent of the radiation 
in the Hp band from 10® to 3X 10°. 

The magnet used in the investigation was 
designed and assembled by William Fretter. A 
detailed description will probably appear else- 
where at some later date. In concluding, the 
author wishes to express his appreciation to 
Professor R. B. Brode for the suggestion of the 
problem and for advice in carrying out the 
research. 





1% T, N. Gautier and A. E. Ruark, Phys. Rev. 57, 1040 
(1940). 
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The specular reflection of plane waves is treated by the 
wave theory method. For the case in which the velocity of 
wave propagation and the radiation impedance vary suffi- 
ciently slowly that multiple reflection may be neglected the 
reflection coefficient is first obtained for monochromatic 
waves incident on a set of m plane parallel plates with 
otherwise arbitrary properties. The method may also be 
used to obtain the reflection coefficient in the case of large 
reflectivity but the results are exceedingly cumbersome for 
a large number of plates. The result for monochromatic 
waves is used to obtain the reflection coefficient for a plane 
wave pulse of arbitrary shape and duration. It is shown 
that the reflected wave train consists of a series of pulses 
of exactly the same form as the incident pulse with the 
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amplitude of each reduced by the value of the appropriate 
Fresnel reflection coefficient. The interference effects are 
also represented. The sufficient condition for reproduction 
of the form of the incident pulse in the reflected wave is 
that the reflection coefficient for monochromatic waves be 
expressible as a sum of terms each of which is periodic in 
the frequency. In contrast to previous statements in the 
literature, periodicity of the coefficient itself is unnecessary. 
The results for the m plate problem are generalized for 
both monochromatic waves and pulses for the case of 
stratified media having reflective properties which vary 
continuously in any manner. The special case of rectangu- 
larly shaped pulses is analyzed and application is made 
to the reflection in sea water and in the atmosphere. 





I, INTRODUCTION 


T is a well-known fact that the propagation of 
wave disturbances is characterized by reflec- 
tion wherever there is a change in the physical 
properties of the medium. In the case that these 
changes are discontinuities localized at simple 
surfaces it is a very simple matter to determine 
the fraction of the incident energy which is 
reflected. Thus, the reflection coefficient for 
monochromatic plane waves incident obliquely 
on a plane parallel plate of homogeneous ma- 
terial separating two media with refractive in- 
dices differing from each other and from that of 
the plate is to be found in most standard texts.'? 
The same problem has been solved for the case 
of plane wave pulses of arbitrary shape by 
Muskat who showed that the reflected train of 
waves will consist of a series of pulses of the same 
shape as the incident one but, of course, with 
diminished amplitude.? The problem of reflection 
* Now at Princeton University. 


1For example, Rayleigh, Theory of Sound (Macmillan 
_ Company, New York, 1929) second edition, vol. 2, p. 


*M. Muskat, J. App. Phys. 9, 275 (1938). Each pulse in 
the reflected train corresponds to an odd number of internal 
reflections in the plate and the amplitude of each individual 
multiply-reflected pulse contains the well-known Fresnel 
reflection coefficient raised to an appropriate power. 
Although the pulses retain their form this would be recog- 
nizable in the reflected wave pattern only if there is no 
overlapping, that is, if the time for traversing the plate 
twice along the directions of the incident and reflected 
waves is less than the duration of the initial pulse. 
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of plane waves, whether monochromatic or not, 
in a medium of continuously varying refractive 
index (and/or density) does not appear to exist 
in the literature. 

Some examples for which the treatment of this 
problem would be of interest may be mentioned. 
Thus, for the propagation of sound signals (or 
electromagnetic signals) through the atmosphere 
the properties of the medium on which the reflec- 
tion depends are, in general, continuously vari- 
able. A similar situation applies in the case of the 
propagation of sound in sea water where the 
acoustic velocity depends on hydrostatic pres- 
sure, salinity, and temperature and hence changes 
continuously with depth. In the two cases cited 
the reflected intensity is, in general, partly due to 
convection currents which are usually present 
but in order to understand the role of each 
mechanism of reflection it is advantageous to 
consider each separately. We therefore leave out 
of consideration any movement of the medium. 
In the following we shall first consider the re- 
flection of monochromatic plane waves in a 
medium of arbitrarily varying reflective proper- 
ties. We may then proceed to the treatment of 
plane wave pulses of arbitrary shape and du- 
ration.® 


* It is perhaps worthy of mention that Rayleigh, in his 
book (reference 1) on page 86 gives the following reference 
which is,here quoted verbatim: “It is interesting and 
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II. REFLECTION OF MONOCHROMATIC 
PLANE WAVES 


Although the following analysis is equally 
applicable to sound waves in fluid media, trans- 
verse polarized waves in general elastic media 
and to polarized electromagnetic waves, we shall 
adopt the language of acoustics. The appropriate 
modifications in terminology for application to 
the other types of waves will be readily apparent 
and needs no further comment at this time. 

Concerning the medium in which the plane 
waves are propagated we make the following 
assumptions: 


(a) The medium is stationary, that is, there 
are no convection currents, nor thermal 
currents, etc. 

(b) The medium is isotropic. This is obviously 
the case in those media for which the 
present problem is pertinent. 

(c) The medium is stratified. That is, all 
physical properties upon which reflection 
depends (velocity and radiation impedance 
mainly) are constant over a set of parallel 
planes and vary only in the direction 
normal to these planes. 

(d) The medium is non-dissipative. This as- 
sumption is particularly well justified in 
the acoustic case or for the longer electro- 
magnetic radiations. 

(e) The medium is non-dispersive. This as- 
sumption is, of course, of particular im- 
portance for the more general case of 
reflection of pulses considered in Section 
III below. 


In accordance with assumption (c) the velocity 
of sound c, and the density and compressibility 
of the medium as well, depends on only one 
variable x. In the planes in which the velocity, 
etc., are constant we use a rectangular coordinate 
system y, 2 and we choose the y axis in the plane 
of incidence. The problem is then two-dimen- 


encouraging to note Laplace’s remark in a correspondence 


with T. Young. The oat analyst writes (1817) ‘Je persiste 
a croire que le probléme de la propagation des ondes, 


lorsquélles traversent différens milieux, n’a jamais été 
résolu, et qu’il surpasse peut-étre les forces actuelles de 
l’analyse.’ (Young's Works, vol. I, p. 374).”” On the other 
hand, quoting from Rayleigh, page 78 of the same reference: 
“The problem of vibrations of a variable medium is 
probably quite beyond the grasp of our present mathe- 
matics’’! 
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sional and 2z is a non-essential coordinate. With- 
out loss of generality we may suppose that there 
exists a region of uniform properties in which 
the incident wave travels. The reflecting medium 
may then extend from x=0 to x=a and we sup- 
pose that for x>a the medium is again uniform 
and non-reflecting. As will be evident, the case in 
which a discontinuity in reflective properties 
occurs at any point is readily obtained from the 
analysis to be given. 

The wave equation for the velocity potential 
may then be set up with the physical constants 
which enter it being regarded as arbitrary func- 
tions of x. In order to solve this equation for 
the most general case we adopt the following 
procedure.‘ The reflecting region between x=0 
and x=a is divided into m subregions by n+1 
planes at x=x,, v=0, 1---m. In each subregion or 
lamina, whose thickness is 6,=x,—%x,-1, the 
velocity c (also the density p) have arbitrary but 
constant values. We thus consider the case of 
reflection from m plane parallel plates of homo- 
geneous material. Since this problem possesses 
some intrinsic interest we may treat the general 
case in which the plate thicknesses 6, are not 
necessarily equal. This generalization is, of 
course, unnecessary for the transition to the case 
of continuous media which will be made later. 

We denote the region of the vth plate by the 
subscript v. The homogeneous medium in which 
the wave is incident (x<0) is then denoted by 
the subscript 0 and the homogeneous medium 
for x>a by the subscript +1. Then if ¢, is the 
velocity potential in the vth region (x,-1 <x <x,) 
and ¢, is the corresponding velocity, the wave 
equation is 

i , f - f _ 
Oey, 
dx Oy? ¢? Of 





where we have used the time dependence e~ 


‘An alternative method for the solution of the wave 
equation with arbitrarily given functions c(x), etc., would 
be the approximate W.K.B. procedure, see for example 
G. Wentzel, Zeits. f. Physik 38, 518 (1926). Here one 
assumes that the relative change in velocity, etc., in a 
distance of one wave-length d is small; i.e., Ad log c/dx<1, 
etc. Although this assumption is justified in most cases of 
interest and the correct order of magnitude of the reflection 
coefficient is obtained, the W.K.B. method neglects inter- 
ference effects, which are in general important. The method 
used here, on the other hand, is exact and automatically 
takes the interference effects into account in a proper 
manner. 
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for monochromatic waves, and 


k,=w/c,, (w=2xeX frequency). (2) 


In each region the velocity potential will consist 
of two terms corresponding to a wave traveling 
in the direction of increasing x (incident wave), 
and a reflected wave traveling toward decreas- 
ing x. Denoting the angle between the direction 
of propagation and the positive x axis by 3, and 
3,’ for these two waves, respectively, we have 
aside from the time factor 


¢,=B, exp [ik,(x cos 3,4 sin 3,) ] 
+B,’ exp [tk,(x cos 3,’+y sin 3,’)], (3) 


where B, and B,’ are amplitude constants which, 
in general, are complex. For x>a we have no 
further reflection ; therefore 


Biyi=0. (4) 


Since we are interested in obtaining the reflection 
coefficient (for the amplitudes) which is given by 


the ratio 
r= By'/Bo, (5) 


we may, without loss of generality, take Byo=1. 
Then r=By’, and the fraction of the incident 
energy flux which is reflected is 


R= |r|*= | Bo'|?. (6) 


There remains 2n+2 unknown cohstants, Bo’, Bi, 
By’, ---, Bn, B,’ and By. which are to be de- 
termined by the boundary conditions: 

(a) Continuity of the normal component of 
velocity 0¢/dx and 

(b) Continuity of the pressure pd¢/dt or of pd, 
at each of the »+1 boundaries. These conditions 
are 2n+2 in number and are therefore just 
sufficient to determine all the constants, and in 
particular Bo’. 

In order that the boundary conditions which 
are to be applied at x=x, be fulfilled at all 
points of the boundary it is evident that the 
exponents 


k,(x, cos 3,+- sin 8,), 
ky41(x, cos 3,41 +Y sin 3,1), 


k,(x, cos 3,’+y sin d,’), 


ky41(x, Cos ity sin 3,41), 


be the same for all y. Hence 
k, sin 3, =k, sin 3,’ =k,4; sin 8,41 

=hp41 Sin 8541, (7) 
d,=r—0,', (7a) 
sin 3,/sin 3,41 = k,41/k, =C,/C,41. (7b) 


Equation (7a) gives the specular law of reflection 
and (7b) gives the usual Snell law for refraction. 
The latter relationship, of course, defines the 
path of a ray in the medium for each angle of 
incidence. 

We may now drop the common factors in (3) 
and apply the boundary conditions to 


or ’ . P 
sin 3,=sin #,’, 
and 


B, exp (ik,x cos 3,) +B,’ exp (—ik,x cos #,). 
Applying these conditions at x, and introducing 


B, —= prB,, B,’ = p,B,’, 


prky cos v, 
= = t,=k,x, cos v,, 





%H+1= 
Pritk, COS 8, ’ 


” = h, > v, ’ 8 
—_e Nr+1 +1X%, COS Vy+1 (8) 


8, exp (tt,) +8,’ exp (—7é,) 
=B41 exp (imy1) +641 exp (—in41), (9a) 
8, exp (it,) —,’ exp (—i&,) 
= ar+i[Br41 Exp (in,+1) 
—Bi41 exp (—in41)], (9b) 
v=0,1---n. 
These equations constitute recurrence formulae 


in the #’s whose definition is completed by the 
auxiliary conditions 


Bo=po, Bn4i=0. (9c) 
The reflection coefficient is now given by 


r= Bo' /Bo=Bo'/ po. 


The recurrence formulae may be placed in 
more convenient form by solving (9a) and (9b) 
for 8, and 8,’. We find 


B, =} exp (—48,)[(1+a41) exp (im41)Br41 
+(1—a,41) exp (—imy:)Be41], (11a) 

B,’ =} exp (<,)( (1 — Gy41) EXP (i741) Br+41 
+(1+-,41) exp (—im41)B,41]. (1b) 


(10) 
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Although these equations may be solved for 
any 8,’ or B, in terms of Bo=po by a straight- 
forward iterated process of substitution, it seems 
impossible to obtain a general formula for arbi- 
trary @ and an arbitrary number of laminae. 
Thus, the above equations may be used for the 
problem of a finite number of plates although 
the result for a large number of plates is ex- 
tremely cumbersome; but it is evident that for 
the case of a continuously varying a, where the 
number of plates must be taken to be very large, 
it is necessary to make some simplification in 
order to obtain a practical solution. We now 
make the simplifying assumption, which will be 
generally well fulfilled, that all the a, are very 
nearly equal to unity. This means that the 
relative change of a [cf. Eq. (8) ], is small com- 
pared to a itself. Thus we write 


a=1+e and «<1. (12) 


The physical meaning of this approximation will 
be made clear in the following. 

It is to be noted that the magnitude of the re- 
flection coefficient, aside from interference effects, 
is determined by a (or better a—1=€), which 
contains the velocity and density (and the angle 
of incidence as well). The velocity and angle of 
incidence are also contained in the exponents & 
and 7 but these exponential factors control the 
interference phenomena only. With the approxi- 
mation (12) it is easy to see that the reflection 
coefficient Bo’/B> will be of order ¢ and all the 
other ratios 8,’/8, will also be of order e. There- 
fore we write (11a) and (11b) in the form 


8, =B,+1 exp (t5,), (12a) 
B,’ = 8,41 exp (—is,) —¥,, (12b) 
where 
Sy = M41— &, (12c) 
and 


Y= 36418941 exp [4(& +41) J. (12d) 
From (12a) we find® 
B,=Bo exp [—1(sotsit-+-+5,-1)], (13a) 


5 This approximation would not be sufficient for the 
determination of the transmitted intensity ~|8,4:|*. For 
this it is necessary to consider the correction terms of 
order ¢*. This may be seen directly from a development of 
the recurrence relations and it is to be expected from the 
conservation of energy flux since the reflected energy is 
of order e as compared to the incident energy, [see Eq. 
(16) below]. 





E. 
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so that 
Yo = 364180 exp [—4(Sotsit---+5,) 
+74(& +41) J. (13b) 
From (12b) we find 
Bo’ = —yo—1 exp (—t50) —y2 exp [—i(so+s:) ] 
—-+++—y, exp [—i(sotsit- ++ +5,-1)] 
+Bra1 exp [—7(sotsit--+++5,)]. 


Setting v= and using (9c) and (13b) we have 
for the reflection coefficient : 


r= —}$[e. exp [27(m1—50) ] 
+e: exp [2i(n2—so—si) ]+--- 


+ én41 exp (27(mn41—So—S1° “ —Sz)) ]. (14) 


This result can be put in more convenient form 
if we note that the exponent in the term with 
€,41 as a factor is 


2x» = 2(m+1—So—S1— + + * — Sy) 
=2Umi-S w+ fi). 
Since §5=0 [xo=0, see Eq. (8) ], we can write 
~=> (&—n)=2 k,6;cos 3;; xo=0. (15) 


The final result for the reflection coefficient for n 
plates in the ease that the density and com- 
pressibility (and therefore the velocity) change 
slowly is 

n+l 


r=—})> « exp (2tx,_1), (16) 


1 
where x, is given by (15) and «, is given by (8) 
and (12). That is, 


pl, COS d,s) p, tan d, 
“= (17) 





1 +641 = = 
Pv+1Cy+1 cos v, Pv+1 tan O41 


by (7b). 


This result for the reflection coefficient may be 


readily understood. Let us consider the case in 
which the first y—1 plates have identical proper- 
ties and the remaining plates are also identical 
but differ in their reflecting properties from the 
first group of plates. Thus, a discontinuity occurs 
only at the interface of the (v—1)st and the vth 
plates and the reflection takes place only at this 
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interface. The reflection coefficient is then 
Z,1—Z, 
Z,1+Z, 
where Z,-1, Z, are the radiation impedances for 


the media to the left and right of the reflecting 
plane. For oblique incidence the impedance is 


given by® 


’,= (16a) 


(16b) 


For our case in which the impedance changes 
slowly we have 


Zy-1 +Z, = 22-1 = 22». 


Hence, using (17), we have 


(\-)--4 
rt ~ — ee 
: Zont 


as in (16). Thus, aside from the interference 
effects contained in the exponential factors of 
(16) the result given for the reflection coefficient 
is simply the sum of the reflected amplitudes 
arising at each of the m+1 interfaces with the 
proper Fresnel coefficient, —$¢,. That the factor 
exp (2ix,-1) correctly represents the effect of 
interference may easily be seen when it is noted 
that 2x,-1 represents the change in phase of a 
wave reflected at the vth interface between 
plates y—1 and v. The factor 2 arises, of course, 
from the fact that the reflected wave is observed 
in the medium of incidence (medium 0) after 
propagation to and from the reflecting plane and 
the ray trajectories to and from this plane are 
symmetrical and give the same phase change. 
The significance of the approximation made 
above [Eq. (12) ] is now obvious. It is assumed 
that the impedance changes so slowly that mul- 
tiple reflections contribute a negligible amount 
to the amplitude of the reflected wave. The 
terms corresponding to multiple reflection will be 
of order ¢*, e--- corresponding to three, five, 
etc., reflections. The terms of order e* must evi- 
dently involve a triple summation since the three 
reflections may occur at any set of three planes.’ 


Z, = Prly SEC V,. 





(16c) 


6 See, for aa. J. C. Slater, Microwave Transmission 
(McGraw-Hill, New York), pp. 103-105. For the electro- 


magnetic case there considered we have the following cor- 
relations in notation: p=, magnetic permeability, 
c=1/(eu)!, € is the dielectric constant (cf. below). 

7 This is subject, of course, to the condition that the 
second plane lies between the first one and the medium 0 
-- hae third lies between the second and the medium 
n+1. 
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The terminology may now be generalized so 
as to include reflection of other types of radi- 
ation : 


(1) For the reflection of transverse waves 
polarized normal to the plane of incidence, 
1/p is the rigidity modulus and r is the 
ratio of displacement amplitudes in the 
reflected and incident waves. 

(2) For the reflection of electromagnetic waves 
polarized normal to the plane of incidence, 
p is the dielectric constant, and for polar- 
ization in the plane of incidence p is the 
magnetic permeability. In each case r is 
the ratio of magnetic or electric field 
vectors for the reflected and incident 
waves. 


In concluding our treatment of the reflection 
of monochromatic waves from m plates we con- 
sider the problem of non-reflecting systems. Such 
considerations are then applicable to the elimina- 
tion of reflection of light from glass or other 
materials by interference of the waves reflected 
from a set of plates, or the elimination of re- 
flection of longer electromagnetic waves in trans- 
mission lines. For these purposes it is advan- 
tageous to select the properties of the reflecting 
system so that the waves reflected from adjacent 
plates (for instance) are exactly out of phase. 
Thus, for the m plates, the thickness and re- 
fractive index of efch plate may be chosen so 
that a phase change of 2/2 takes place in travers- 
ing the plate once in the propagation direction, 
(quarter wave-length plates). For short waves it 
is convenient to use larger thicknesses and in 
this case the phase change may be made equal 
to x/2 plus any integral multiple of (integral 
multiple of half wave-lengths). This additional 
thickness does not change the interference effect 
and therefore gives the same reflection. Then for 
each plate 


q=k,6, cos 0,=-+mr, m=0, 1-- *» (18) 


which may be fulfilled for a particular angle of 
incidence. 

In this case it is possible to solve the recurrence 
relations (11) in a rigorous manner. We introduce 


F,=exp (t¢,)8,, F,’=exp (—ié,)B,’, (19) 
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so that (11a) and (11b) become 


F,=$((1+a41) exp (—igr41) Fr41 
; +(1—ay41) exp (4gr41) Fri] 
= (—)-t +41) R41 — (1 —ay41) Frat], (20a) 


F,’=3[(1 —a,41) EXp (1441) Fr4i 
; +(1+a,41) exp (iqr41) Fri] 
= (—)=C(1 — Qy41) Fui- (1 +ay+1) Fest], (20b) 


where we have used (18). The solution of these 
recurrence formulae gives for the reflection coeffi- 
cient : 
(1) The number of plates n is even: 
Oi2Ql4* * * An — A1A3* * * An+1 


f= ‘ (21a) 


G2Q4* * * An OL1G3* * *On+1 





In terms of the impedances of each plate [cf. 
Eq. (16b) ] this becomes 


Zo ZrZa- + -Zn)*—Zns(ZrZa: Zn)" 








r= .. (21b) 
Zo(Z2Z 4: + *Zn)?+Zn41(Z1Z3° + *Zn—-1)” 
(2) The number of plates m is odd: 
Ogl4* * *An+1— A1A3* * * Ay 
ae "Anz i tags: ** ay 
rs (2:23: - Zn)? —LoZnsi(Z2Za° ++ Zn-1)” (22) 





* (2123: Ln) +Lobnsi(ZrZe- ; -Zn-1)? 


Both (21a) and (22) reduce to the proper limiting 
value in the case of small «, [cf. Eq. (16) ]: 


r= —}(€1—€2-+e3— +++ (—)"€n41). 


It is obvious that there are many ways in 
which the impedances may be adjusted so as to 
give no reflection. One convenient way is to 
choose : 

(1) for m even 


A1y=A2, Ag=— Aq" **An-1 = An 


or, 
Z2=Zels, Ze=ZZq°+Zn-1=ZeZ0-2; (23a) 
(2) for n odd 
Qj= Aq, Ag=—Aa4°** An = An+1 
or, 


ZY=Lol2, 22=L2Z4-++Zn?=ZpiZny1. (23b) 
Thus, the impedance of each plate. is the geo- 
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metric mean of the impedances of the neighbor- 
ing plates. In the case that m is odd it is seen 
that the reflection coefficient is zero. In this way 
the reflection from each pair of adjacent inter- 
faces (at xXox1, X2X3, ***, Xn-1%,) vanishes by 
interference. In the case that is even, the values 
of the impedances given by (23a) gives a reflec- 
tion coefficient [from Eq. (21b) ] 


Zn— Zant 


oo, (24) 
“Z]atZnti 


This result is also obvious since there are an 
odd number of interfaces and the choice of 
impedances has been made so that the reflection 
from all but the last interface vanishes (in pairs) 
by interference. To obtain no reflection in this 
case it is only necessary to match the impedances 
Z, and Z,4; so that no reflection takes place at 
the last interface. 


III. REFLECTION OF PLANE WAVE PULSES. CASE 
OF N PLATES 


In his treatment of the reflection of a plane 
wave pulse from a single plate Muskat? has 
shown that if the reflection coefficient for the 
monochromatic case is exactly periodic in the 
frequency, the reflection from the plate of a 
pulse will consist of a series of wave trains of 
exactly the same form as the incident pulse. 
The reflection coefficient for two or more plates 
is not, in general, periodic in the frequency as is 
evident from the result (16) of the previous 
section. Nevertheless, it can be shown that the 
theorem demonstrated by Muskat still holds in 
our case. The periodicity condition is therefore 
sufficient but not necessary. 

The treatment of the reflection of the plane 
pulse for the case of m plates proceeds in just the 
same manner as for the case of a single plate. 
The incident pulse has an amplitude #; given by 


(25) 
®;=0; ¢>b and ¢<0O, 


where 


Here s; is the distance in the medium 0 measured 
along the wave normal, the time ¢ is measured 
from the moment of incidence and s; is measured 
from the point of incidence. The form of the 






f=t+s;/co. (25a) | 
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pulse f(¢) is entirely arbitrary as is the duration 
of the pulse b. 

The pulse may be regarded as a superposition 
of monochromatic plane waves given by the 
Fourier integral 


60) =—— f de f “Flweie Pd, (26) 


Since each component plane wave of frequency 
w is reflected with an amplitude r(w) times the 
corresponding incident amplitude, the reflected 
pulse is given by 


1; ab 
#,(5")=— f rede J fluye*-tdu, (27) 


where {’=t—s,/co and s, defines the phase planes 
in the reflected wave system. Inserting the re- 
flection coefficient r(w) from (16) we have 


1 n+1- b © 
@,(¢’)=—-— De duftu) { dw 
4n 1 Wy ae 


Xexp [iw(u—¢’+2p,-1)], (28) 


wherein the order of integration has been inter- 
changed. In (28) we have introduced 


v—l1 6; 
Psri=>- — cos 0; (28a) 
fef. Eq. (15) ]. 
The integral 


1 @ 

— f dwe'** = D(x), 
2r J _« 

which occurs in (28) has the following properties: 

D(x)=0 except for x=0 where D has a singu- 

larity of such strength that 


f Deas = D(x)dx =1. (29) 


From (29) it follows that 


+0 


f(x) D(x)dx = f(0) (29a) 

—0 
provided f(0) is non-singular. The limits of in- 
tegration in (29) and (29a) need only include 
the singular point. Now the reflected pulse may 
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be written 


n+l b 
$,=—} X of duf(u)D(u—¢’+2p,_1) 


n+1 


=—3 2 efi —2p,-1), (30) 
where f=0 unless 0 < ¢’ — 2,_1 <b. Equation (30) 
shows that the reflected pulse is given by a 
series of pulses each corresponding to the re- 
flection of the incident pulse at one of the n+1 
interfaces. Each component of this train of 
reflected pulses has an amplitude equal to the 
incident value multiplied by the Fresnel coeffi- 
cient — }«¢,. The interference is represented by the 
change in the argument (phase) of f and this 
change 2,_; is, as was explained before, just 
the expected phase change (divided by w) for 
reflection at the vth interface. Thus, the re- 
flected pulse consists of a series of wave trains 
identical in form with the incident pulse, aside 
from overlapping. The duration of each pulse in 
the reflected train is 6 and at s,=0 the vth pulse 
starts at t=2p,_,. Overlapping occurs whenever 


24, 
b> At, =2(p,— p,-1) = — cos #,. (31) 


Otherwise the pulses are separated by a time 
interval At,—b. It is obvious that for the case of 
continuously varying velocity and impedance, 
overlapping must occur and the form of the 
incident pulse will not be recognizable in the 
reflected wave pattern. 

The present results may be compared with the 
case of a single plate.? It may be noted that the 
assumption of small reflection (no multiple re- 
flection) corresponds in Muskat’s procedure to 
neglecting all terms in the Fourier series expan- 
sion of the monochromatic reflection coefficient 
beyond the second [cf. Eq. (7) of reference 2). 
The first term of this expansion (frequency in- 
dependent) corresponds to direct reflection at the 
first interface so that the wave does not enter 
the medium proper. The second term corresponds 
to a single reflection at the far side of the plate. 
It is easily verified that Muskat’s result for the 
reflected amplitude [Eq. (12) of reference 2] 
reduces to the first two terms of (30) in the 
limiting case of small reflection. 











IV. REFLECTION OF PLANE WAVES IN 
CONTINUOUS MEDIA 


The transition to the case of a medium with 
continuously varying velocity and impedance is 
readily made. For this we have only to consider 
the Fresnel reflection coefficient and the phase 
change 2),_1w. For the former we have 


p, tan 3,— pr41 tan 3,41 





€y41 >= 
Pv+1 tan O41 
In the limit of continuous variation this becomes 
(with 6,-dx) 


d 
€41—€(X) = —ée—- log p tan #. (32) 
ne 


For the phase change we have 
= dx 

pr ro(e)={ cosd. (33) 
e ¢€ 


For monochromatic waves the reflection coeffi- 
cient now becomes 


° d 
raf ae(= log p tan a eee, (34) 
0 dx 


where a is the thickness of the reflecting medium ; 
i.e., 


a= &,. (34a) 


For plane wave pulses the reflected wave is 
from (30) 


d 
&,(¢’) =3 fa(< log p tan ») 
0 dx 


Xf(¢’—2e(x)). (35) 


It may be noted that while the departure of the 
radiation impedance from homogeneity is essen- 
tial for the magnitude of the reflection this in- 
homogeneity is not essential in the interference 
factor g(x). In fact, to take the variation in 
velocity into account in evaluating ¢ would lead 
to correction terms of higher order in the param- 
eter of smallness (¢). Therefore, in view of the 
approximation of small reflection already made 
it is consistent and necessary to evaluate ¢ for 
a homogeneous medium. Then the reflected wave 
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(35) becomes 


* d 
$,=} { ax(— log p tan ») 
0 dx 
2x 
xs( = cos 0»). (36) 
Co 


This constitutes our final result. The reflection 
then depends only on the radiation impedance. 
It is evident that when the variation in im- 
pedance is given and the form and duration of 
the initial pulse is known it is a fairly simple 
matter to obtain the reflected wave pattern 
from (36). 

A particularly interesting example to which 
these results may be readily applied is that of 
the rectangular pulse. Then for arguments be- 
tween 0 and b the pulse f is constant, and since 
this constant enters only as a scale factor it may 
be taken equal to unity. Then 


7 ad p2 tan 3 
&,=} { ‘dx— log p tan d=} log = : 
z dx 





pi tan 3; ( ) 


where x; and x2 are the values of x for which 
the argument of f in (36) has the values 6 and 0, 
respectively; p2, #2 and pi, 3: are the values of 
p and # at x2 and x, respectively. Since we may 
take s,=0, ¢’=t=0 and x22x,20. For other 
values of s, we need only make a shift of time 
scale to obtain the reflection pattern from that 
at s,=0. In order to evaluate the limits x; and x, 
as functions of the time it is necessary to distin- 
guish between two cases. Introducing 7 accord- 
ing to 
2a 


T=— COS Bo (38) 
Co 


we must consider separately the cases b>7 and 


b<r. 
In case (1), )>+7, the duration of the pulse is 


greater than the propagation time r. Then it is 
seen that the limits are: 


x1=0, x2=4cot sec Vo; 


x1=0, x.=a; (39a) 


for 0<t<r, 
r<t<b, 
b<t<b+r, x1:=4c0(t—b) secdo, x2=a. 
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Fic. 1. Limits x: and x2 [Eq. (39)] plotted against ¢. 


In case (2) b<7 we find 
for 0<t<d, 
b<t<r, 


X1=0, xX2=}Cot sec Vo; 
X1=4$¢9(t—b) sec do, 
X2=4}Cot sec Vo; 

r<t<b+r, x1=}¢09(t—b) secdo, Xx2=a. 


(39b) 


If the limits x, and x2 are plotted against ¢ in 
the same diagram the figure obtained is a paral- 
lelogram as shown in Figs. 1A and 1B, which 
correspond to cases (1) and (2), respectively. 
In Fig. 1 the pulse duration 6 is considered to be 
the same in each case but for case (2) the re- 
flecting layer is thicker. Similarly, for a fixed 
reflecting layer, case (1) corresponds to a pulse 
of longer duration than does case (2). In either 
case the parallelograms have the same base } 
and the same angle @. Finally, the duration of 
the reflected train is always r+). 

The meaning of the foregoing is evident. At 
any time ¢ the reflected wave is made up of 
contributions due to reflections at points extend- 
ing at first from 0 to x2, where x2 is the maximum 
depth to which the wave may penetrate in order 
that after reflection it reaches a detector at x=0, 
at time ¢t. As ¢ increases x2 obviously increases 
until it reaches the maximum value of a. Even- 
tually x; departs from 0 because for larger values 
of t reflection from the layer at x=0 reaches the 
detector too soon. This must always occur for 
t=b, the duration of the pulse. The only differ- 
ence between cases (1) and (2) is that in the 
former case the reflection from the deepest 
point x2=a is attained before x; departs from 0, 
while in the latter case the opposite is true. 
Therefore there is a time interval r<t<b in case 
(1) during which the entire reflecting medium is 
operative in contributing to the intensity at the 
detector at a given time. 

From the results obtained in (39a) and (39b) 
the reflected wave is obtained at once from (37): 
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Case (1) 


p2 tan v2 
for 0<t<r, &,=}3 log ———_;; 
po tan do 
pa tan 3, 
rT<t<b, ,=}3 log ————_ 
po tan Vo 
(40a) 
=constant; 
p2 tan db, 
b<t<b+r, ,=} log ———. 
pi tan 3; 
Case (2) 
po tan de 
for 0<i<b, &,=} lo 
po tan Bo 
p b 
b<t<r, ,=} log — *. (40b) 
p, tan v; 
pa tan J, 
tT<t<b+r7, ,=} log 
pi tan 3; 


where p; tan #; and p2 tan 32 are evaluated at 
X1=}Co(t—b) sec Jo, x2=4}cet sec Jo. (40c) 


It is apparent that the reflected amplitude 
vanishes at the beginning and at the end of the 
time interval during which the reflected pulses 
reach the detector, (x1 =x,att=Oand att=b+7). 
Hence the absolute value of the amplitude or the 
intensity must reach a maximum during this 
time interval. If the impedance is a monotonic 
function of depth this maximum will occur during 
the second time interval r<t<b in case (1). In 
case (2) the reflected amplitude must reach a 
maximum in the second interval (b<t<7) since 
from t=0 to b, @, must increase while during the 
interval r to b+7, &, must decrease. The maxi- 
mum in this case will not be as great as in case (1) 
however. When the impedance does not vary 
monotonically it is not possible to predict the 
time interval within which the maximum occurs 
but in any case this maximum cannot exceed the 
value 


(p tan d) max 


innit (p tan 9) min 


(41) 

It is of interest to consider a numerical ex- 
ample in order to obtain a quantitative idea of 
the amount of reflection taking place in various 
media. For purposes of simplicity we consider the 
case of normal incidence which will suffice to 
give an order of magnitude estimate for the re- 
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flection. For normal incidence the impedance is 
HH simply pc and in the above formula the factor 
i" tan 3 is omitted. The reflected wave is then 


{ given by 






,=}$ log (02:/pz;). (42) 


Consider the case of reflection in sea water in 
which the density is a function of depth alone. 
The temperature and salinity of sea water are 
fairly constant for small depths (up to about 500 
meters at certain latitudes and at certain seasons). 
The temperature and salinity then decrease com- 
paratively rapidly over a distance of order 1000 
meters and then remain roughly constant, with 
the temperature approaching a value in the 
neighborhood of 4°C. Since the relative change of 
salinity is not very great (from about 0.036 to 
0.035 from the surface to 3000 meters depth) we 
may consider only the temperature variation. 
Then the maximum reflected amplitude is given 
by 


















ie Ateliend. Cf 


where pmax is the density at great depths and pmin 
is the surface density. For a surface temperature 
of 15°C and a salinity content of 0.036 the surface 
density is 1.0255 g/cm’. For a temperature of 
5°C and the same salinity the density at great 
depths is 1.0286 g/cm*.* Hence 


1.0286 
1.0255 











=0.0015. 





(#,) max — ; log 






Therefore at most the reflected energy flux is 
2X10-* percent of the incident flux. We may 
conclude that unless much greater variations of 
density may occur, reflection in sea water is 
extremely small and provides a relatively un- 
important mechanism of absorption. 

A much more important case of reflection 
occurs in the atmosphere where relatively greater 
density variations take place. As a rough approxi- 
mation we may take the barometric formula for 
the variation of density with height x above the 
surface of the earth. 


p= poe, (44) 


where 1/u=8000 meters. Considering again the 
















8 International Critical Tables (McGraw-Hill, New 
York), first edition, vol. III, p. 100. 
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case of normal incidence we have for the reflected 
wave a formula of the type 


}, = $u(x1—X2). (45) 


Since the reflecting layer will be, in general, 
sufficiently large to make the propagation time 
longer than the pulse duration, we must use the 
results found above for case (2). In this case the 
maximum reflection occurs for b<t<7r where 
Xe—X, attains a (constant) maximum value. 


Then 
},= —tyucob. (45a) 


For a pulse duration measured in seconds, and 
with co= 330 meters per sec. we find 


&,= —10~%. (45b) 


Thus for 1-sec. pulses, 1 percent of the energy 
flux is reflected. 
The ratio of reflected energy to incident 
energy is 
E, 1 2,2 
—=-{ od¢=—— 
E; bdo 16 
from (45) and (39b), and when we use (38) for 
3)=0. Since r>b we may write 


E,/Ey= wer? /64 = ye (ua)*. (46a) 


Equation (46) is valid only if uwa<1. For a= 2000 
meters (which implies a pulse duration )<12 
sec.), E,/E;~1/256~0.4 percent. Since we have 
used the barometric formula which underesti- 
mates the variation of density with elevation, 
the reflection should be somewhat more intense 
than is here indicated. 

We may note that the reflection of mono- 
chromatic waves in the isothermal atmosphere, 
the density being given by (44), is also quite 
small. From (34) we find (since the sound velocity 
is constant and therefore 3 is constant), 


i7?+3br—3b*) (46) 





tb 
r= —--— sec Vo. (47) 

0 
For a frequency of 1000 sec.~', |r| =1.6x10-* 
sec 3». This small reflection is essentially due to 
the slow rate of variation of the density. The 
distance 1/y, in which the density changes 
appreciably, is much larger than any other 
length involved in the problem. 
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The lead tetra-alkyls, Pb(CHs;), and Pb(C2Hs)«, were decomposed to yield free methyl and 


ethyl radicals in a specially designed furnace built into the ionization chamber of a 180° mass 
spectrometer tube. From the initial breaks of the ionization efficiency curves of the ions, 
CH;* and C:Hs*, we have found the ionization potentials of the methyl and ethyl radicals to be 
Tvert(CHs) = 10.09+0.1 ev and J yert(C2Hs) = 8.69+0.1 ev. These values are in excellent agree- 
ment with those found by indirect calculations from other electron impact data on hydro- 


carbons. Some of the implications of these results are briefly discussed. 





INTRODUCTION 


ESIDES their importance to the formulation 

of the theory of the thermal decomposition 
of hydrocarbons,' the energies of the methyl and 
ethyl free radicals have considerable theoretical 
significance. The calculation of the relative en- 
ergies of CHy, CHs +H, CH2+2H, CH+3H, and 
C+4H is a very suitable test of any quantum 
mechanical theory of valence.2 Thus accurate 
values of these energies are desirable for checking 
the calculations or for the evaluation of param- 
eters in the theoretical equations. A knowledge of 
the ionization potentials of the various radicals 
is important for the same reasons. 

In the recent past a considerable body of in- 
formation has been built up pertaining to the 
ionization and dissociation of various hydro- 
carbons by electron impact. In order to interpret 
these data in an unambiguous fashion and to 
check those interpretations which have been 
made, the energetic quantities mentioned above 
were needed. 

Several years ago Fraser and Jewitt’ reported 
values for the ionization potentials of the methyl 
and ethyl radicals, determined by the molecular 
beam technique. Their well conceived experi- 
ments suggested to us the desirability of re- 
peating this work with the mass spectrometer as 


* Presented at the Physical Society Meeting at State 
College, Pennsylvania, June, 1942. 
Westinghouse Research Fellow. Present address—Shell 
ayx.y = nee, Emeryville, California. 
1F, QO. Rice and K. K. Rice, The Aliphatic Free Radicals 


(Johns om Press, Baltimore, 1935). 

*H. H. Voge, J. Chem. Phys. 4, 584 (1936). The most 
elaborate quantal calculation. 

*R.G. ). Fraser and T. N. Jewitt, Phys. Rev. 50, 1091 
(1936); Proc. Roy. Soc. A160, 563 (1937). 
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a detector of the ionization. The repetition 
seemed particularly desirable, inasmuch as a 
great deal of the more recently obtained data 
were very difficult to reconcile with the ionization 
potentials reported by Fraser and Jewitt. 


EXPERIMENTAL 


The mass spectrometer used in this research 
is to be described by one of the authors (J. A. H.) 
in a separate publication, hence only pertinent 
details will be reported here, 

The electrons are obtained from an oxide- 
coated platinum cathode, their accelerating po- 
tential being controlled by means of a wire 
wound drum potentiometer. The accelerating 
voltage is read with an accuracy of +0.1 percent 
on a calibrated voltmeter. A separate pumping 








Fic. 1. The furnace and ion source assemblage. The lead 
alkyls were admitted to the tube through the quartz tube 
—the end of which is heated by the tungsten spiral. 
Emerging from the quartz tube, some of the molecules and 
free radicals are ionized by the electron beam. The ions 
are drawn to 2 by a small field between 1 and 2; passing 
through 2, they are accelerated through the slit in 3 by a 
difference of potential of about 500 volts between 2 and 3, 
entering the analyzer region where the magnetic field (H) 
sorts out the ions of different masses (180°-analyzer). 
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Fic. 2. The effect of furnace temperature on the relative 
abundance of the ions Pb(CH;)* (open circles), CHs* 
(closed circles) and CH,* (half closed circles) in the 
spectrum obtained from Pb(CHs),. 


lead to the arm of the tube, containing the 
cathode chamber, assures the complete removal 
of any pyrolysis products formed on the filament 
itself and on the adjacent parts of the tube 
heated by radiation from the filament. Argon 
admitted simultaneously with the substance 
being studied serves to calibrate the electron 
accelerating voltage scale. 

An isometric view of the ion source and the 
furnace in which the lead alkyls were thermally 
dissociated is shown in Fig. 1. The furnace con- 
sists of a quartz capillary tube (I.D.~1 mm) 
closely wound with 20-mil tungsten wire for the 
last 1.5 cm of its length. The wound zone is sur- 
rounded by a nichrome V radiation shield. This 
radiation shield, welded to one of the tungsten 
heater leads (80 mil) is connected to the high 
voltage electrode 1. The heater current is con- 
trolled by a variac, isolated from the high voltage 
electrodes by means of an insulating transformer. 
We have no means of measuring the temperature 
in the present arrangement, but five amperes 
bring the furnace to incipient red heat (~ 500°C). 
The end of the quartz tube is about 5 mm from 
the electron beam, which is indicated in the 
figure. 
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: Fic. 3. The variation of the relative abundance of the 
ions C2H,* (A), C:H5* (B), Pb(C2Hs)3* (C), and C2H¢* (D) 
in the spectrum of Pb(C:H;), with furnace temperature. 
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Fic. 4. Ionization efficiency curves: J and JI are CH;*, 
II, C2Hs*, and JV, C2H,*. See text, Table I. 


The lead alkyls were admitted at low pressure 
to the lead connecting to the furnace by means 
of a fine capillary leak, the pressure in back of 
the leak being adjusted by fixing the temperature 
of the liquid (Pb(CHs3)4~—10°C-Pb(C:Hs), 
~+30°C). The argon was admitted simultane- 
ously through a parallel leak. The vapor (and 
argon) flow through the quartz furnace into the 
path of the electron beam where the ionization 
takes place. The products of the thermal dissoci- 
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ation in the furnace can then be studied with the 
mass spectrometer. 

The specially purified samples of the lead 
alkyls were kindly supplied to us by the Ethyl 
Gasoline Corporation. In Fig. 2 the relative 
currents of the positive ions Pb(CH3)3+, CH,*, 
and CH;* produced from Pb(CHs),4 with 20-volt 
electrons are shown as a function of the heating 
current. The Pb(CHs3)3* is used as an indicator 
of the completeness of the thermal decomposition 
because the parent ion is absent. The striking 
feature of these curves is the sharp resonance 
found for the production of the methyl radicals 
as a function of the temperature. It should be 
noted that, as the temperature is increased 
beyond that shown in the figure, the CH;* 
current passes through a minimum and then 
increases at the same rate as the CH,* current; 
this indicates that in the higher temperature 
range the CH;* observed is the product of the 
dissociation of methane by electron impact. In 
addition to the ions mentioned above, there is a 
very large current at mass 28, increasing with 
temperature, which can be attributed to ethylene 
with confidence. 

At first it was hoped that a sufficient number 
of free radicals would be formed by dissociation 
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of the lead compounds on a filament close to 
the electron beam. The failure to observe any 
free radical formation by this technique with 
either tungsten or platinum filaments was ex- 
plained by the sharp maximum for the CH;* 
current as obtained with the final arrangement. 

As may be seen from Fig. 3, the behavior of 
lead tetraethyl is similar to that of lead tetra- 
methyl except for the lack of a pronounced reso- 
tiance maximum in the current of the radical ion. 

From the observations reported in the pre- 
ceding paragraphs and the fact that the pressure 
of the lead alkyls in the furnace was very low 
(about 10-° mm) it is clear that these compounds 
decompose on a hot quartz or tungsten surface 
to give primarily stable hydrocarbons, e.g., 
methane, ethylene, ethane, etc., at all tempera- 
tures. At first sight this result may seem to 
contradict the results of Paneth type experiments 
on the preparation and properties of free radicals." 
The latter experiments are carried out at much 
higher pressures and the decompositions are 
probably homogeneous thermal gas reactions. 
Our results are those of purely heterogeneous 
dissociation of the lead alkyls, and may well 
depend on the past history of the quartz or 
tungsten surfaces involved. 


Curves I and IJ of Fig. 4 are the ionization efficiency curves for mass 15 (CH;*) obtained with 
3.5 and 5.5 amperes through the furnace. Comparison with the curves of Fig. 2 indicates the following 


I. CH;-CH;t+e 
IT. CHyo>CH;t+H+e 


assignment of reactions to these ionization efficiency curves: 





A1(CH;*) = 10.09+0.1 ev‘ 
Ar(CH3*) = 14.44+0.1 ev. 


















The numerical values of the appearance potentials are based on the value J*(A) = 15.69 ev for the 
ionization potential of argon.5 This value is dependent on the old value of the conversion factor 
relating cm to ev. If the new value,® 1 ev=8066 cm-, is accepted J*(A) = 15.76 ev and the appear- 
ance potentials should be increased 0.07 ev. Since the dissociation energies which we deduce from 
electron impact data involve differences between appearance potentials, the error in the absolute 
voltage scale cancels. The appearance potentials were taken from the initial breaks of the curves.’ 
The appearance potential A; is to be associated with the vertical ionization potential of the methyl] 










* The uncertainties indicated for the appearance potentials in this paper are meant to show the reproducibility and 
the accuracy of differences. The range of voltages within which the absolute values of the appearance potentials lie is 
aie about twice as great as the indicated uncertainty, e.g., Az(CH3*) = 10.0+0.2 ev on an absolute scale. 

*R. F. Bacher and S. Goudsmit, Atomic Energy States (McGraw-Hill, New York, 1932). 

°G. Herzberg, Molecular Spectra I (Prentice-Hall, New York, 1939). 
7D. P. Stevenson and J. A. Hipple, J. Am. Chem. Soc. 64, 1588 (1942). 
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TaBLeE I. Summary of results. 




























Furnace current Fig. 4 Ion A(X*) Reaction 
PRICE: ‘4 3.5 amp. I CH;* 10.09+0.1 ev CH;~CH;t+e— 
Pb(CHs), 5.5 amp. II CH;* 14.4,+0.1 CH,~CH;*+H+e- 
POH 4.0 amp. III C:H;* 8.69+0.1 C2H;—~>C2Hs*+ «— 
Pb C:Hs). ; 4.0 amp. IV C:H,t 10.8,+0.1 C:H.>C.Hyt+e- 
















radical. Az; agrees exactly with the value reported by Smith’ for the appearance potential of CH,;+ 
in the mass spectrum of methane. 

Provided there is no activation energy for the formation of a bond between a methyl ion and a 
hydrogen atom, the difference, Az; —A;=4.44+0.2 ev, is equal to the strength of the first carbon- 
hydrogen bond of methane. This value of D(CH;—H) is in much better agreement with the recent 
indirect estimate,’ 4.3s+0.2 ev,® than the estimated uncertainty would lead us to expect. 

The reason that the value reported by Fraser and Jewitt? for the ionization potential of the methyl 
radical, J(CHs)=11.2+0.8 ev, is so much higher than our value probably lies in the inability of 
their molecular beam technique to distinguish between the various thermal dissociation products of 
the lead tetramethyl. As there was no evidence to the contrary, they assumed the predominate 
product of the thermal decomposition to be the methyl radical. We found that under all conditions 
the yield of ethylene far exceeded that of methyl radicals, even at the very low pressure of lead 
tetramethyl we employed, as has been discussed above. Thus their ionization efficiency curve was 
probably a composite one representing the sum of ionization efficiency curves for the various ions 
arising from ionization and dissociation of ethylene, methane and the methyl radical. 

Some rough measurements on the relative abundances of the lead-methyl ions, made with the 
furnace cold, showed the trimethyl lead ion to be about fifty times as abundant as the tetramethyl 
lead ion for 15-volt electrons. For this reason the ionization potential assigned to lead tetramethyl by 
Fraser and Jewitt should probably be associated with the appearance potential of Pb(CHs)3* in 
the Pb(CHs)«4 spectrum. This absence of the parent ion from the mass spectrum has also been noted 
for tertiary butyl chloride.'® 























II. The Ethyl Radical 


The furnace temperature was not as critical for the formation of ethyl radicals as for methyl 
radicals. Curves JJJ and IV of Fig. 4 were taken with 4 amperes through the furnace. From their 
initial breaks we deduce 


ITl. C.H;-C2H;++e«- Arn(C2Hs*) = 8.69+0.1 ev 
IV. C.H,--C.Hygtt+e- Ary(C2H4*) = 10.84+0.1 ev. 















Thus we have 8.69 ev for the vertical ionization potential of the ethyl radical. Equation 7V was 
assigned to the appearance potential of C.H,* on the basis of the agreement of the value, 
Arvy(C2H,4*) = 10.8, ev, with that reported by Kusch, Hustrulid and Tate" for the appearance poten- 
tial of this ion in the spectrum of ethylene. 

Since Fig. 3 indicates that the C;H,* ion far exceeds the C.H;* ion in the spectrum obtained from 
the decomposition products of Pb(C2Hs),, it is likely that Fraser and Jewitt were actually measuring 
the ionization potential of ethylene by their molecular beam technique rather than that of the ethyl 
radical. This is the most reasonable explanation of their value, 10.6+0.8 ev. 

As is the case with lead tetramethyl, the abundance of the parent ion in the Pb(C2Hs5)4 spectrum 












8 L. G. Smith, Phys. Rev. 51, 263 (1937). 
* D. P. Stevenson, J. Chem. Phys. 10, 291 (1942). 
10D, P. Stevenson and J. A. Hipple, J. Am. Chem. Soc. 64, 2766 (1942). 
4 P, Kusch, A. Hustru id, and J. T. Tate, Phys. Rev. 52, 843 (1937). 
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IONIZATION AND DISSOCIATION 








is almost unmeasurably small. Thus the value given by Fraser and Jewitt for the ionization potential 
of Pb(C:Hs).4 should probably be associated with the appearance potential of Pb(C2Hs)5*. 
From data on the appearance potentials of the ethyl ion in the spectra of ethane and n-butane, 



























te we deduce 8.63+0.2 ev’ for Iver.(C2Hs). The direct and indirect determinations are in excellent 
agreement. 
a. The results described in this section, and the conditions under which they were obtained are 
summarized in Table I. 
13* DISCUSSION 
la The indirect method® of determining dissociation energies, D( ), from electron impact data 
on- involves taking differences between appearance potentials of similar processes. Formally we have 
-nt R—R’-Rt++R’/+e-;_ Ai(R*+)=AE,, (1) 
yl R—H-Rt++H+e-;  A2(R+)=AE2, (2) 
. R’—R+2H-(R’—H)+(R-H); AEs, (3) 
0 
ate R’—H=R’+H; D(R’—H) =AE,—AE,;—AE35 Ai:—A2—AEs, (4) 
_ where the AE’s are the minimum energies for the various processes. The equality sign in Eq. (4) 
ad for D(R’—H) holds for two conditions: (1) The equalities 4; =AE, and A2=AEz, or (2) the appear- 
yas ance potentials A, and A» exceed AE; and AE; by the same amount. 
_ The direct method of determining D(__)’s from electron impact data is based on the equations 
‘he (R’—H)—>R’*+H+e;_ As(R’t)=AE;, (S) 
- R’R'*+e; A (R')=F(R), (6) 
y 
Sas (R’—H)=R’+H; D(R’—H)=AE;—J*(R’) SAs5—Ae. (7) 
ed Because process (6) involves the removal of a non-bonding electron from a radical R’ it seems 
reasonable to assume the equality of A. and J*. Thus the direct method should give an upper limit for 
D(R’—H). If R of Eq. (2) and R’ of (5) are similar, e.g., aliphatic radicals, one would expect A» 
and A, to deviate from AE, and AE; by the same amount. Under these conditions the only way 
yl D(R’—H) calculated by Eq. (4) can agree with that calculated by Eq. (7) is for equality signs to 
eir hold for Eqs. (1), (2), and (5) and thus for (4) and (7). The excellent agreement found to obtain 
between the direct and indirect estimates of D(CHs;—H) and D(C:Hs—H), indicated above, offers 
strong evidence for the absence of an activation energy in such reactions as 
R++H-—-RH}?, (8) 
- R++CH;—RCH;"*, (9) 
a, where R* is an aliphatic radical ion. 
- The datum of Cummings and Bleakney” on the appearance potential of CH;* in the methyl 
alcohol provides another check on the direct method of calculation. 
m 
ng CH,+H,O=CH;OH+H AF®=5.65 ev, (10) 
yl CH;,OH—CH;++OH+e- A(CHs*+) =14.0+0.4 ev, (11) 
- H+OH=H;0 AFE®= —5.17 ev, (12) 
CH;-CH;++e A(CH3;*) =10.09+0.1 ev, (13) 







CH,=CH;+H; D(CH;—H)=5.65+14.0—5.17—10.0=4.43+0.4 ev. (14) 


2 C. S. Cummings and W. Bleakney, Phys. Rev. 58, 787 (1940). 

















Hipple’s!* datum on A(CH;*) in the ethane spectrum does not lead to quite as good agreement. 
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2CH,=C:H.+2H AE*=5.18, (15) 
C.H.>CHst+CHs+e- A(CH3*) =14.240.4, (16) 
CH;-CH;++e A(CHs*) = 10.00+0.1, (13) 
CH,=CH;+H 2D(CH;—H) =5.18+14.2—10.0 (17) 


D(CH3;-— H) = 4.65+0.3 ev. 


The deviation of this value of D(CH;—H) from 
that which we accept, namely, 4.49+0.1," is 
almost large enough to be real. If the deviation 
is real, one must assign an activation energy of 
~0.3 ev to the reaction CH;++CH;—-C2Het, 
but no activation energies to the reactions, 
CH;+H-CH,", CH;++OH-—-CH;0OH?, C.H;* 
+CH;—C3Hs", etc. It seems more likely that the 
value of A(CHs;*) in the ethane spectrum is in 
error by 0.2 to 0.3 ev. 

Data acquired by one of the authors on 
butene-1, on the other hand, clearly indicate that 
reactions involving the vinyl and allyl radical 
ions adding to H, CHs3, or C2Hs to form the 
corresponding olefinic ions require an activation 
energy of ~0.4 ev. The addition of a vinyl 
radical to an ethyl ion, however, seems to require 
no activation energy. 

In a recent paper’ it was pointed out that for 
exactly similar processes in various molecules, 
simple ionization for instance, the differences 
between the initial breaks of the ionization 
efficiency curves are equal to the differences 
between the intercepts of the extrapolated linear 


18 “> Hipple, Phys. Rev. 53, 530 (1938). 

4 This value D(CH;— H) =4.49+0.1 seems to be a good 
representation of the indirect value, 4.33 ev (reference 9), 
the direct value from J#(CH3) and A(CH;"*) in the methane 
s rum 4.44, and the value recently deduced by Kis- 
tiakowsky and co-workers (J. Chem. Phys. 10, 305 (1942)) 
from kinetic data, 4.3. ev. Similarly we take D(C2.H;—H) 
=4.29+0.1 from the indirect value, 4.1¢,° the direct value, 
4.24, and the kinetic value 4.22. We have modified the values 

uoted by Kistiakowsky and co-workers to account for the 
yet between that which they took for D(H—Br) 
and that given by Herzberg (reference 6). 

15 It should be noted that on the basis of the conclusions 
reached in the pores paragraphs it is no longer ible 
to reconcile the electron impact values of DiC 3—H) 
with the estimate made by Burton (J. Chem. Phys. 6, 818 
(1938)) from photochemical data. See reference 9. 
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portions of the ionization efficiency curves. To 
state this in another manner, for such processes 
as M@—>Mt+-+e-, the difference Vz~— Vg~ = 1.06, 
+0.060 ev where V7 is the linear intercept and 
Vg~ the initial break of the ionization efficiency 
curve. For eleven reactions of the type R—R’— 
R++R’/+e-, Vi~— Va~=1.54+0.09 ev. If the 
linear intercepts are used, the combination of 
the appearance potential of C.H;* in the ethane 
spectrum with D(C:H;—H) leads to 9.17 ev for 
I#(C2Hs), while the initial break gives 8.63. Both 
methods yield 8.6 ev when applied to the ioniza- 
tion efficiency curve of C2H;* in the ethyl radical 
spectrum. The fact that the value of J*(C2Hs) 
deduced from the initial break of the ionization 
efficiency curve of C.H;* in the ethane spectrum 
agrees with the directly measured value offers 
further substantiation of our conclusion that the 
linear intercepts are not applicable to the deter- 
mination of appearance potentials.'® 

It is interesting to note that the decrease in 
the vertical ionization potential brought about 
by the substitution of CH; for H in the methyl 
radical, 10.0 to 8.6 ev is parallel to the difference 
previously found between methane® and ethane,’ 
ethylene" and propylene,” chlorine atom and 
methyl chloride,'’ and other similar pairs. Price!” 
has discussed the reasons for these decreases. 

In conclusion, we should like to express our 
sincere thanks to Dr. H. A. Beatty and his col- 
leagues at the Ethyl Gasoline Corporation Re- 
search Laboratories for their generous cooperation 
in supplying us with the specially purified samples 
of the lead tetramethyl and lead tetraethy]. 


1942), P. Stevenson and J. A. Hipple, Phys. Rev. 62, 237 
17 W. C. Price, J. Chem. Phys. 4, 539 (1936). 
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The general theory of the magnetron tubes, as developed in a former paper, showed that 


the magnetron should sustain oscillations in an outer circuit, when the frequency was nearly 
v2 times the Larmor’s frequency wy. The present paper contains a more detailed study of 
the behavior of a magnetron with one cylindrical anode, when the radius a of the filament 
and 6 of the anode are both taken into account. It is shown that the magnetron is able to 


sustain oscillations on the frequency 


8 w = wy (2+ 2a*/b*)) 


which lies between V2wy and 2wy, according to the dimensions of the electrodes. The agreement 
of this theory with some interesting results obtained by Blewett and Ramo is shown, and 


the limiting case of the plane magnetron is briefly discussed. 





1. INTRODUCTION. PROPER OSCILLATIONS OF AN 
ELECTRON CLOUD IN A MAGNETIC FIELD 


HE aim of the present paper is to give a 

more precise and accurate discussion of 
some results obtained in a former paper! and to 
show their agreement with a similar theory 
developed by Blewett and Ramo? and with the 
very interesting experimental proofs brought 
forth by these authors. 

All the known effects, as the apparent change 
in dielectric power, and the possibility of 
sustaining oscillations depend on the fact that 
the rotating electron cloud has a proper fre- 
quency of oscillations (provided the radius of 
the filament can be neglected) given by 


w=V2wy, (1) 


where wy = Larmor frequency, as noticed in the 
previous paper. It can be shown by a very 
direct calculation how such oscillations of 
cylindrical symmetry will take place in the 
electron cloud. Equations have been obtained 
[see Eq. (12) of reference 1] for motions of this 


type 
n=C/r, #/r=(e/m)2xpt+7?—wx*, (2) 
where 7 is the angular velocity with respect to 


an axis rotating with Larmor angular velocity. 
Critical conditions, as observed in a magnetron, 


1 L. Brillouin, ‘‘Theory of the magnetron. I,” Phys. Rev. 
60, 385 (1941) and Elec. Commun. 20, 112 (1941). 
adie Blewett and S. Ramo, J. App. Phys. 12, 856 
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correspond to 
C=0, 7=0, po=(m/2re)wy’. (3) 


Let us now assume small oscillations of cylindri- 
cal symmetry to be excited in this cloud. An 
electron, originally located, at the instant /, at a 
distance ro will be displaced to 


r=ro(1+&) in which € is small. (4) 


This means an expansion of the whole cloud, 
the density of which becomes 


p=po(1 — 2£) (5) 


so that the total electric charge pr? remains 
constant. Introducing (4) and (5) in the equation 
of motion (2) one obtains 


#/r=E/(1+€) =wn*(1 — 28) —wn? 
or 
t= — 2wné, (6) 


which yields very directly the proper frequency 
(1) of these cylindrical oscillations. It is now 
easy to understand how a medium possessing 
such a proper frequency of oscillations will show 
an abnormal behavior of its dielectric power, 


€= 1 — (2wy?/w?) e.s.u. (7) 


as was proved both theoretically and experi- 
mentally by Blewett and Ramo. 

Turning back to the problem of the magnetron 
and assuming these elementary calculations to 
yield a first approximation we come to the 
following conclusion: A magnetron, with a 
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filament of radius a and an anode of radius 6 
should behave like a cylindrical condenser filled 
with a medium of dielectric power ¢« [Eq. (7)]. 
Its capacity would then amount to 


C=—_—_——_ 8 
2 log (b/a) @) 


per unit of length of the filament. This yields 
an impedance, at the frequency w 
1 2 log (b/a) 2w log (b/a) 


Z=-—= 


Cw we 


(9) 





2wH? — w? 

Z would be positive for w<V2wy and negative 
for w>vV2wy, so that the magnetron would 
behave like a self-inductance below its frequency 
of resonance V2wy, and like a capacity above 
this proper frequency. 

Such a treatment is, however, just a crude 
first approximation; a more accurate discussion 
shows that the internal frequencies of oscillation 
are not constant throughout the electron cloud, 
but vary from 2wy near the filament to V2wy at 
a distance. This will result in a modification of 
formula (9) as shall be shown in the next sec- 
tions, where the problem of the positive or 
negative internal resistance will be discussed 
simultaneously. 


2. INTERNAL IMPEDANCE OF A MAGNETRON 


The theory developed in ‘‘Magnetron. I”’ for 
a one-anode magnetron is correct so far as 
Eq. (57) is concerned, but some approximations 
made in the deduction of Eq. (58) proved to be 
unreliable, and the calculations have to start 
again from that point. 

It should be first noticed that all potentials 
V(r) are with reference to the filament, taken as 
zero potential point. The current J per unit of 
length of the filament has been defined by a 
formula (38) and is thus taken as positive if 
running from the filament to the anode. This is 
the opposite of the usual assumption in elec- 
tricity, where a current running from a point at 
potential zero to a point at a positive potential 
V(b) should be counted as a negative current.’ 
We shall now revert to these usual definitions 


* The same correction was also applied in ‘‘Magnetron. 
II,” Phys. Rev. 62, 166 (1942), p. 175, before Eq. (77). 
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and change the sign of the currents, writing 


Be = Je! 
t 
I dt =i(J/w)e'{ 1 —e-#7] 


to 


t=to+r (10) 


instead of formula (54) in ‘“‘Magnetron. I.” As 
already noticed in Section 7 [Eq. (60)] of the 
earlier paper it is necessary, for physical and 
mathematical reasons, to add in Eq. (49) a small 
damping term s. This will result in the following 
small change ih formula (57): 


OVe m 2iw 1—e-“" 
=—w’r,= Je - ° (11) 
Ore w*(r)—w*+isw 








This formula gives the alternative potential 
V.(r) at a distance r; it contains the transit 
time 7 from the filament (r =a) tor. The damping 
coefficient s is assumed to be very small, and its 
absolute value will play a very small role 

The assumption is that the direct current J, 
flowing through the magnetron is small enough 
not to perturb appreciably the potential distribu- 
tion V(r) inside the tube. The average potential 
should thus be the one obtained under critical 
conditions 


m a*\* 
V(r) = ~~ anr(1 ~~) (12) 
2e r? 


as shown in “Magnetron. I,” Eq. (51), while 
the space charge density is 


Mwy a‘ 
alr)= ( 1+) (13) 
2re r! 


and the proper frequency, Eq. (52), 





a‘ 
v'(7)=2an'(14+—), (14) 
r 


The value of w(r) thus decreases from 2wy near 
the filament down to V2wy at a distance. We 
furthermore need the relation between transit 
time 7 and distance r, which we obtain from 
Eq. (43), corrected for the change in the sign of 
the current: 


2 dV, 


. m a‘ 
a —"anr(1-). 
r or e r! 








SS ee ..9> 


=_— Sa of 
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Fic. 1. Curve described by 8 in the imaginary plane, 
when w increases. 


MOH 
=e A(t -)=— _(1 -“), (15) 
21 .e 2waL? 


where L is the characteristic length defined in 
“Magnetron. I,’’ Eq. (33). 


Hence 


[L?= —el./mwy'*. 


This enables us to give a more precise statement 
of our assumptions: It is assumed that the 
direct current J, is small enough to make L 
smaller than the radius of the filament: 


L<a, =wa?/2wyLl*>1. (16) 


Expressions (14) and (15) for w? and r must be 
used in formula (11) from which we obtain the 
potential V,(7): 


V.(r) = 2iwJ .e'** 


S( 1—e-* \- 
xX — 
a \2wq?—w*?+iswt+2wy?(at/r*)] 4 





(17) 


with 


(r) =wr(r) = «(—-*). 


The problem is now to perform the integration, 
and to obtain the potential V,(b) on the anode, 
from which we deduce the internal impedance 


of the magnetron: 
V. a(b) _ 
2iw fon (-- 


18 
“a (18) 
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It is convenient to use a new variable ¢ and to put 
=r2/a? 
—2wy? w*? — 2wy?+isw 


B= = 2wy* » (19) 


2wu? —w*+isw (2w? — w?)? + $2? 








which yields 


tw 
Z= 





2wy? —w*? +isw 





fd. 


oe —ito(t—1 
x f exp [—ifo(¢—1/5)] (20) 
is 


1 ‘ae 


This is the integral to discuss, and it should be 
stated that no approximations have been used, 
but for (15). These equations differ from those 
of ‘‘Magnetron. I,” only for the introduction 
of the damping term s. The physical necessity 
for this damping term was already emphasized 
in “Magnetron. I,’”’ but one more word of 
explanation should be added: There is no 
damping term for a magnetron of infinite length, 
which does not emit any radiation (no more 
than an infinite solenoid or a toroidal one) but 
for a magnetron of finite length, damping is 
unavoidable, as any motion of the electrons 
inside the magnetron will result in electro- 
magnetic fields emitted from both ends, and 
radiative dissipation of energy. There is also a 
possibility for energy dissipation inside the 
magnetron itself resulting from the non-linear 
terms which have been neglected in the equations. 


3. INTEGRATION AND DISCUSSION 


It is essential to discuss the variation of 8 
and # as functions of the frequency w. The last 
expression must be of importance in the calcu- 
lation of the integral (20), which has two poles 


fp=+6}, fp*= (21) 


This condition has a clear physical meaning, as 
it corresponds to the resonance on one of the 
internal frequencies of the electron cloud. If we 
neglect, for one moment, the damping term s 
(supposed to be very small) then condition (21) 


means 
w? = w?(r) = 2wy?(1+a*/r*) 


which indicates resonance on the layer at dis- 


(21a) 














L. 





Sp had 
“J 


Xp 





o| i-7-- BE 
a 


Fic. 2. Curve described by the poles P for increasing an. 
tp=fhaar+iyp. 





tance r. Here is the place where the approxima- 
tions made in Section 1 prove unreliable: The 
radius a of the filament may be extremely small, 
but nevertheless, the inner proper frequencies 
always range from V2wy (r large) to 2wy (r=a). 

Figure 1 shows the motion of the point 8 in 
the complex plane, when w varies from zero to 
infinity. According to Eq. (19) 


1 
w)?— 1 —1i(w1/Q) 
w= w/V2wy = y/V?2, Q=V2wu/s. 


The curve starts froom B= —1 for w,;=0 with a 
vertical tangent, and then runs upwards, with 
an almost circular shape. The upper point 8=7Q 
corresponds to w:=1 while w:=~ brings the 
representative point near the origin. The dotted 
curve corresponds to negative w; and has no 
meaning in our problem. The damping term s 
is supposed to be very small, so that the Q 
factor will be extremely large. Figure 2 shows 
how the position of the poles 6! varies with the 
frequency w. For very low frequency the poles 
are in +7, then describe large loops, reach a very 
large distance Q! in the 45° direction (for w:=1), 
and come back to 0 along the real axis, for w: 
infinite. 

The integral (20) is taken along the real axis, 
from ¢=1 to }?/a? and this path runs below the 
trajectory followed by the pole. The pole may, 
however, be very near to it (just above it) when 
its abscissa x lies between 1 and 6?/a? which 





B= 


with 
(21b) 
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means approximately 
2wH 2 w 2 wa(2 + 2a*/b*) i, 


In case } is much greater than a the lower limit 


is very nearly V2wz. 
The integral (20) may conveniently be split 


into two parts 


f=b?/a? 
f + -df=Nt+l, 
f=1 


b?/a? 2¢de 
n=3f " 


1 —_ 


(22) 











=}|log (¢?—8) |”, 
23) 





M/a* exp [—tto(S — 1/5) ] 
In= -f ¢d¢. 
1 ao 
logarithm of imaginary 


Integral J, is the 


expressions 
21, =log [(b*/a*) —8 ]—log (1—8). 


Now, as well known, for an imaginary g= ¢,+7¢; 


(24) 


log g=log | ¢| +i arctg ¢i/¢,, | ¢|?=¢,?+¢7. 


Both terms in (24) are of the following type, 
with X real (1 or b*/a*) and 8=8,+78; asin (21b), 





log (X —8) =log |X —8| +7 arctg 


=log |X —B| +iy. 
The y angles are indicated on Fig. 1, hence, 


b*/a*—B 





27,=log +16, 6=yYo-—y1. 








As noticed before, the damping coefficient s is 
very small, which means that Q is a very large 
quantity, and we assume 


Q>b*/at; B~B,, 


Figure 1 shows then that the angle @ is very 
small, except when the point P comes very 
near the integration path (conditions 22), in 


Bi =0, (24a) 





or 
RE 








(22) 
mit 


split 


(23) 


- 1p. 


s is 
urge 


4a) 
ery 


ery 
in 








which case @ is practically equal to x. Finally 


b*/a* —B, 


21, =log +i[ xr]; 








(24b) 


r 


[x ] for conditions (22) only. 


Turning now to the J; integral, in (23), we notice 
that the coefficient & is large, so that the 
exponential oscillates very rapidly, which will 
practically result in cancelling the integration 
wherever the ¢/(¢*?—8) factor is a slowly varying 
function of ¢. Hence the most important part 
of the integral comes from the neighborhood of 
the poles {p= +6}. As shown on Fig. 2, the pole 
with the plus sign may come near the integration 
path, which runs from 1 to 6*/a? along the real 
axis. 

Two cases must be distinguished: 
A—pole near the integration path, but outside 

of it. 


B=6,+78; and 6,<1 or £6,2b*/a‘. 


B—pole near the integration path, and inside of 
it, which corresponds to conditions (22). 


1 <6, <b*/a‘. 


Case B needs a special discussion before it can 
be reduced to a similar situation as in case A. 
We shall, in case B, divide the integration into 
3 parts: 


4 1<¢°<B,—a_ B; negligible, 


Il. B--a€f*<8,+a exponential prac- 
tically constant, (25) 


Ill. B,+acKg*<Cb*/a* B; negligible. 


a is a small quantity, which is chosen such as to 
satisfy the two following conditions: 


2B, 
B<a<—_—__—__. (25a) 


The first inequality 8;<a enables one to neglect 
8; in both intervals I and III where one has 


f°—B=f°*—B8,—18; = f?—8, 
| ¢?—£,| 2 a>Bi. 


The second inequality results in the fact that 
the exponential is practically constant inside the 


(26) 
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second interval where 
f=VB(1+e), 1/f5=(1—6)/V68,, 
£? = B,(1+2e), 
fol — (1/5) ]= bo VB-— (1/8) ] 
+ foel /8-+(1/ /B,) ]. 


The exponential is practically constant and 
the ¢ term can be dropped if 


efoL. V8.+ (1/V/B,) <1 ’ 


but the interval II is defined by 


(27) 


lg?—B,|<a@ or 28,|e| <a, 
and the inequality on ¢ is certainly satisfied if 


fol V8 +(1/V/B,) \(a/28,)<1, 


which is the second condition stated in (25a). 
These conditions (25a) can be satisfied only if 


2 
Fw (28) 
VB, €ol1+(1/8,) ] 





and this will result in a condition to be fulfilled 
by the damping coefficient s or the Q factor. 
According to (21b) and (16) 





w;?—1 1 
(wi:?— 1)?+w,?/Q? w?>— 1 
w/Q @1 


(w:2—1)?+012/Q®  Q(w:?—1)' 
in which 1/Q? is neglected. 





w a? wa? 





{=— —=-—_>1 
Qwn L? VIL? 


Hence condition (28) yields 
w) 4a? 


> > 
(wi?—1)!2v2L? 





(29) 


When condition B (that the pole be inside the 
interval of integration) is satisfied, the reduced 
frequency w; lies between 1 and v2. On the 
other hand, & is large, which means L?<a?, and 
condition (29) proves that Q must be very large 
in order to make our approximations reasonable. 
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This was always our general assumption. At any rate, conditions (24b) and (29) both postulate a 
very large Q value and extremely small damping coefficient s, which corresponds to actual condi- 


tions in magnetrons. 
As a result of this discussion, we are left with the following integrals to calculate (case B, condi- 


tions 22): 


(B,—a) 4 b?/a? 
= J -f  expl-iter- VN 
+a)4 
" (ta) ede 


—exp (—itelve.—1/VB)1f : pono (30) 
a) r tpi 


The first two integrals correspond to cases I and III and the last one to case II. Let us first notice 
that we can add to the first two integrals an additional one, taken from (8,—a)! to (8,+a)! (on 
the interval II) and which is practically zero. 








(6,+a)4 tdt (+04 ode 
f exp [—ito(¢—1/t)]_——~exp (—ifoL /Be— (1/8) ) J 


(8,—a)4 f?—£, ,—a) eH B, 


=} exp (—tfol V8,— (1/V/8,) ]) log |a/a). 


According to (27), the exponential is constant in this interval, hence the integral reduces to a loga- 
rithm, which is now taken on real quantities and yields log |a/a| =log 1=0. Integral (30) can thus 


be written 
(8, +a) 


rds 
exp (itl v8.-(/VaID J | = GN 
a) vr tpi 





b?/a? tdt 
a f exp (ie — (1/0) 
1 2 


The last integral is again a logarithm, but contains an imaginary part 8;; it is similar to J, [Eq. 
(23) ] but with different limits and yields 


(B,+a)4 ¢dt 
f ed be | 
(s,-a)§ £?—B,—1B; 


a—1B; 
a+1B; 











1 
2 a 
The angle @,, as shown on Fig. 3, is given by 

tg(0a/2)=a/Bi>1, 30.~}n, 


according to the inequalities (25a), while the first log is zero. Summarizing these results we may write 





b*/a* ¢d¢ in 
h=-f exp (—itsl- 0/9 )——-| = exp (- seve (1/v8dD) 2) 
1 


r 


where the last term inside the brackets has to be used only when conditions (22) (case B) are fulfilled. 
We are left now with the problem of estimating the contribution of the first integral, where all 
the difficulties related with the imaginary part in the denominator have been excluded. This first 
integral now keeps the same type in both cases A and B. © 
As already noticed, the exponential oscillates very rapidly, and this practically cancels most of 
the integral, but for the neighborhood of the pole {p=8,4. Near this ) point we may conveniently 


use the following expansions: 


f=V/6,(1+¢), f-1/f= V8-(1+¢) -(1 —€)/V/Br, 





es 
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* de 
I,= —} exp (— itl ve (1/V8)) f exp (—d£oel./Be+(1/+/8,) ])— 
€ € 


"2 


dn 
= exp (—iteV8—(1/Va)) f 
with ioe foe /8,+(1//8-) ]= Sol 1+ (1/8,) ](¢ — VB,). . . 


The new integral in 7 is 


» dn f* ., an koe 
f cnn f (cos y—1 sin »)— = — Cini t+Cine+i siny—i sing 
AT n n n 


according to the usual definitions‘ of integral sines and cosines for positive x 


; 7 dt 7 dt 
Cix = -f cost—, six= -f sin t—. (35) 
=z t z t 


When x is negative, the same definition may be used for six as the point x=0 is no pole for this 
integral, and the relation holds 
si(—x) = —x—six. (36) 
For Ci(—x), our integration path runs through the pole, which means that we must take (see refer- 
ence 4, p. 4) 
Ci(—x) = Ci(x). (37) 


Here we shall notice that the additional — term in (36) comes into play only when the lower 
limit 4; is negative, which means that conditions (22) (case B) are realized. But then we have also 
an additional term in (32), and these two terms exactly cancel each other, hence, we obtain the 
following formula, as a good approximation for all cases: 


I,= —} exp (ito /8-— (1/+/8,) ]) — Ci} m| +Ci| n2| isi | m1| Fisi| n2| }. (38) 
the lower signs, for the si functions, correspond to the case of negative arguments 7 or np. 
As well known, the Ci and st functions become rapidly very small for large values of the argument 7 


sing cos 7 
n>1, Cin= , sine . (39) 
n n 


These asymptotic expressions can be used as soon as 7 is larger than 2x. Now we must keep in 
mind that the approximate expression (38) for Jz has been obtained only for small values of ¢«. But 
a small ¢ may mean a rather large n, according to (33), because £ is assumed to be large. Furthermore, 
we must notice that both J; and formula (38) become very small for large » and so we may, for 
practical purposes, assume formula (38) to hold for any value of n, even large. The limits are 


fi=1, m= fol 1+(1/8,) ](1— V8), 
f2=6?/a?, n2= fol 1+(1/8,) ][(b?/a*) — VB]. 
As soon as 9; and m2 are none too small (larger than 10), the asymptotic expansions (39) may be used 
Ci| n| Fist] | ~(1/|9|)(sin || +7 cos ||) =(1/)(sin 1 +4 cos ) = (é/n)e*. (41) 


Hence, the following result may be used only when 7; or 2 are not too small. 


tiff (vant) -2ol-a(vi-Z)-a] 
g=—)| — CE —1é0 oo $0, i —— —1éo p= ==§ ° 
a. Val ole 179 ae | 


‘E. Jahnke and F. Emde, 7ables of Functions (Teubner, Berlin, 1938), third edition, p. 3. 


(40) 
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But, according to the approximations used in (33), 


fol V/Br— (1// Br) ]+ 02 = Eo f2— (1/f2) J = Sol (b?/a*) — (a?/b*) J, 
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fol VBr— (1// Br) J m © Sol $1 — (1/51) J=0. i. 
Hence, 
1 [ 1 exp [ —£o { (b?/a?) — (a*/b?) } J 
I» =~ = | (44) 
2eo[1+(1/8,) L1—-V/B, (62/a*)—V/8, 
Returning now to Z from Eq. (20) and using (24b) and (44) we finally obtain 
3w b*/a* —B, 1 
Zz i log | —- 
2wy? — w*-+1Sw 1—8, foL1+(1/8,) J 
exp | —ifo[ (b?/a*) — (a*/b*) J} 
x{—- -CeI]}. as) 
1—/8, (b?/a?) —V/B, 


The —[] term comes in only when conditions (22) are fulfilled. One should be reminded that 
formula (45) is valid only for large n:n2 which means 4/8, is very different from 1 or b?/a?. Let us 


here notice the following: 


1. s may now be neglected, as very small, its introduction was only necessary for the calculation 
of residues near the pole. Q is practically infinite, hence, 


2. B=B,= —2wx?/(2wn?—w’), 
3. fol (b?/a*) — (a?/b*) ]= rar. 


(46) 


Tab represents the transit time for electrons running from the cathode a to the anode 6. 


4. PRACTICAL RESULTS 


Let us first assume the magnetron to be 
operated below its critical conditions, with no 
current flowing to the anode, which makes £ 
infinite and J; naught. Formula (45) then 
reduces to 


1 w 
Z= | i log 
2 2wy? —w? 
the [x] term comes in only when conditions 
(22) are fulfilled, namely: 


wa[ 2+ (2a1/b*) Jig we 2wy. (22) 


This proves that such a magnetron behaves like 
a pure inductance for low or high frequencies, 
and like an inductance plus a positive resistance 
for frequencies defined in (22). The following 
cases may be considered: 


wXVlwy, — ~<B<—1. (48a) 


For such low frequencies the [2] term does not 
come in and 


1 tw — 
Z~-————- lo 
, 1+ |B! 


2 2wy*?—w 
w>2wy, 0 <B<1. 


(b*/a*) —B 


| -Cr}]. (47) 

















For very high frequencies the [aw] term also 
drops and 


1 tw b4 
Z=~-— ———— log (-). (49b) 
2 2wy?—w? a‘ 


This later case corresponds to the result pre- 
dicted in Eq. (9) and checks the conclusions of 
Blewett and Ramo as explained in Section 1. 

The formula (47) needs also a special discus- 
sion of the behavior near V2wy, when the 
denominator is zero, as shown below. Let us 
assume w?=2wy*+e, 8B = 2wy?/e, € small. 


(b4/a4)—B 1—(€b*/2wy*a*) € b! 
=— =1+ ( 1— -). 
1-— (€/2wy?) 2wy" a‘ 

















1—8 
Hence, 
(b*/a*) —B € [ b4 
log = 1~ (— | 
1-—,s Quonl a‘ 
and 











t w € b4 1 b4 
Z=- (:-—)- (-.-1) (50) 
2 (—e) 2wy? a‘ 2vV2wH a‘ 


which remains finite. This shows that these 
conditions do not lead to any difficulty. 
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The resistance term in [2], which comes in 
when conditions (22) are fulfilled is positive, as 
the minus sign before x is compensated by the 
negative denominator 2wz?—w*. It should be 
explained by the following physical considera- 
tions, based on the remarks made at the end of 
Section 2: When the frequency lies between the 
limits (22), a certain layer of the electronic 
space charge oscillates very strongly on resonance 
conditions and radiates energy outside of the 
magnetron through both ends of the magnetron. 

Let us now assume that the magnetron is 
operated just above critical conditions, with a 
small current flowing to the anode and a large 
but finite & value. The complete expressions 
from (38) must now be taken into consideration 
but for practical discussions we may use the 
simplified formula (45). The internal impedance 
of the magnetron will yield a real part R or 
resistance, with the following approximate 


values: 


b2 
Case I. ove dan <e<onf2+(2° “)] 

















8 very large, 
w cos £o[ (b?/a*) — (a?/b?*) ] (1) 
~ Yon? —w? 2EoL1 + (1/8) IL(b?/a*) —/8] 
Case II. 1<1/8<6?/a? (conditions 22), 
w 7 1 
x |. + 
w?—2wy2t2 2&[1+(1/8) ] 
1 ol (b?/a*) — (a?/b?) 
x cos £o[ (b?/a*) — (a? =". (52) 
1-/e —_(B*/a®)-/B I" 
Case III. VW/8<1, w>2wy, 
w 1 
(53) 





~ @— en? 26o[1+(1/8) 11 — 8) 


Let us now discuss the sign of the resistance: 
In cases II and III the resistance is always 
positive. In case II, Eq. (52), & is large and 1/8 
lies between 1 and 6?/a*® so that the first term 
in 2/2, inside the brackets, is by far the most 
important. In case III, Eq. (53) the result is 
obviously positive. Hence the discussion is now 
centered on case I, which shows a certain 
number of frequency bands with negative 
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resistance, according to the sign of the cosine. 
The denominator is positive, and we obtain 


R<0O when cos £9(b?/a?—a?/b?) <0. (54) 


The negative resistance may reach high values 
when the frequency is such as to make the 
denominator D very small. By use of (46.2) 
we have 


D=2£0(2wn*—w*)[1+ (1/8) JL (6?/a*) —--V8] 


- vom (5)-(S2a) | 


This amounts to 
D= fo(w?/wu*)(w?— 2wn*)! 
X [V2wn — (b?/a*)(w?— 2wy*)*] (55) 
which is zero for 
w= Vw (561) 
w=wy(2+2a*/b*)!. (S611) 


or 


The first case does not seem of real importance, 
first because the corresponding zero would be 
avoided by keeping the damping term s as in 
Eq. (45) and second because the calculations 
and approximations made in Eq. .(33) and 
following are no more reliable, as for such 
frequencies the pole runs far away from the 
integration path, on the line at 45° (Fig. 2). 
Hence it does not seem that the negative 
resistance might reach any high value near the 
frequency (56l). 

The second root (56I1) is much more important 
as it means that the frequency w is just the one 
which sets into resonance the space charge 
density near the anode. The conclusion of all 
this discussion is the following: 

One anode magnetrons should yield a number of 
very narrow frequency bands with negative re- 
sistance [Eq. (54)] for frequencies lying just 
below the limit (5611). 

This checks the general conclusion of a former 
paper (Magnetron. I). One should notice here 
that the negative resistance R does not become 
infinite on the limit (56II), as Eq. (51) seems 
to indicate. Near this limit, the approximate 
formula (45) should not be used any more, as 
it rests on the use of asymptotic expansions (39). 
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One should revert to the original formula (38) 
with integral sines and cosines, thus avoiding 
infinite R values. 


5. PLANE MAGNETRON 


As an example of application of the general 
theory, the case of the plane magnetron will be 
briefly discussed. This case is obtained as the 
limit of a cylindrical magnetron when the radius 
a of the filament is increased to infinity, while 
the distance d= 6—a between filament and anode 
is kept constant. Hence 


a2, b=a+d=a(1+p~). (57) 


A point at a distance r from the center will be 
defined by a parameter p: 


r=a(i+p) (58) 


in which p=(r—a)/a is small. The angular 
velocity of the electrons is given by Eq. (15) of 
“Magnetron. I’’ which now reads 


6=wn[1 — (a?/r*) ]=an(1—1/(1+)?) 
=wa(2p—3p?---). (59) 
This means that the tangential velocity becomes 
£=76=2awyp(1+p)(1—$p) + 2wyxap(1—4$p). 


Denoting y the distance r—a from the surface of 
the cathode, we obtain 


£=2wyy[1 — (y/2a) ]. (60) 


For the plane magnetron (a—~) the x velocity 
increases proportionately to the distance y from 
the cathode. 

Under critical conditions (no current on the 
anode), the space charge density is given by 
Eq. (25) of ‘Magnetron. I’’: 


mw? a‘ Mw? 
ei (1+=)~ (1-29), (61) 
rT 


2xe we 








which on the limit of a plane magnetron (p=0), 
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is a constant. In a plane magnetron (of infinite 
length) operated under critical conditions, there 
is a uniform space charge density (61) between 
the plane cathode and the plane anode. These 
charges have no perpendicular component of 
velocity (y=0) and move parallel to the elec- 
trodes with a velocity « given by Eq. (60). 
The voltage V is obtained from Eq. (23) of 
“Magnetron. I.” 


Mw 17” a? 2 
SS P-(5)] 


Mun? 


2 2mwx? 
=~ at(1+2-——) =F arprcr 9 
é 


2e 
2mwy? y 
Vol) = »( 1-"), (62) 
é a 


This yields the anode voltage V(b) when y=d is 
inserted in the formula. When the total electric 
charge Q bétween the electrodes is computed, 
it comes out to be 


Q=—VoL/2xd, (63) 


L is the length measured along the x axis, and 
the calculation refers to a unit length of the 
filament. The plane magnetron thus yields a 
capacity twice as large as a plane condenser of 
same dimensions. 

Such a plane magnetron can be used for 
sustaining oscillations, and the discussion of the 
preceding section shows that it should oscillate 
in the neighborhood of the frequency 











or 





b 
w=2wy, ——I, (64) 
a 


which corresponds to the proper frequency of 
vibration for the electronic layers near the 
filament, as emphasized in Section 1. 

It is a pleasant duty for me to thank Pro- 
fessor Le Corbeiller for many helpful criticisms 
and suggestions. 
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ROMPT publication of brief reports of important dis- 

coveries in physics may be secured by addressing them 
io this department. The closing date for this department is the 
third of the month. Because of the late closing for the sec- 
tion no proof can be shown to authors. The Board of Editors 
does not hold itself responsible for the opinions expressed by 
the correspondents. Communications should not in general ex- 
ceed 600 words in length. 





Thermal Equilibrium Between 
Elementary Particles 


GLEB WATAGHIN 
Department of Physics, University of Sao Paulo, Brasil 
January 2, 1943 
ECENT investigations by C. F. Weisziacker,' S. 
Chandrasekhar,’ and L. R. Henrich on the prestellar 
stage of the universe, based on the actual distribution of 
chemical elements and on the theory of thermal reactions 
between nuclei and elementary particles, induce me to 
publish further remarks on an earlier study® of thermal 
equilibrium between elementary particles at extremely high 
temperatures. 

The chief feature of the statistical equilibrium at tem- 
peratures T>8X 10°°C (xT >2mc*, where m is the mass of 
the electron) is the production of pairs by photons having 
energy hy >2mc*. In these conditions (as was stated in the 
above-mentioned papers) the energy of the rest mass mc* 
takes part in the exchanges of thermal energy and, by 
increasing the temperature, the total inertial or gravita- 
tional mass of the created pairs increases also propor- 
tionally to T*. An exact calculation of the constants of 
degeneracy which appear in these new statistical conditions 
leads us to the following conclusions: At temperatures 
T>2mc*/x, the concentrations m* and n~ of positrons and 
electrons and N of photons with hy>2me* are nearly equal 
(as in the energetic cascade showers) 


nt~n-~N~(168/Ch) (xT)? 


The constant of degeneracy a which appears in the for- 
mulae: 


ia 8x © (x? — x9?) bxdx 
n= pS, ateq 
_ 8x © (x2 — x92) 'xdx 
iar ae aro 


where xo= mc?/xT, x = E/xT, and E is the total relativistic 
energy component, is given approximately by: n+/n-e~** 
or, in the case aX1, a~(n-—n*)/(n-+n*). 

In the case of thermal equilibrium between protons, 
neutrons, electrons, and photons, the general formulae 
obtained for a constant number of heavy particles lead to 
the result that at temperatures =(1/«x)2mc?, at which aS1 
and n*2=n™~ is greater than the concentration of heavy 
particles, the number of protons is nearly equal to the 
number of*neutrons. The statistical laws above-mentioned 
suffer an alteration at temperatures T=6X10"° if the 
definition of distinguishable and independent quantum 
states is altered in a manner proposed by the author,‘ 


VOLUME 63, NUMBERS 3 AND 4 


137 





FEBRUARY 1 AND 15, 1943 






so that the asymptotic number of states in an interval 
(p, p+dp) is (8xV/h*)g(p)p*dp, where g(p) is a “‘cut off 
factor.” 

An entirely new situation appears in the conditions in 
which the total energy of the assembly is so high that the 
average kinetic energy EF of a proton or neutron is 
E>>uc’~108 ev, where yu is the rest mass of the meson. In 
these conditions, as seems to follow from the experimental 
evidence on cosmic rays, the high energy primaries give 
rise to the creation of particles by explosion showers, 
producing groups of mesons, electrons, photons, and 
neutrinos, and in each of such collisions a considerable 
amount of primary energy can be transformed into neu- 
trinos (through meson decay or directly). On the basis of 
the theory of these collisions proposed recently by the 
author,5 we must expect that no equilibrium condition is 
possible in these cases and that the idea of a temperature 
xT>>2ye* cannot be introduced in the usual way in an 
assembly of heavy particles. One must expect a rapid 
transformation of the energy of the heavy particles into 
the energy of the created ones. The kinetic energy of the 
proton results divided among a great number of particles 
of average kinetic energy ~2yc* and is simultaneously 
transformed into the energy of the rest mass (uc*) of the 
created particles and into the gravitational energy. 

One of the most important features of the collisions in 
which many particles are simultaneously created (including 
mesons and neutrinos) is the irreversibility, due to the 
vanishing probability for the inverse processes. This 
probability depends obviously on the cross section for a 
nearly simultaneous collision of m particles of energy 
~2yuc*, which decreases rapidly with increasing n. If the 
number of created particles is proportional to the primary 
energy measured in the frame of the center of masses, the 
inverse transition to the high energy states has a proba- 
bility which tends to O with increasing energy EZ. So the 
equilibrium is practically impossible, because the condi- 
tions of detailed balance cannot be fulfilled in an assembly 
of elementary particles having plausible concentrations. 

If the property of creation of groups of particles in the 
high energy collisions is not only limited to the heavy 
particles (as we want to assume), we are led to the con- 
clusion that the value 7.~(1/x)2uc?~6X 10" degrees of 
the temperature indicates a separation between two differ- 
ent types of thermal processes, the usual one well below 
T. and the others in which a statistical equilibrium is 
practically impossible and in which the creation of mesons 
and neutrinos constitutes the dominant feature of the 
collision process. Instead of the neutrino hypothesis one 
can abandon the conservation laws of energy and mo- 
mentum. Then, for instance, the disappearance of a great 
amount of energy in 8-processes and decay processes would 
produce a decrease of the intensity of gravitational forces 
and thus an increase of gravitational energy (in accordance 
with a well-known suggestion of Landau). 


1C. F. Weisziicker, Zeits. f. Physik. 38, 176 (1937); 39, 633 (1938). 
2S. Chandrasekhar, Astrophys. J. 95, 288 (March, 1942). 

= Mag. 17, 910 (May, 1934); Comptes rendus, 909 (March 11, 
4G. Wataghin, Nature 142, 393 (1938). 
5G. Wataghin, Symposium on Cosmic Rays, Rio de Janeiro (August, 

1941), An. Acad. Brasil. Sci. (Dec. 1942). 
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Proceedings of the American Physical Society 
MEETING AT Los ANGELES, CALIFORNIA, DECEMBER 29, 1942 


HE 252nd meeting of the American Physical Society was held in the 

University of California at Los Angeles, on Tuesday, December 29, 1942. 

The attendance at the two sessions varied from 40 to 45, which is greater than 

anticipated, in view of the trying conditions under which the meeting was held. 

Most of those in attendance had luncheon together, Tuesday noon, in the 
faculty dining room of Kerckhoff Hall. 

Of the 18 papers whose abstracts are given below, only one, No. 7, was read 
by title. Professor S. J. Barnett gave a detailed account of his own extensive 
researches in the field of gyromagnetism. This was an invited paper, given at 
the close of the morning session. By previous arrangement, and with the 
consent of the society, Professor Millikan expanded paper No. 11 into a half- 
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hour invited paper. 


1. The Polarization of the Absorption Spectra of Neo- 
dymium Ions in Crystals. RoperT W. KRUEGER, Douglas 
Aircraft Company, AND E. L. Kinsey, University of Cali- 
fornia at Los Angeles.—It has been discovered that many 
of the sharp electronic lines in the absorption spectrum of 
crystals containing neodymium ions are polarized. Photo- 
graphs were taken of the absorption spectra of neodymium 
ions in crystals of NdCl;-6H:O and Nd(BrO;)3-9H;0 cut in 
various directions with respect to the crystallographic axes, 
using polarized light. The polarization effect has been suc- 
cessfully explained by the assumption of a perturbing field 
having the symmetry of the arrangement of the water 
molecules surrounding the neodymium ion. This field 
produces a splitting of each unperturbed energy level into a 
number of component levels equal to the number of 
irreducible representations of the symmetry group con- 
tained in a reducible representation presented by the 2/+1 
unperturbed eigenfunctions corresponding to the unper- 
turbed energy level. When the mixing of unperturbed 
eigenfunctions to produce a perturbed one is simple, 
polarization effects become pronounced. Knowing the 
correct mixing from group theory, one can predict polariza- 
tions, using the correct selection rules which are found to be 
those for forced electric dipole radiation. The splitting and 
polarization effects observed in at least four groups of lines 
in the bromate spectrum substantiate the theory very well. 


2. On the Heat Transfer from Radiating Gases to Water 
Tubes of a Boiler. R. M. LANGER AND G. POTAPENKO, 
California Institute of Technology.—The transfer of heat 
from radiating gases to water tubes of a heat exchanger 
depends upon the arrangement of tubes, because this 
arrangement determines the “radiant path length,” i.e., the 
average thickness of gas layer which a tube faces. The 
theory developed allows one to calculate the radiant path 
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length for any arrangement of tubes. The application of 
the theory to the case of water tube boilers permits a com- 
plete physical analysis of heat exchange in such boilers and 
leads to the following general formula showing the tempera- 
ture of gases at any point along the pipes: 


raf (E2-i)-u(E2-1)+--} 


where x is the distance along the boiler under consideration, 
F is the fraction of the heat content of the entering gases 
which is absorbed in the first unit of length. T(x) and 7(0) 
are the gas temperatures at a point x and at x=0, re- 
spectively. A is the ratio of the convective to the radiative 
heat transfer coefficient. In the part of the boiler where the 
radiative heat transfer predominates, A is small and only 
the first term in the brackets needs to be considered. If AT 
is the temperature rise due to combustion of the fuel gases 
and if the efficiency of boiler of a length x is E(x) then 
T (0) — T(x) = E(x)AT. The quantities F, A, and AT can be 
calculated for any given set of operating conditions and so 
it is possible to calculate the length of a boiler for any 
desired efficiency E, and vice versa. 








3. The Production of Photons Relative to Ionization by 
Collision in a Townsend Gap. RONALD GEBALLE, Uni- 
versity of California.—Measurements were made in a 
Townsend gap to determine the average number of photons 
liberated by an electron for each secondary electron it 
liberates. Gases used included Hg, Nz, A, and air. In Hg, this 
number was found to decrease rapidly with increasing ratio 
of field strength to gas pressure, in the range from 50 to 
150 volts/em/mm, and at pressures of the order of one 
millimeter. Determination of the absolute number of 
photons per electron depends on the unknown efficiency of 
the photoelectric cell used to count the photons. However, 
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this number is estimated as of order unity. In He, the 
effective radiation has an absorption coefficient of 0.55 +0.2 
cm~ at a pressure of one millimeter of mercury. This large 
coefficient indicates that only those photons of rather high 
energy were counted. There is corroborative evidence that 
the brass surface of the photoelectric cell had a high work 
function. Thus only high energy photons would be expected 
to register. Photoelectric currents in the other gases were 
too small to be measured with the available apparatus. It 
thus appears that these gases have absorption coefficients 
for the effective radiation even higher than that found 


in He. 


4. The Magnetic Properties of Iron-Cerium Alloys. 
J. R. Crark,* S. T. Pan, anp A. R. KaurMann, Massa- 
chusetts Institute of Technology—The susceptibility and 
magnetization of iron-cerium alloys containing 1.07 to 37.8 
atomic percent (0.44 to 19.8 weight percent) iron have been 
measured in fields up to 35,000 gauss and at temperatures 
from 77°K to 1200°K. The structure of the alloys has been 
investigated by means of x-ray diffraction patterns recorded 
at 77°K, 195°K, and 300°K. The alloys containing up to 
about 30 atomic percent iron have a ferromagnetic Curie 
point at about 275°K; the alloy containing 37.8 atomic 
percent iron has a ferromagnetic Curie point at about 
475°K. Below 275°K the intercept magnetization increases 
with decreasing temperature and increasing iron content, 
but there is an abrupt change in the magnetic properties 
between 150°K and 100°K. At 150°K the alloys are readily 
saturated; but below 100°K the susceptibility is field- 
dependent in fields up to 35,000 gauss, the field-dependency 
being most pronounced in the higher iron-content alloys. 
The appearance of new lines and a change in the relative 
intensities of the face-centered cubic lines at 77°K indicate 
a structural change in the 1.07 atomic percent iron alloy. 
This transformation is not apparent in the alloys containing 
more iron. The Curie points observed are believed to be 
due to intermediate phases of the Fe-Ce system. The 
relation between the structural and magnetic transforma- 
tions at low temperature requires further investigation. 


* Now at Stanford University. 


5. Theoretical Dependence of Continuous X-Ray In- 
tensity and Polarization upon Several Experimental 
Parameters. PAUL KIRKPATRICK, Stanford University.— 
The Sommerfeld theory of the continuous x-ray spectrum 
from thin or atomic targets has been put by Weinstock into 
a form which makes computation feasible. Encouraged by 
certain agreements between this theory and our experi- 
ments, computations have been conducted of the depend- 
ence of intensity and polarization upon target atomic 
number, radiation frequency, bombardment energy, and 
spatial distribution. The computations are still laborious, 
requiring about six hours of the work of a skilled computer 
per calculated point. They make possible for the first time 
the visualization of the principal features of the spectrum. 
Curves and space models are used to display the functional 
relations disclosed. 


6. Direct and Fluorescence Excitation of the 11; Level 
in Thick Targets of Thorium. Burr G. BurBANK, Stanford 
University. (Introduced by Paul Kirkpatrick.)—X-rays from 
an oil-cooled thorium target were analyzed by a Bragg 
spectrometer, the intensities at the peak of the La line and 
in its continuous background at both sides being measured. 
The line intensity was partly produced by direct cathode 
electron ionization of the thorium atoms and partly by 
ionization produced by absorption of the thorium con- 
tinuous spectrum. To determine the ratio of these com- 
ponents auxiliary measurements were performed in which 
the thorium target was covered by a foil of palladium thick 
enough to prevent direct ionization of the thorium and to 
provide a palladium continuous spectrum sufficiently 
intense to excite a measurable intensity of Th La by 
photoelectric ionization. From measurement of this purely 
indirect line excitation the indirect component with pure 
thorium targets was computed by a procedure taking 
account of the relative continuous background intensities 
in the two experiments. The ratio of the number of La line 
quanta produced by direct electron collision to the number 
produced by fluorescence for normally emergent x-rays 
from a thick thorium target was found to be 2.8 for 40-kv 
electrons. Determinations were made in the range from 40 
to 100 kv. 


7. Inelastic Scattering of Slow Neutrons. Rosert 
WEINSTOCK, Stanford University.—The first-order (single 
phonon emission or absorption) inelastic scattering cross 
section for neutrons by polycrystals has been calculated. 
The Born approximation and the Fermi (é-function) inter- 
action between slow neutrons and bound nuclei are em- 
ployed. In the preliminary work, attention is limited to 
spin-free nuclei, and both isotopic disorder and magnetic 
interaction are neglected. The Debye approximation of a 
constant sound velocity is used. The ratio of the probability 
of a first-order process to zero order (elastic) is found to be 
roughly of the order of neutron mass divided by scattering 
nucleus mass. At temperatures high compared with the 
Debye temperature of the scatterer, however, the inelastic 
cross section may become larger than this ratio indicates. 
Computations carried out for Fe scattering room tempera- 
ture neutrons show a rise from three percent at 151°K to 
27 percent at 1000°K for the ratio of (first order) inelastic 
to elastic scattering cross section. The sum of these cross 
sections, however, varies only slightly over this temperature 
range. Future work will take into account both isotopic 
disorder and magnetic interaction. 


8. Dimensional Characteristics of the Large Scale 
Condensations of Matter in the Universe. F. Zwicxy, 
California Institute of Technology.—lIt is one of the most 
fundamental facts of observation that the large scale dis- 
tribution of matter is very far from uniform. Clusters of 
nebulae, nebulae, and stars the dimensions of which are of 
the order 10** cm, 10” cm, and 10" cm, respectively, repre- 
sent the most conspicuous condensations on a cosmic scale. 
A complete theory of the structures of these condensations 
is lacking because it involves the statistical mechanics of 
cooperative gravitational systems, a discipline which has 
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not yet been developed. General considerations coupled 
with dimensional analysis show that the “structural length”’ 
a for clusters of nebulae and for nebulae must be of the 
order k(&/I'po)* where & is the average kinetic energy per 
unit mass of all the matter in the universe, po its average 
density, I the universal gravitational constant, and k a 
dimensionless constant. For individual clusters of matter 
which are in stationary equilibrium k assumes the value 
(12x)? if po denotes the central density. These conclusions 
are checked through a comparison with the data obtained 
on the three large clusters of nebulae in Coma, in Perseus, 
and in Hydra. 


9. Static Charges Produced on Moving Rubber-Tired 
Vehicles. S. S. MACKEOWN AND VicToR WoukK, California 
Institute of Technology.—The current flowing from an auto- 
mobile whose rear wheels were turning on a dynamometer 
was measured. These results show that the total current 
generated is proportional to speed, independent of load or 
tire pressure. Measurements made when the tires were 
made semi-conducting and with a series resistance suggest 
a simple equivalent electric circuit that is consistent with 
the experimental results. Approximate values of the con- 
stants of this equivalent circuit were obtained. Constant 
currents as high as 4 microamperes were obtained at a 
speed corresponding to 40 miles per hour. These results 
were checked by road tests. The tires are sufficiently con- 
ducting to discharge the automobile in a short time after 
the car is stopped. 


10. Preliminary Energy Distribution Curve of Cosmic 
Rays. L. Katz, R. V. ADAMs, JR., AND W. E. DEEDs, 
California Institute of Technology. (Introduced by C. D. 
Anderson.)—Several test runs have been made with the 
twenty-four inch, vertical cloud chamber at a field strength 
of approximately 5000 gauss. A preliminary energy distri- 
bution curve has been made, on the basis of measurements 
of 107 tracks, with energies ranging from 100 Mev to 
31,000 Mev. Altogether, there were 54 tracks of negative 
particles, 51 of positive particles, and 2 of doubtful sign. 
In the region above 5000 Mev there were 5 negatives, 9 
positives, and 2 doubtful cases. The tentative curve shows 
several secondary maxima, including one at 1700 Mev, 
which may or may not be due to statistical variation. 


11. The Origin of Cosmic Rays. RoBert A. MILLIKAN, 
H. Victor NEHER, AND WILLIAM H. PICKERING, California 
Institute of Technology.—The authors had predicted that 
since their hypothetical silicon-annihilation rays should 
have enough energy (13.2 bev) to get through the earth’s 
magnetic field at the equator in Peru, though not in India, 
there should be found both at sea level and at all altitudes 
in the Americas a very long plateau of uniform, vertical 
cosmic-ray intensities extending north from Peru clear up 
about to the magnetic latitude of Victoria, Mexico (32.8° 
N Mag.). There the strong band due to the annihilation 
of the oxygen atom (computed energy 7.5 billion electron 
volts) should first begin to be able to break vertically 
through the earth’s magnetic field. This prediction and 
others deducible from the hypothesis were experimentally 
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tested in Mexico and the United States in December, 1940, 
with no results as yet found anywhere at variance with the 
atom-annihilation hypothesis. 


12, An Improved Cosmic-Ray Radio Sonde. WILLIAy 
H. PicKerinG, California Institute of Technology.—A de- 
scription will be given of the cosmic-ray radio sonde used 
for the measurements described in the preceding paper 
(No. 11). The chief difference between this and the pre- 
ceding instrument lies in the method used to scale down the 
counts. This is now done at the receiving station. The 
transmitted pulses are made sufficiently short to avoid any 
loss of efficiency at high counting rates. 


13. Mesotron Mass and Heavy Tracks on Mt. Evans, 
Cari E. NIELSEN AND WILSON M. PowELL, University of 
California.—A cloud chamber in a 2500 oersted field was so 
operated that diffuse tracks permitting ionization measure- 
ment by droplet count were photographed. A few slow 
mesotrons with ionization and curvature suitable for mass 
determination were observed in the approximately 7000 
photographs taken at 14,100 ft. and the 1500 at 12,700 ft. 
Preliminary measurement of two of the best tracks indi- 
cates a mass of about 200 electron masses; the mass is 
determinable from each track to within +15 percent, 
There were observed also many heavily ionized tracks 
made by protons or heavier particles—several such appear 
in each hundred photographs. A count is being made of the 
relative frequency of these heavy tracks at the two eleva- 
tions; the decrease with elevation is not conspicuous. 


14. Collision Electrons Accompanying Mesotrons in 
Lead at 10,000 Feet. Wayne E. Hazen, University of 
California. (Introduced by R. B. Brode.)—A cylindrical 
cloud chamber (30 by 30 cm) containing eight 0.6 cm lead 
plates was expanded at random. Two thousand of the 
photographs have been studied for the production of 
knock-on electrons in the lead. In these photographs 4540 
particles traversed at least three plates without multiplica- 
tion but only 4570 traversals were classified as suitable for 
the observation of knock-on electrons. Since 340 electrons 
were observed, of which only 2.5 percent passed through 
another 0.6 cm of lead, the plates were thick enough to 
establish equilibrium. Ten percent of the knock-on electrons 
were noticeably scattered in the gas and therefore had 
energies less than 105 ev. There were two cases of small 
showers accompanying mesotrons out of the lead, one of 
three and the other of four particles. From the above data 
the probability that a collision electron accompanies a 
mesotron in lead is 7.5+0.5 percent. Seren' found 9.7 per- 
cent for the probability that a collision electron with 
E>4X10* ev accompanies a mesotron in lead at sea level. 
The difference is probably a consequence of the higher mean 
energy for mesotrons at sea level. 


1L. Seren, Phys. Rev. 62, 204 (1942). 


15. Optical Constants of Electron Microscopes. L. 
MARTON AND R. G. E. Hutter, Stanford University.— 
In electron microscopical work it has been previously 
assumed that the field of the objective lens which extends 
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beyond the specimen into the region between it and the 
condenser lens is not strong enough to exert an appreciable 
converging effect on the illuminating beam of electrons. 
Recent calculations of Glaser' and Dosse?® show that if the 
object is placed at a distance from the lens where the 
spherical aberration becomes a minimum, the converging 
effect is far from negligible. Taking into account the neces- 
sary corrections for such an extended objective field and 
with a minimum of assumptions, several optical constants 
of the electron microscope are recalculated. For the new 
values thus obtained, the resolution of the instrument is 
calculated and the role of the limiting factors discussed. 


1 Zeits. f. Physik 117, 285 (1941). 
2 Zeits. f. Physik 117, 316 (1941). 


16. The Influence Function for a Cubic Lattice. H. 
BATEMAN, California Institute of Technology.—In 1912 
Born and K4rma4n studied the periodic solutions of a 
differential difference equation analogous to the wave 
equation for a homogeneous isotropic medium. An ex- 
pression is now found for the disturbance produced when 
one member of an elastic lattice system is given an initial 
displacement or velocity. 


17. Theory of Elastic Oscillations in Plates and Shells. 
Paut S. Epstein, California Institute of Technology.— 
The assumptions forming the foundation of the classical 
theory of plates and shells are open to the objection of 
being arbitrary. A new theory meeting this objection is 
offered and carried through for the following cases: (1) 
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with respect to terms of zero order in the thickness of the 
shell (4), for shells of any shape; (2) with respect to terms 
of second order in h, for circularly cylindric shells (includ- 
ing, as a special case, plane parallel plates). The zero-order 
terms turn out to be identical with those resulting from 
the classical theory, while the second-order terms are 
different. The new theory gives a satisfactory account of 
torsional oscillations which cause difficulties in the classical 
theory. 


18. Continuous Radiation in Prismatic X-Ray Spectra. 
PAUL KIRKPATRICK, Stanford University.—In determining 
e/m by x-ray reflection and refraction' Stauss observed in 
his prismatic x-ray spectrograms unexpected refraction 
between the Ka lines of Mo and the undeviated beam. He 
was inclined to attribute this radiation to imperfection of 
the surfaces of his quartz prism. Pardue* while investigat- 
ing the dispersion of x-rays in calcite noticed the same 
effect and attributed it to refracted continuous spectrum 
radiation of wave-length less than that of the molybdenum 
lines. The band possessed a sharp edge on the long wave 
side which he thought might be a result of a deposit of 
tungsten on the molybdenum target. It is shown in the 
present work, using a glass prism, that the band is indeed 
continuous radiation since it responds appropriately to 
voltage and filtration variations. The sharp edge on the 
long wave side is found to be a photographic effect caused 
by the K discontinuity of the silver in the film. 


1H. E. Stauss, Phys. Rev. 36, 1101 (1930). 
2?L. A. Pardue, Phys. Rev. 38, 1808 (1931). 
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